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ABSTRACT 



INCLUSIVE AND EXCLUSIVE COMPTON PROCESSES 
IN QUANTUM CHROMODYNAMICS 

A. Psaker 
Old Dominion University, 2008 
Director: Dr. Anatoly Radyushkin 

In our work, we describe two types of Compton processes. As an example of an inclu- 
sive process, we consider the high-energy photoproduction of massive muon pairs off 
the nucleon. We analyze the process in the framework of the QCD parton model, in 
which the usual parton distributions emerge as a tool to describe the nucleon in terms 
of quark and gluonic degrees of freedom. To study its exclusive version, a new class 
of phenomenological functions is required, namely, generalized parton distributions. 
They can be considered as a generalization of the usual parton distributions mea- 
sured in deeply inelastic Icpton-nuclcon scattering. Generalized parton distributions 
(GPDs) may be observed in hard exclusive reactions such as deeply virtual Compton 
scattering. We develop an extension of this particular process into the weak inter- 
action sector. Wc also investigate a possible application of the GPD formalism to 
wide-angle real Compton scattering. 
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CHAPTER I 
INTRODUCTION 

Strong interaction physics, the study of hadron structure in general, is one of the most 
fascinating frontiers of modern science. The underlying theory is the universally 
accepted nonAbelian gauge field theory called quantum chromodynamics (QCD). 
In principle, QCD embraces all phenomena of hadronic physics. It postulates that 
hadrons (baryons and mesons) are composite objects made up of quarks, and that 
the color interaction between quarks is mediated by gluons as the gauge bosons. The 
main difficulty of this elegant theory lies in the fact that it is formulated in terms 
of colored degrees of freedom (quarks and gluons) while the physical hadrons are 
colorless. How the hadrons are built out of quarks and gluons has yet to be answered, 
and it represents a challenging task. 

Historically, QCD originated from the constituent quark model PQ I21- In the 
1960's, the development of high-energy accelerators made it possible to resolve the 
structure of hadrons for the first time. The studies of the hadron spectroscopy led to 
the concept of quarks as the fundamental building blocks of hadrons, in particular, 
three quarks for baryons and a quark-antiquark pair for mesons. The natural step that 
followed, was the search for quarks in the experiments with the momentum transfer 
large enough to look inside the hadrons. The first series of such high-resolution ex- 
periments was performed at SLAC by probing the proton structure. The process was 
known as deeply inelastic electron-proton scattering (DIS). It is probably the most 
studied QCD process. The data on DIS have strongly supported the parton model, a 
physical picture given by Feynman OIHEl, in which the proton, when observed with 
high spatial resolution, is built out of almost-free and point-like constituents. They 
were called partons. Thus the dynamics of partons described by the QCD parton 
model should have the property of asymptotic freedom. In other words, the coupling 
constant of interaction between partons is vanishing at small distances (or large mo- 
mentum transfers). In addition, the DIS experiments also suggested that charged 
partons should be spin-1/2 particles that were later identified with quarks. Another 
important result from these experiments was the first evidence for the presence of 

This dissertation follows the style of The Physical Review. 
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electrically neutral gluons. It turns out that quarks carry only about a half of the 
total nucleon momentum. 

The fact that the theory of quark dynamics should exhibit the desired property of 
the asymptotic freedom led physicists to consider the nonAbelian gauge field theory, 
originally constructed by Yang and Mills [F, as the best candidate [71 |HJ inj- In 
order to resolve several difficulties in the constituent quark model, a new quantum 
number, color, was introduced to quarks jTnilIIl[I2]- By identifying color as the SU(3) 
symmetry of the nonAbelian gauge theory fSlE]; QCD was finally established. 

The description of all inclusive hard reactions in QCD is possible due to the so- 
called factorization. According to this property, the cross section of a particular QCD 
process splits (factorizes) into the hard (short-distance) and the soft (long-distance) 
parts. Introducing the factorization scale as a point of separation can only be possible 
by assuring the presence of a large invariant, such as the virtuality of the probe or the 
momentum transfer jTHl UHl 113 CHI HSl 1201 1211 122] • The clear separation of scales is of 
the crucial importance. The asymptotically free nature of QCD allows to discuss the 
short- distance interactions by means of perturbation theory. This is, however, not the 
case for the soft (nonperturbative) stages of interactions. They are expressed in terms 
of the hadronic matrix elements of well-defined quark-gluon operators, i.e. the QCD 
operators taken on the light-cone and sandwiched between the hadronic states. These 
matrix elements accumulate information about long-distance dynamics. They emerge 
as a result of the description of the hard hadronic processes using a powerful tool 
called the light-cone operator product expansion (OPE) [231211 ■ Originally applied 
to DIS [23 I2ni I21j; the technique is based on the statement that, in the particular 
kinematical regime, the asymptotic behavior of the relevant scattering amplitudes is 
governed by the singularities on the light-cone which, in the framework of a light-cone 
expansion, can be described in terms of contributions of definite twist (see [21] for 
a recent review). In other words, one performs a systematic expansion in the inverse 
powers of a characteristic momentum scale, and extracts the leading and higher twist 
contributions, the latter containing information about quark-quark and quark-gluon 
correlations inside the hadron. 

The hadronic matrix elements are process-independent nonperturbative objects 
measured with the help of different probes, such as photons and weak interaction 
bosons. For that reason, these objects are parametrized by universal functions, al- 
lowing to relate various light-cone dominated scattering processes to each other. The 
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well-known example of such phenomenological functions is given by the usual parton 
distribution functions (PDFs). They parametrize forward hadronic matrix elements 
and, since these elements are related to inclusive cross sections, PDFs enter the de- 
scription of hard inclusive reactions. Additional information on quark and gluon 
structure of hadrons is provided by measuring the nonforward hadronic matrix ele- 
ments. They correspond to nondiagonal transitions between the hadronic states in 
momentum, flavor and spin spaces, and hence are described in terms of more general 
phenomenological functions known as the generalized parton distributions (GPDs) 
[211 EOl inn 1^21 IM] • These distributions can be accessed in hard exclusive reac- 
tions jSSlESI, such as deeply virtual Compton scattering (DVCS) pUl 1^ 13^ IH^ ^^ 
and deeply exclusive production of mesons (DMP) [331 ^^ IS7j . In particular, the 
DVCS process attracted a lot of attention, both theoretical and experimental. 

The thesis is organized as follows. In Chapter II, we give a short review on the DIS 
process and discuss its relation to the forward virtual Compton scattering amplitude 
(VCA). In addition, we introduce parton distribution functions within the context 
of the parton model. The commonly used phenomenological functions, such as form 
factors, PDFs, distribution amplitudes (DAs) and GPDs are the subject of study 
in Chapter III. We list the definitions and discuss some of theoretical aspects, i.e. 
the basic properties of these functions and their relations to one another. Typically, 
GPDs can be measured in hard exclusive leptoproduction processes. DVCS, as the 
simplest process in this respect, is discussed in Chapter IV. The VCA is calculated 
at the leading twist within the framework of the nonlocal light-cone expansion of 
the product of currents in QCD string operators. Moreover, a simple model for 
the nucleon valence GPDs (i.e. those, which do not include the contribution from 
sea quarks) is introduced. In Chapter V, we apply the parton model to study the 
inclusive photoproduction of lepton pairs. Chapter VI is devoted to the exclusive 
version of the same process. It is related to time-like Compton scattering (TCS), that 
is the inverse of the DVCS process. With the help of a different set of GPDs, which 
cannot be accessible in the standard electromagnetic DVCS process, the extension 
into the weak interaction sector is made possible. This work is presented in Chapter 
VII. In Chapter VIII, we turn our attention to wide-angle real Compton scattering 
(WAGS). Our approach is based on the handbag dominance, in other words, the 
light-cone singularities of the Compton scattering amplitude again play an important 
role EH HO] . Here we use the formalism of double distributions (DDs) to describe 
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nonforward matrix elements of light-cone operators. We study both the leading and 
next-to-leading twist contributions. Finally, we draw conclusions and outline future 
research plans in Chapter IX. 
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CHAPTER II 
DEEPLY INELASTIC LEPTON SCATTERING 

11.1 INTRODUCTION 

We start with phenomenologically most important of the hard processes, inclusive 
scattering of high-energy (charged and neutral) leptons on hadrons (to be specific, 
we will consider here and throughout this thesis, except in Chapter VIII, a spin-1/2 
hadronic target, e.g. a nucleon) with the exchange of vector bosons (photons, 
or Z^). It played the key role in revealing the quark structure of hadrons, and it 
is, in addition to electron-positron annihilation into hadrons, the simplest process 
involving strongly interacting particles. Since leptons do not possess a resolvable 
internal structure, the reaction cross section depends solely on the internal structure 
of hadrons. 

In Section II. 2, we discuss the kinematics and the relation between the hadron 
tensor and the forward VGA. Both the electron and neutrino-induced DIS cross sec- 
tions are given in terms of the relevant structure functions. The parton model is 
introduced in Section II. 3. With the help of the parton model master formula, we 
calculate the DIS cross section, and express the structure functions in terms of PDFs. 

11.2 KINEMATICS 

Consider the electron scattering process off the nucleon target through a single photon 
exchange, 

e-ik) + N{P) e-(A;')+X, (1) 

as illustrated in Fig. ^ In other words, we consider only the quantum electrodynamic 
(or shortly QED) interaction between the electron and the nucleon, and further keep 
only its lowest order. We denote the nucleon four-momentum with P and the initial 
and final momenta of electrons with k and k', respectively. In the laboratory frame, 
in which the target is at rest, the nucleon four-momentum is P = (m, 0^, where 
M stands for the nucleon mass. The electron four-momentum in the initial state is 
k = (uj,k), and in the final state is k' = (u',k'). The momentum transfer to the 
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FIG. 1: Deeply inelastic electron-nucleon scattering in the one-photon exchange ap- 
proximation. 

target, carried by the virtual photon, is q = k — k' = {u, q) . In this scattering process, 
due to the large q, the nucleon breaks up and forms an infinite number of possible 
hadronic final states labeled by X, which remain unobserved. Thus only the outgoing 
electron is detected. One measures its energy u' and the scattering angle 6 relative 
to the incident beam of fixed energy u. 

There are two characteristic Lorentz-invariant quantities of the process, namely, 
the virtuality of the photon (since it is space-like, < 0, one conventionally defines 
a positive quantity), 

q' ^ -Q\ (2) 

which fixes the spatial resolution of the scattering process, A ~ and the 

invariant mass of the final hadronic state, 

= {P + qf = M^ + 2{P-q)-Q\ (3) 

In the laboratory frame, they are given by 

Q2 = Aujuj'sm^{e/2) (4) 

and 

W"^ = M"^ + 2Miy-Q^. (5) 

Note that in Eq. (0)) we have neglected the electron mass compared to its energy. 
Unlike elastic scattering, for which W"^ = and the momentum transfer squared 
is connected by the relation u = Q'^/2M to the energy loss of the electron u, we 
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now have two independent variables, v and . The inclusive differential cross section 
(see, e.g. El 1121 ESj) factors into the electron tensor and the hadron tensor 

where dQ is the solid angle into which the electron is scattered and a = e'^/An ~ 
1/137 is the electromagnetic fine structure constant (in the following we use the 
Heaviside-Lorentz convention together with the natural choice of units, frequently 
used in particle physics, in which h = c = 1). The electron tensor can be computed 
exactly in QED. For unpolarized electron scattering, the expression for the spin- 
averaged electron tensor is 



s,s' 



= 2[k^k', + Kk'^~g^,{k-k')\. (7) 

The hadron tensor, on the other hand, describes the response of the nucleon and 
therefore, it should include all possible transitions of the nucleon from its ground state 
\N{P,S)) to any hadronic final state \X{P')). It is worth noting that at present, 
the nonperturbative nature of strong interactions prevents us from calculating the 
hadron tensor directly within the framework of QCD. Formally, it can be written as 

x(X(P')U^Af(0)|iV(P,5)), (8) 



where Jem is the quark electromagnetic current. Using the translation invariance, 

{N{P,S)\JUA^)\X{P')) = {N {P,S)\JU,mX{P'))e-^^^'-^y\ (9) 
and the completeness condition of the states \X {P')), 

j:\xin) {x{p')\ = 1, (10) 

X 

the hadron tensor can be presented in a more compact form, namely, 



^Jd'ze^^^-^lY.{N{P,S)\J^EMi^)JEMiO)\N{P,S)). (11) 
Adopting a simplifying notation. 
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{N{P)\0\NiP)) ^ W{NiP^S)\0\N{P,S)), (12) 

^ s 
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where O is an arbitrary operator we have 

= ^ld'ze^^'^-HN{P)\r^,Az)rEMiO)\N{P)). (13) 

Interchanging the currents inside the matrix element leads to a delta function 
5(4) (^p' _ p which, in the laboratory frame, requires that E' = M — u < M. 
Note that such a state \X [P')) does not exist since the nucleon \N{P^S)) is the 
ground state baryon. Hence the delta function cannot be satisfied and the expression 
with the interchanged currents vanishes. By adding this vanishing matrix element, 
one can also rewrite the hadron tensor as a current commutator, 

= ^Jd*ze^(^-^{N{P)\[r^,Az),JEMm\N{P)). (14) 

Next we discuss the relation between W'^'^ and the object known as the forward 
virtual Compton scattering amplitude (VGA). The latter, shown in Fig. |21 is defined 
as a Fourier transform of the correlation function of two electromagnetic currents, 

T"-" = ijd'ze^^'^-HN{P,S)\T{J^E^{z)r^^,m\N{P,S)). (15) 

The time-ordering symbol T inside the matrix element instructs us to place the 
operators into chronological order with the operator having the later time argument 
to the left. The amplitude (fT^ is averaged over the nucleon spin. It corresponds 
to the scattering amplitude for the Compton process off the proton in the forward 
direction, when contracted with the polarization vectors of the off-shell photon of 
momentum g, 7^ 0, 

or 

T(7P^7P) = ehl{q)tAq)T^''. (17) 

Here e denotes the electric charge of the electron. With the help of the optical 
theorem, the hadron tensor can be written as 

W^""" = -^QT^"'. (18) 
27rM ^ ^ 

Lastly, it is customary to express the hadron tensor, see Eqs. (fT^ and (IHj), in 
terms of the so-called structure functions. The tensor expression is Lorentz- invariant. 
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FIG. 2: Forward virtual Compton scattering amplitude. The dashed line denotes the 
cut, which gives the imaginary part, and it is further related to the hadron tensor 
through the optical theorem. 



and it should also be invariant under parity and time reversal. Together with the 
electromagnetic current conservation, d^J^^^j = 0, implying the electromagnetic gauge 
invariance. 



q^Wf"" = and W^^'q^ = 0, 



(19) 



one arrives at the most general form for W'^'^ in the case of unpolarized scattering. 
It is determined only by two independent scalar structure functions (or alternatively 
the response functions). 



M2 



9"" + 



pu 



(20) 



The hadron tensor (|2(Jj) is symmetric with respect to Lorentz indices /i and v. In the 
case of polarized scattering the tensor has an extra term, which is antisymmetric with 
respect to /i and v and contains two additional structure functions. We notice that 
both functions Wi and W2 depend on two invariant variables, v = {P ■ q) /M and Q^. 
Substituting Eqs. ((Tj) and ((201) i'^to Eq. (jHI) and using Eq. (jH), gives the differential 
cross section in the laboratory frame in the form 



(fa, 



(eiV) 



a 



dQdu' 



4a;2sin^ (9/2) 



2Wi (u, sin^ (9/2) + W2 U Q^) cos^ (9/2) 



(21) 



Usually, u is replaced by another Lorentz-invariant quantity, the Bjorken scaling 
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variable, 

and further the dimensionfuU VF-structure functions by their dimensionless partners 
defined as 

Fi(xB,g') = MWi{u,Q^), 

F2{xb,Q') = uW2{u,Q'). (23) 

Since the invariant mass of the unobserved hadronic final state is always larger than 
the nucleon mass, it follows from Eqs. (jH)) and ()22|) that 

In particular, the kinematical point = 1 corresponds to elastic electron-nucleon 
scattering. In addition to x^, there is another dimensionless variable frequently used, 
the so-called inelasticity parameter, 

^ ^ my 

It specifies, in the laboratory frame, the fraction of the electron energy that is trans- 
ferred to the nucleon, y = 1 — uj' /u, and hence it is also kinematically constrained to 
the region (0, 1]. The invariants Q"^, xb and y are related to each other through 

Q' = XBy{s-M^), (26) 

where s = {k + P)^ is the usual Mandelstam invariant of the scattering process. The 
last formula is obtained by writing the Bjorken scaling variable in the laboratory 
frame, xb = Q^/2Mz/. Sometimes it is convenient to present the cross section in the 
invariant form. After some algebra, the cross section results, in terms of dimensionless 
structure functions, into 

(P(JieN) _ / 2nMujy \ (Pa^eN) 
dxBdy \ uj' / dVtduj' 



2 



47ra 



(27) 



For the sake of completeness, let us investigate the weak version of inclusive lepton- 
nucleon scattering. In general, the weak boson is used as a probe to study the nucleon 



11 



structure. To be specific, we shall consider only the charged current case using muon 
neutrinos. The latter convert into muons by the emission of a boson. The 
reaction is 

iy^{k)+N{P) fi-{k')+X. (28) 

Recall that the incident neutrino is always left-handed and for that reason, there is 
no averaging over the initial spin. Moreover, having the initial lepton (i.e. neutrino) 
state polarized yields an extra term in the lepton tensor that is antisymmetric in 
Lorentz indices and u, 

■^[iJ ~ S ^^J^u ~^ ^i^^fi Q^u {k ■ k ) -\- i^^jyfjT-k k . (29) 

In constructing the corresponding hadron tensor we have to respect that the parity in 
the weak interaction is no longer conserved. Thus, due to the presence of the parity- 
violating term ie^'^°'^ PaQp in the hadron tensor, a third response function {u, Q"^) 
arises in the expression for the cross section. In analogy with Eq. ©, along with the 
factor replacement a/Q"^ Gpl^'n^, where Gp = g^/4:V2M^ ~ 1.166-10-^ GeV"^ 
is the so-called Fermi constant with g being the weak interaction constant and 
the mass of the VF-boson (note that one chooses the coupling constant of the W- 
boson to the nucleon as (7/2 v^, and that the factor 1/M^ in the definition of Gp 
comes from the VF-boson propagator considered in the limit <^ in other 

words, one deals with the Fermi contact interaction), the laboratory differential cross 
section is 



dVtdu' 27r2 

UJ + Uj' 

+ 



2M 



2iy (z/, Q2) sin2 (^/2) + Wlf""^ (z/, Q^) cos' {6/2) 
Wt""^ {u,Q') sin' (9/2) 



(30) 



Finally, in terms of the Lorentz invariant variables Xb and y, the invariant cross 
section reads 



dxBdy 2tx 



XBv'Ft''^ {xB, Q') + {l-y) Ft""^ {xB, Q' 



(31) 



Similarly to the electromagnetic case, the dimensionless structure functions in the 
weak interaction sector are defined as follows 



Fr^{xB,Q') ^ Mwr\^,Q') 
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Fr^(xB,Q^) - M'^'^l^.Q')- (32) 

II.3 BJORKEN SCALING AND THE PARTON MODEL 

The nucleoli structure with the size of the order A < 0.2 fm can only be resolved at 
sufficiently large >1 GeV^. Having large also implies large values of i/, more 
precisely v"^ ^ Q^, since the ratio xb = jlMv is finite and bounded between 
and 1. Accordingly, a massive hadronic state X is produced with the invariant mass 
equal to 

W'^ = Q^1~xb)/xb + M^:^M^. (33) 

It lies well above the resonance region, W > 2 GeV. In this kinematical regime, 
known as the deeply inelastic region, the structure functions F{xbiQ'^) extracted 
from the measured inelastic cross sections do not depend significantly on Q"^ but 
rather only on xb- This interesting feature is termed Bjorken scaling [HI I45| 1^ I47j . 
It can be stated in the following way: in the Bjorken limit, where both ^ oo and 
— > oo whereas is kept fixed, one experimentally observes 

f(xb,Q2) F{xb). (34) 

The fact that, for sufficiently large values of Q^, the structure functions are inde- 
pendent of implies that the nucleon, or any other hadron, is made of point-like 
constituents. Recall that a finite size object must have a form factor and hence intro- 
duce some dependence on Q^. Furthermore, since the structure functions are Lorentz 
invariant, one can study the scattering process in any reference frame. However, the 
description is considerably simplified, if we look at the nucleon in a very fast moving 
system with its momentum approaching to infinity along the z-direction, i.e. in the 
so-called infinite momentum frame (the physics of the process is, of course, indepen- 
dent of this choice). Then the transverse momenta, the rest masses of the nucleon 
constituents and, for consistency, the nucleon mass can be neglected and accordingly, 
the nucleon momentum is = (-Pj, 0^, -Pz)- In other words, the structure of the 
nucleon is described only in terms of the longitudinal momenta of its constituents. 
This is the basis of the parton model of Feynman j3] , which gives the clearest physical 
interpretation of scaling. In the parton picture, the nucleon is viewed as a collection 
of noninteracting, point-like constituents, the partons. The interaction of the electron 
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with the nucleon can be then viewed as the incoherent sum of interactions between 
the electron and the individual partons. This approximation is valid as long as the 
duration of the electron-parton interaction, which is regarded as elastic scattering, is 
so short that the interaction between the partons themselves can be safely neglected. 
Thus we consider the scattering process in the impulse approximation by picking up 
only the lowest-order electromagnetic contribution, and neglect all the QCD correc- 
tions associated with the exchange or emission of gluons. 

In the infinite momentum frame, a given parton is characterized by the longitu- 
dinal momentum fraction x e [0, 1] of the total nucleon momentum P, 

p = xP. (35) 

In addition, for each parton species a, we define the parton distribution function 
(PDF) fa/N {x). It describes the probability of finding a parton of type a at the lon- 
gitudinal momentum fraction x inside the target nucleon N. Note that these functions 
cannot be computed using QCD perturbation theory. All fractions x have to add up 
to 1 and hence the normalization 

[ dx xfa/N (x) = 1 (36) 

holds. The scattered parton has the momentum p' — p + q, given by the four- 
momentum conservation, where q is the momentum of the virtual photon. Since the 
parton is on its mass shell, it follows that 

p'2 = {xP + qf 

= x'^M^ + 2x{P-q)-Q^ 

fa 0. (37) 

Neglecting x'^M^ compared to and u (recall that {P ■ q) = Mi>) yields 

^ = 2(^^^^- 

The longitudinal momentum fraction is found to be identical to the Bjorken scaling 
variable. It means that the parton must have the fraction xb of the nucleon momen- 
tum in order to absorb the virtual photon. It should be emphasized, however, that 
the variable xb has this simple meaning only in the infinite momentum frame. 

Due to the fact that the partons are point-like and noninteracting, the DIS cross 
section for nucleon-electron scattering in the parton picture is simply given by the 
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e-(k) 



N(P) 




FIG. 3: Kinematics of deeply inelastic electron-nucleon scattering in the parton 
model. 



incoherent sum over all the contributing partons (i.e. the partonic cross sections 
for elastic scattering of an electron from the individual partons), weighted with the 
proper distribution functions. The fundamental relation, known as the parton model 
master formula, reads 

a \e- [k) N (P) ^ e {k') x] = E /' dx fa/N {x) 

xa [e- (k) qa (xP) ^ {k') {p')] , (39) 

where the sum runs over quarks and antiquarks, in order words, over all charged 
partons because they are the ones that interact with the virtual photon. Schemati- 
cally, see Fig. El the hadronic interaction is broken into a PDF part, represented by 
the blob, and hard (parton-electron) scattering part. According to the factorization, 
PDFs do not interfere with the hard scattering part. For that reason, they are uni- 
versal in a sense that they are same for all inclusive hard scattering processes, not 
only for electromagnetic DIS j3H]- 

Let us work out the explicit leading-order formula for the DIS cross section using 
Eq. ()39|) . First we calculate, at the most elementary level, the cross section for 
electron scattering from a single parton. For <^ — where u = {k' — P)^ 
is the Mandelstam variable, one has 

Q"^ = xbVS and y = , (40) 

s 

and the invariant form of the partonic double differential cross section is (see, e.g. 

miiini) 



l + {l-yf xbS{x-xb). (41) 
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The electric charge of the parton of type a in units of |e| is denoted by Qa- Substituting 
Eq. (gH) into Eq. (jSHI) results in 



dxsdy yl + {l-y)\i:Ql^Bfa,N{xs). 

Comparing now Eqs. (jl^ and (j2ZI); under the assumption that <^ s, gives 



xbV^Fi (xb, Q2) + (1 - y) F2 (xb, Q' 



We notice that the right-hand side of Eq. 



(42) 



1 + (1 - 2/) J H QlxBfa/N (xb) ■ 
a 

(43) 

depends only on xs and y. Therefore, 



or, 



XBV^Fi (xB, Q') + il-y) F2 [xB, g') = const. for all (44) 
Fi[xB,Q^) = Fi{xb), 

F2{XB,Q^) = F2{xb). (45) 

Hence the structure functions in the parton model exhibit the scaling property, i.e. 
they do not depend on Q^. Conversely, Bjorken scaling suggests the asymptotically 
free quark dynamics, the property that the quark interaction gets weaker at short 
distances. In addition, by comparing powers of y in Eq. (|43|) . we can express the 
F-structure functions in terms of PDFs, namely, 

1 



Fi (xb) 



F2ixB) = J2Ql^Bfa/N (xb) 



(46) 



Finally, from Eq. PU)) one can read off the following simple relation specific to the 
scattering of electrons from massless fermions 



F2 (xb) 



2xbFi (xb) ■ 



(47) 



It is known as the Callan-Gross relation jHO!, and it is an important evidence that 
the partons detected in DIS are indeed the spin-1/2 quarks of hadron spectroscopy 
[15] . It is worth noting that a more descriptive notation is commonly used, in which 
fa/N (x) is replaced by the flavor labels Un (x) , ^at (x) , Sat (x) , Mat (x), etc. Then the 
proton structure function Fi is 



1 r4 _ 1 r 

Fipix) = -<^-[upix)+Up{x)] + -[dp{x) + dp{x) 
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(x) + Sp (x)] 



(4J 
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and similarly for (x) . 

In summary, the parton model serves both to explain scaling and to identify 
partons as quarks. However, it is important to add that ultimately, the success of 
this simple model must be justified in quantum field theory, in particular within 
the framework of the QCD operator product expansion (OPE). It turns out that, in 
the Bjorken regime, the dominant contribution to W^''', see Eq. (jl3|) . comes from 
the integration region < < const. /Q^. Accordingly, the DIS cross section is 
dominated by the light-cone region, — 0, of the space-time integration in Eq. (fTSj) 
|51j . In fact, the scattering amplitudes of other hard inclusive and exclusive processes 
are also governed by the product of currents (z) J'^ (0) near the light-cone. Thus 
OPE is a powerful tool which, not only recovers the parton model predictions but 
also allows to analyze, in a systematic way, the terms in the light-cone expansion 
contributing to a given power of 1/Q^. 
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CHAPTER III 
PHENOMENOLOGICAL FUNCTIONS 

111.1 INTRODUCTION 

The fundamental particles, from which hadrons are built, are known. They are quarks 
and gluons. The interactions between them are described by the QCD Lagrangian, 
which is also established. Unfortunately, knowing the first principles is not sufficient. 
We still need to understand how QCD works, in other words, how to translate the 
information obtained from experiments on the hadronic level into the language of 
quark and gluon fields. One may, for example, consider projecting these fields onto 
hadronic states. The resulting matrix elements can be interpreted as hadronic wave 
functions |^ I53j . An alternative approach to describe the hadronic structure is to 
use different phenomenological functions. The well-known examples are form factors, 
usual parton distributions functions and distribution amplitudes. Since they have 
been around for a long time, they are termed the old phenomenological functions. 
We discuss them separately in Sections III. 2, III. 3 and III. 4. On the other hand, 
the concept of generalized parton distributions [23 ISOl HH IS21 ISSl IMI (for reviews, 
see |54[ I36[ 133] and recently j|56j) is new. These new phenomenological functions are 
hybrids of the old ones and therefore, provide a unified and more detailed description 
of the hadronic structure. Different species of generalized parton distributions are 
presented in Section III. 5 together with some of their general properties. 

111.2 FORM FACTORS 

Form factors are defined through the matrix elements of electromagnetic and weak 
(neutral and charged) currents between the hadronic states. In particular, the matrix 
element of the electromagnetic current (for the weak currents, see Section VII. 3), 

^^m(0) = E^?/^/ (0)7^^/(0), (49) 
/ 

where Qf is the electric charge (in units of |e|) of the quark of fiavor /, between the 
nucleon states N {pi, si) and N (p2, S2) is parametrized in terms of two independent 
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nucleon electromagnetic form factors known as the Dirac and Pauli form factors. 
Namely, 



{N{p2,S2)\J^MiO)\N{pi,Si)) = U{P2,S2) 

xu{pi,si), (50) 



rF,it)-ta'^-^^F2it) 



where u {p2, S2) and u {pi, Si) are the Dirac spinors, r = pi — p2 is the overall mo- 
mentum transfer and the invariant t = r'^. The elastic form factors should not be 
confused with the dimensionless structure functions Fi and F2 {xb,Q^) in- 

troduced in the preceding chapter. Note also that for elastic scattering, we have 
= —t. The nucleon mass M in Eq. ()5U|) is introduced only for dimensional con- 
venience. Similarly to the flavor decomposition of the electromagnetic current given 
by Eq. (j49p . the Dirac and Pauli form factors can also be expressed in terms of their 
flavor components. 

Flit) = EQfFifit), 

f 

F2{t) = Y.QfF2f{t). (51) 
/ 

Their limiting values at t = are known. The Dirac form factor gives the total 
electric charge of the nucleon N, Fi (t = 0) = Qn, and the Pauli form factor gives its 
anomalous magnetic moment, Fi [t = 0) = kn- For the proton, we have Hp = 1.793, 
and for the neutron = —1.913. The Q^- dependence of form factors is a clear 
evidence for the extended structure of the nucleon, in particular the charge and current 
distributions. 

Note that writing the matrix element (A^ {p2, S2)\ Jem (0) 1^ (Pi) -^i)) in the most 
general form, one might expect to include, in addition to terms 7^^ and a^^r^, also 
terms like {pi + P2Y ■, and a^'^ {pi +P2)y- However, with the help of the Gordon 
identity, 

u{p2,S2)-i^u{pi,si) = ^u{p2,S2)[ipi+P2T -i(r^''r^]u{pi,si) , (52) 

one can express the terms {pi + ^2)'^ as the linear combination of terms 7^ and a^^Ty. 
Moreover, the electromagnetic current conservation, r^J'^j^^ = 0, implies that the 
term should vanish and hence the most general form reduces to the expression 

dsni). 
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Instead of Fi and F2, we often introduce the so-called Sachs electric and magnetic 
form factors, Ge and Gm, respectively. Then the matrix element takes the form 

^GE{t)-{t/4M')GM{t) 



{N{p2,S2)\J^MW\NiPus^)) = uip2,S2 



r 



-la 



1 -t/4M2 
r. Gm (t) - Ge (t) 



2M 1 - t/4M2 



u {pi, s 



1) , 



where 



and 



(53) 



(54) 



Guit) = Fi(t) + F2(t). 



These form factors are normalized at = in the following way 



GEpit = 0) = 1, 

GMpit = 0) = /ip = 2.793 



for the proton and 



Geu {t = 0) 
Gmu (^ = 0) 



0, 



-1.913 



(55) 



(56) 



(57) 



for the neutron. Their magnetic moments are given in terms of the nuclear magneton, 
fx = \e\ /2Mp = 5.051 ■ lO'^^ Am^ in the SI units. 

The nucleon electromagnetic form factors can be measured through elastic 
electron-nucleon scattering. 



e' ik) + N{p^) 



e- {k') + N{p2). 



(5^ 



The process is shown in the one-photon exchange approximation in Fig. EJ 



III.3 PARTON DISTRIBUTION FUNCTIONS 

Parton distribution functions (PDFs) are defined through the forward matrix elements 
of the light-like correlation functions, i.e. the quark and gluon fields separated by 
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e-(k) 



Nip,) 




N(pJ 



FIG. 4: Elastic electron-nucleon scattering in the one-photon exchange approxima- 
tion. 



hght-hke distances. For the unpolarized case we have 



{N (P, S)\^Pf {-z/2) rijf {z/2) \N (P, S)) 



z^=0 



u (P, S) ^^u (P, S) 



X / dx 
lo 



(59) 



and for the polarized one 



{N (P, 5)1 ijf {-z/2) Yl^i^f {z/2) \N (P, S))^,^, = u (P, S) Yl^u (P, S) 



X dx 
Jo 



e"(^-")A/^ (x) + e-*^^^-"^A/^ (x) . (60) 



ix{P-z) 



To make connection with GPDs, which are usually discussed in the region — 1 < x < 1 
(the variable x runs from to 1 for quarks and from —1 to for antiquarks), it is 
convenient to introduce new distribution functions, 



In {x) 



Jn {x) X > 
-In {-x) X < 



and 



A/jv {x) 



AfN {x) X > 
A/at (-x) X < 

and alternatively, write the integrals over x in Eqs. fl59p and ()60|) as 



(61) 



(62) 



dx 



^i<P-z)f^ (x) - e-"^^-"7iv {x) 



_ -ix{P-z) . 



dx e^^(^-")/ 



N [X) 



dx 



^ix(P-z) 



AfM {x) + e-*^(^-^) A/^ (x) 



dx e'"(^-")A/ 



N [X 



(63) 
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Furthermore, we observe that the definition of PDFs has the form of the plane- 
wave decomposition. Thus it allows to give the momentum-space interpretation. For 
example, (x) (^Jn (x)^ is the probability to find the quark (antiquark) of fiavor / 
carrying the momentum xP inside a fast-moving nucleon N having the momentum 



Parton distribution functions have been intensively studied in hard inclusive pro- 
cesses for the last three decades. The classic example in this respect is the deeply 
inelastic scattering process shown in Fig. ^ Its structure functions are directly ex- 
pressed in terms of PDFs. By substituting Eq. (fTHj) into the expression (0), it is easy 
to see that the DIS cross section is given, via the optical theorem, by the imaginary 
part of the forward virtual Compton scattering amplitude. 

Let us consider this amplitude in the specific kinematics known as the Bjorken 
limit. Here the invariant momentum transfer to the nucleon system is sufficiently 
large, —q^ = Q"^ ^ oo, together with large total center-of-mass energy of the virtual 
photon-nucleon system, s = {P + q)^ —>■ oo, while the Bjorken ratio Xb = <5^/2 {P ■ q) 
is finite. In the deeply inelastic region discussed in Section II. 3, we have, for instance, 

> 1 GeV^ and W"^ = s > A GeV^. Hence in this particular regime, the behavior 
of the forward VGA is dominated by short distances, i.e. when the separation be- 
tween the two point-like photon-quark vertices (note that the photons couple to the 
quarks of the nucleon) in the amplitude is light-like. As a result, QGD factorization 
works and the amplitude factorizes into a convolution of a perturbatively calculable 
hard scattering process at the level of quarks and gluons, and process independent 
matrix elements, containing the soft nonperturbative information about the nucleon 
structure. These matrix elements are parametrized in terms of PDFs. Schematically, 
factorization allows to write the leading-order amplitude in the form of the so-called 
handbag diagrams. In these diagrams, two photons couple to the same quark line, 
as illustrated in Fig. |3| It is worth noting at this point that throughout this thesis 
we work only to the leading order in the strong coupling as- Perturbative corrections 
produce logarithmic dependence of PDFs on the scale Q^, in other words, they define 
the evolution of parton distributions, which can be calculated in QGD [HZl EHJ EH] ■ 

Now taking the imaginary part of the forward VGA generates the delta function. 



X. Thus in the DIS process we measure parton distribution functions fjs (x) and 



P. 




(64) 



which selects two points, x = ±xb, after the integration over the momentum fraction 
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FIG. 5: Deeply inelastic electron-nucleon scattering cross section as the imaginary 
part of the forward virtual Compton scattering amplitude. The summation over X 
reflects the inclusive nature of the nucleon structure description by parton distribution 
functions. 



^Jn {x) at X = Xb corresponding to the quark PDFs, and at x = —xb for those of 
antiquarks. Unlike the form factors, one deals, in the case of the parton distribution 
functions, with a light-like separation instead of a point vertex, and also the initial 
and final nucleon momenta are equal. 

In summary, the unpolarized DIS experiments have mapped out the quark and 
gluon distributions in the nucleon while the polarized DIS experiments have shown 
that quarks carry a small fraction of the nucleon spin. As a result, new investigations 
to understand the nucleon spin became necessary. 



III.4 



DISTRIBUTION AMPLITUDES 



Distribution amplitudes (DAs) (sometimes also referred to as the hadronic wave func- 
tions) [nni IMl in21 inni ini] describe hadrons in hard exclusive scattering processes and 
therefore, in addition to usual parton distribution functions, provide complementary 
information about the hadronic structure. They are defined through the vacuum-to- 
hadron matrix elements (0| ... |P) of light-cone operators. For example, in the pion 
case we write 



iP^U f da e^"(^-")/2 



V9^+ {a) 



(65) 
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FIG. 6: The handbag contribution to the process 7*7 tv^. The hght-cone domi- 
nance is secured by the high virtuahty of the incoming photon 7*. 



The fractions of the pion momentum carried by the quarks are (1 ± a) /2. One can 
interpret ^9^+ (a) as the probabihty amphtude to find a positive fast-moving pion vr"*" 
in a quark-antiquark state ud, with the longitudinal pion momentum P shared in 
fractions (1 + a) /2 and (1 — a) /2. The distribution ip^+ (a) is an even function in 
the relative fraction a |54j . 

The simplest and cleanest process, in which the pion DA can be accessed, is the 
transition |M1 lUH] 

l*{q)+l{q') vr°(P). (66) 

For large virtuality Q^, the leading-order contribution to the amplitude is given by 
the handbag diagrams depicted in Fig. IHl 



III.5 GENERALIZED PARTON DISTRIBUTIONS 

The key idea of the GPD approach is the so-called hybridization. As hybrids of form 
factors, PDFs and DAs, generalized parton distributions provide the most complete 
information about the hadronic structure. Parton distributions parametrize the lon- 
gitudinal momentum distributions (in an infinite momentum frame) of partons in 
the hadron while the Fourier transforms of form factors in impact parameter space 
(impact parameter measures the transverse distance of the struck parton from the 
hadron center) describe the transverse coordinate distributions of the hadron's con- 
stituents jnniEZI- Generalized parton distributions, on the other hand, encapsulate at 
the same time the longitudinal momentum and transverse coordinate distributions. 
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thereby providing a much more detailed and comprehensive, three-dimensional snap- 
shot of the substructure of the hadron. In addition, the universality of GPDs allows to 
develop a unified description of wide variety of different hard processes, both inclusive 
and exclusive. 

Generalized parton distributions can be, in general, divided into two groups, 
namely, the skewed parton distributions (SPDs) and the double distributions (DDs). 

III. 5.1 Skewed Parton Distributions 

There are two implementations of the SPD formalism, as illustrated in Fig. [7| The 
nonforward parton distributions (NFPDs) ■34^ have the advantage of using the vari- 
ables similar to those of the usual PDFs. The distributions depend on X, the fraction 
of the plus component of the light-cone momentum P+ of the hadron carried by the 
parton (the plus and minus light-cone components are defined by = (a° ± a^) /a/2 
for any Lorentz four- vector a), or alternatively, the longitudinal momentum fraction 
with respect to the initial hadron momentum P; on the skewness parameter ( spec- 
ifying the difference between the initial and final hadron plus momenta, r"*" = CP~^] 
and on the invariant momentum transfer t = r'^. For instance, the nonforward par- 
ton distribution (X, t) is the probability amplitude that the initial fast-moving 
hadron, having longitudinal momentum P"*", emits a parton of flavor / carrying the 
momentum XP~^ while the final hadron, having longitudinal momentum (1 — (^) 
absorbs a parton of flavor / carrying the momentum (X — () P'^. In this particular 
scheme, the initial and final hadron momenta are not treated symmetrically. The 
off-forward parton distributions (OFPDs) [SDIE]! on the other hand, use symmetric 
variables expressed in terms of the average hadron momentum, p = {pi + P2) /2, with 
Pi being the momentum of the initial hadron and p2 the momentum of the final one. 
In the symmetric scheme the hadron longitudinal (the plus component) momenta are 
(1 ± and accordingly, those of the active partons become {x ± ^)p'^. Similarly 
to NFPDs, the off-forward distributions are defined for each quark flavor, and are 
the functions of three variables, namely, the light-cone momentum fraction x, the 
skewness, ^ = r+/2p+, (here the skewness is introduced as the coefficient of propor- 
tionality between the light-cone plus components of the momentum transfer and the 
average hadron momentum, and like ( varies between and 1) and the invariant t. 
In fact, one should bear in mind that both the nonforward and off-forward parton 
distributions also depend weakly (i.e. logarithmically) on the probing scale Q^. 
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p (1-VP (i+^)p (i-^)p 

FIG. 7: Comparison between the descriptions in terms of nonforward and off- forward 
parton distribution functions. 



The two pairs of variables {X, () and (x, ^) are related to each other. The conver- 
sions between them are 

X ^ ^ (67) 

and 

We emphasize that these variables solely characterize the plus, or the longitudinal, 
components of the momenta of the partons involved, however, their transverse mo- 
menta are related to the t-dependence of GPDs. The latter is driven by the t- 
dependence of the corresponding elastic form factors. Thus it is possible to access 
simultaneously the longitudinal momentum and transverse position of the parton in 
the infinite momentum frame [Uni EZ] • Moreover, by removing the parton with the 
light-cone momentum fraction x + ^ from the hadron and replacing it, at some later 
point on the light-cone, with the parton of the momentum fraction x — one can say 
that SPDs (or GPDs, in general) measure the coherence between two different parton 
momentum states of the hadron whereas usual PDFs yield only the probability that 
a parton carries a fraction x of the hadron momentum. In addition to the parton 
momentum as well as the spin correlations, we can even consider the matrix elements 
corresponding to different hadrons in the initial and final states, e.g. the proton-to- 
neutron transition accessible through the exclusive charged pion electroproduction, 
the proton-to-A transition in the kaon electroproduction, the nucleon-to-delta transi- 
tion, and hence study flavor nondiagonal GPDs. 

Let us focus on the off-forward parton distributions. At the leading, twist-2 
level, the hadron structure information can be parametrized in terms of two un- 
polarized OFPDs, Hf {x,C,,t) and Ef (x, and two polarized OFPDs, Hf {x,^,t) 
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and Ef (x, Since — 1 < x < 1 in Fig. [7J the momentum fractions x ± ^ of the 
active partons can be either positive or negative. Positive and negative momentum 
fractions correspond to quarks and antiquarks, respectively. Therefore, each OFPD 
has three distinct regions. When ^ < x < 1, both partons represent quarks while for 
— 1 < X < —C,, they are both antiquarks. In these two regions, the distributions are 
just a generalization of the usual PDFs. In the central region, < x < ^, which is of- 
ten referred to as the mesonic region, the parton with a positive momentum (x + ^) p"*" 
is going out from the blob and represents a quark. The returning parton has a negative 
momentum and therefore, should be treated as an outgoing antiquark with the mo- 
mentum — x)p+. The total momentum of the quark- antiquark pair, r"*" = 2^p"'", is 
then shared in fractions (x + ^) = (1 + x/^) /2 and {C, — x) = (1 — x/^) r+/2. 
In this region of x, which is not present in deeply inelastic scattering, OFPDs behave 
like meson distribution amplitudes with a = x/C,- 

Clearly, in the nonforward kinematics, SPDs uncover much richer information 
about hadronic structure, which is not accessible in the DIS process. This new infor- 
mation can be extracted with the study of hard exclusive processes, such as deeply 
exclusive photon or meson electroproduction (one refers to the former as deeply virtual 
Compton scattering or shortly DVCS), which turns out to be a much more difficult 
task due to the small cross sections. Nevertheless, high-energy and high-luminosity 
electron accelerators combined with large acceptance spectrometers give a unique 
opportunity to perform precision studies of such reactions. 

The factorization into short and long distance dynamics is more general. Having 
large space-like virtuality of the initial photon is sufficient for QCD factorization to 
work [211 EH EHl inni- In particular, in the DVCS process, which will be studied in 
detail in Chapter IV, the initial photon is highly virtual, —qf — > oo, while the final 
photon is on shell, g| = 0. In the leading-twist handbag approximation (DVCS is a 
handbag-dominated process for —qf as low as 2 GeV^) illustrated in Fig. |S1 the hard 
short- distance part of the so-called nonforward virtual Compton scattering amplitude 
factorizes from the nonperturbative long-distance part. The latter is represented by 
the lower blob and contains now the nonforward matrix elements (A^ {P2) \ ••• |^ (Pi)) 
of the same quark and gluon operators as in the forward case. These matrix elements, 
describing the nucleon structure, are parametrized in terms of GPDs. Even though 
the parton picture of DVCS, see Fig. |S1 looks similar to that of DIS, see Fig. 
there are three crucial differences. In DVCS, one deals with the skewed kinematics. 
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FIG. 8: Handbag contribution to the nonforward virtual Compton scattering am- 
plitude. The hard quark propagator (this particular quark propagates along the 
light-cone) in the s- and u-channel diagrams is convoluted with generalized parton 
distributions. 



in which the plus momenta of the initial and final hadrons as well as partons are 
not equal. Furthermore, the invariant momentum transfer t in DVCS is small but 
not zero. These two extra degrees of freedom, ( (or ^) and t, make the dynamics of 
DVCS rich and diverse. Finally, in DVCS (or in any other hard exclusive process) the 
virtual Compton scattering amplitude described by GPDs appears at the amplitude 
level whereas in DIS (or in any other inclusive process) the amplitude described by 
PDFs enters through the optical theorem at the level of the cross section. 

Skewed parton distributions have interesting properties linking them to usual 
PDFs and form factors. In the forward limit, pi = p2 and r = 0, ^ = 0, t = 0, 
they reduce to PDFs obtained from the DIS process. In particular, the OFPDs Hf 
and Hf coincide with the quark density distribution (x) and the quark helicity 
distribution A/^r (x) given by Eqs. (|HT|) and (jU^ . We write the so-called reduction 
formulas, 

\ (x) X > , , 

Hfix,0,0) = { f)\ (69) 
[ -In {-x) x<0 

and 

Hf{x,0,0) = i (70) 
[ A/iv i-x) X <0, 

while the OFPDs Ef and Ef have no connections to PDFs. They are always accom- 
panied with the momentum transfer r and therefore, simply not visible in DIS. One 
can say that Ef and Ef encode completely new information on the hadron struc- 
ture, which is not accessible in inclusive measurements. Even though they have no 
analogue in the forward limit like Hf and Hf, their limits do exist. 
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In the local limit, z = 0, off-forward parton distributions reduce to the standard 
vector and axial vector form factors. In other words, the first moments of OFPDs, 
obtained by integrating the distributions over x, are equal to the nucleon elastic form 
factors (i.e. the Dirac, Pauli, axial and pseudoscalar form factors). For any C, one has 
the following relations for a particular quark flavor 

j'^dxHj{x,^,t) = F,f{t), 

I'^dx Ef{x,^,t) = F2f{t), 

J ^dx Hf {x,^,t) = gAf{t), 

j'^dxEf{x,^,t) = gpfit). (71) 

We call these relations the sum rules. It is important to note that these sum rules 
are model and ^-independent (the dependence on ^ drops out after integration over 
x). 

Generalized parton distributions are also relevant for the nucleon spin structure 
and have received considerable attention in recent years in connection with the so- 
called proton spin puzzle (for reviews, see, e.g. [701 EH IZ21)- Namely, certain low 
moments of GPDs can be related to the total angular momentum carried by quarks 
and gluons (or generically, partons) in the nucleon jSU]. In particular, in terms of 
the off-forward distributions, the second moment of the unpolarized OFPDs at t = 
gives the quark total angular momentum, 

Ji = \Y. f dxx[Hs{x,i,t = Q) + Ef{x,i,t = Q)]. (72) 

The above equation is independent of ^ (again the ^-dependence of Hj and Ef is 
removed by integration over x). The quark angular momentum, on the other hand, 
decomposes into the quark intrinsic spin AS and the quark orbital angular momentum 

Jg = ^AS + L„ (73) 

where AS is measured through the polarized DIS process. Substituting Eq. (|7^ into 
Eq. (j73|) we can determine Lg. Moreover, the total spin of the nucleon comes from 
quarks and gluons, 

I = J, + J,. (74) 
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where Jg is the total angular momentum carried by gluons. Thus we can also extract 
the gluon contribution to the nucleon spin. Note that AS/2 accounts only for ap- 
proximately 30 % of the nucleon spin. In summary, by measuring GPDs we obtain 
information about the angular momentum distributions of quarks and gluons in the 
nucleon. Therefore, the DVCS process, combined with measurements of the quark 
helicity distributions from inclusive deeply inelastic scattering, can unravel the orbital 
angular momentum carried by partons, on which little or no information is currently 
available. 

III. 5. 2 Double Distributions 

Two approaches are used to model generalized parton distributions: 

• A direct calculation of GPDs in specific dynamical models, such as the bag 
model j^, the chiral soliton model [73], the light-cone formalism [321, ^^c. 

• A phenomenological construction [7^ 17^1 177] based on reduction formulas re- 
lating GPDs to PDFs /at (x) and A/^v (x) and form factors Fi {t), F2 [t), (t) 
and gp{t). 

The most convenient way to construct models in the second approach, and further 
study the interplay between x, ^ and t dependencies of GPDs is performed using the 
formalism of double distributions (DDs) f{l3,a,t). Here we only consider the so- 
called a-DDs corresponding to the symmetric description with respect to the initial 
and final hadron momenta. In the parton picture of double distributions shown in Fig. 
ini the active parton momentum, = (3p^ + (1 + a) r"^/2, is represented as the sum 
of two components (3p^ and (1 + a)r^ /2. The former specifies the momentum flow 
in the s-channel due to the plus component of the average hadron momentum p"*", and 
the latter specifies the momentum flow in the t-channel due to the plus component 
of the momentum transfer r"*". Despite the proportionality between r"*" and , they 
correspond to the momentum fluxes in two different channels. Their superposition is 
the main feature of double distributions. In addition, it is important to note another 
characteristic feature, i.e. the absence of the ^-dependence in / 

Thus double distributions are hybrids, which look like usual parton distributions 
with respect to the variable /?, and like distribution amplitudes with respect to a. 
Therefore, when modeling DDs we usually represent a double distribution in the 
factorized form as the product of a usual PDF in the /3-direction and a distribution 
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^p+(l + a)r/2 



p+r/2 




^p-(l-a)r/2 



p-r/2 



FIG. 9: Description in terms of double distributions. 




P 




a 



FIG. 10: Conversion from a double distribution to an off-forward parton distribution 
function with a zero and nonzero skewness. 



amplitude, which drives the a-profile. The connection between the DD variables 
a and (3 and the OFPD variables x and ^ can be established through the formula 
r+ = 2^^"*", namely, x = (3 + ^a. 

The support area of double distributions is the rhombus |a| + \(3\ < 1. Moreover, 
due to the hermiticity properties and time-reversal invariance of the nonforward ma- 
trix elements, they are even functions in a, f {P,a,t) = f —a,t). This property 
is termed as the Munich symmetry [ZH!- Now, in order to obtain the usual PDFs 
from DDs (recall that the former correspond to the forward limit, ^ = and t = 0, 
of OFPDs), one should simply integrate them along the vertical lines (3 = x, see the 
rhombus (a) in Fig. ^| On the other hand, to get OFPDs H (x, ^, t) (we skip the 
flavor index for convenience) with a nonzero skewness (e.g. for ^ = 1, they behave 
like meson distribution amplitudes), one should integrate DDs f {P,a,t) along the 
parallel lines a = {x — j3) /C, with a ^-dependent slope. 



rl .1-1/31 

H{x,C,t) = dp daS{p + ^a-x) f(p,a,t) 

J-l 7-1+1/31 



(75) 
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as illustrated in the rhombus (b) in Fig. ^| We call this process of integration 
or scanning, the DD-tomography. To summarize, double distributions are primary 
objects producing SPDs after appropriate integration. 
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CHAPTER IV 
DEEPLY VIRTUAL COMPTON SCATTERING 
IV.l INTRODUCTION 

In recent years, significant effort was made to access GPDs through the measurement 
of hard exclusive leptoproduction processes, such as deeply virtual Compton scat- 
tering or deeply exclusive meson production. The simplest process in this respect is 
DVCS. It can be accessed through the following reaction 

e-(fc) + iV(j9i) ^ e'{k')+N{p2)+l{q2), (76) 

and it is illustrated in Fig. ^2 There are three relevant diagrams. The nucleon blob 
with two photon legs, see the diagram (a), represents the virtual Compton scattering 
amplitude, which will eventually become the subject of our study. This diagram is 
referred to as the DVCS diagram or the Compton contribution. Unfortunately, the 
final real photon can be emitted not only by the nucleon, but also by the electron. 
The latter is presented by the remaining two diagrams. They are referred to as the 
Bethe-Heitler contribution. Here lower part, the nucleon blob, stands for the electro- 
magnetic form factor while the upper part can be exactly calculated in QED. Despite 
this disadvantage (in measuring of the VCA by extracting it from the cross section, 
the pure DVCS process is always in competition with the Bethe-Heitler process), in 
addition to a small cross section, deeply virtual Compton scattering is still regarded 
to be the cleanest tool to access the underlying GPDs. 

At this point, we focus only on the Compton part. Hence an electron (or muon) 
scatters off a nucleon via the exchange (in the leading-order QED) of a space-like 
photon with virtuality qf = {k — k'Y < 0, producing an intact nucleon (with altered 
momentum) and a real photon in the final state. Since we turn a virtual photon 
into a real one there is always a nonzero momentum transfer. At the quark level, 
in leading twist, the electromagnetic current couples to different quark species with 
strength proportional to the squares of the quark charges, selecting specific linear 
combinations of GPDs. Flavor-specific GPDs can be reconstructed by considering 
DVCS from different hadrons (protons and neutrons, for instance), and using isospin 
or flavor symmetry to relate GPDs in the proton to those in the neutron. 
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a be 



FIG. 11: DVCS (a) and Bethe-Heitler (b and c) diagrams contributing to electropro- 
duction of a real photon. 



In the Bjorken regime (recall that in addition to the large invariant mass of the 
photon-nucleon system, {pi + gi)^ —>■ oo, the initial photon is also highly virtual, 
nevertheless, their ratio xb = — ?i/2 (pi ■ gi) is finite), the VGA is dominated by 
light-like distances. The dominant light-cone singularities, which generate the leading 
power contributions in l/|q'i| to the amplitude, are represented by two handbag 
diagrams shown in Fig. in which the (hard) quark propagator is convoluted with 
the soft function parametrized in terms of GPDs. In addition, keeping the momentum 
transfer squared to the nucleon, t = {pi — ^2)^, as small as possible, one arrives at the 
DVGS kinematics, s > —qf ^ — t. This particular kinematics implies, on one hand, 
that —qf should be large enough to ensure scaling regime for the amplitude and, on 
the other hand, it implies small t. 

One of the methods to study the virtual Gompton scattering amplitude in the 
DVGS kinematics is the approach based on the nonlocal light-cone expansion of the 
product of currents in QGD string operators in coordinate space [TH]. It will be 
employed in the present work. 

A detailed derivation of the leading-twist (and to the lowest order in a^) VGA 
is provided in Section IV.2. In Section IV.3, we discuss the kinematics, common to 
all DVGS-like reactions. After introducing the simple model for nucleon GPDs, see 
Section IV.4, both the Gompton and Bethe-Heitler cross sections are estimated in 
Section IV.5. 

IV.2 VIRTUAL GOMPTON SCATTERING AMPLITUDE 

We begin with an analysis of some general aspects of the VGA. In the most general 
nonforward case, the amplitude is given by a Fourier transform of the correlation 
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FIG. 12: Handbag diagrams {s- and u-channel) in the virtual Compton scattering 
amplitude. 

function of two electroweak currents. In particular, for the standard virtual Compton 
process on the nucleon, where both currents are electromagnetic one has |543 

(77) 

Due to the current conservation, the amplitude is transverse with respect to the 
incoming virtual and outgoing real photon momenta, 

TemQiu = and g2/.T^M = 0- (78) 

It will be convenient in the analysis to use symmetric coordinates, defined by intro- 
ducing center and relative coordinates of the points x and y, X = (x + y) /2 and 
z = y ~ X. Accordingly, the amplitude takes the form 

T^l^ = t j (fx j (fz e-^('?i-<?2)-^+i('?i+92)-/2 

X {N (p,, S2) I T { J^,, (X + z/2) r^,, {X - z/2)} \N (pi, s,)) . (79) 

Furthermore, in order to treat the initial and final nucleons in a symmetric manner, 
we introduce, as independent momentum variables, the averages of the photon and 
nucleon momenta, q = [qi + ^2) /2 and p = [pi + P2) /2, and the overall momentum 
transfer to the nucleon, t = pi — ^2 = ?2 ~ Then we can write 

g2 = g2^2 - t/A with t = (80) 

From the on- mass-shell condition, pf = P2 = M'^, it follows that 

p"^ = - t/A and {p ■ r) = 0, (81) 

where M denotes the mass of the nucleon. After performing the translation in x in 
Eq. (fTUj). (P2I {x) \pi) = (P2I (0) \pi) e~**^Pi~^2)-a;^ g^^^ integrating over the center 
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coordinate, one finds 

= {2n)U^'Hpi + qi-P2-q2)T>^lf, (82) 



where 



■'■EM 



; j d'z e'('^-) (AT {p - r/2, s,) \ T { J^,, (^/2) J^,, {-z/2)} \N {p + r/2, .0) • 

(83) 



We call the amplitude (j83|) the reduced virtual Compton scattering amplitude. It ap- 
pears in the invariant matrix element (or alternatively the T-matrix) of the standard 
electromagnetic DVCS process. 

Having introduced the reduced VGA we demonstrate its calculation at the twist-2 
level in the DVCS kinematics, which amounts to neglecting contributions of the order 
M^/q"^ and t/q"^ in the hard part of the amplitude, and keeping the t-dependence only 
in the soft part. It was noticed that a straightforward use of the twist-2 result for 
t leads to inconsistencies, e.g. the VGA is not electromagnetic gauge invariant 
The invariance is restored through inclusion of the twist-3 corrections to the 
amplitude. They are power suppressed in g^, and have been calculated by several 
groups [TSUHOllHIllHlEllHlllHnilHnilHZllHHlEni using different approaches. 

Since the final state photon is on shell, g| = 0, it follows in this particular kine- 
matics that {r ■ qi) ^ — Q'i/2 = Xb {pi ■ qi)- Hence the momentum transfer should 
have a large component in the direction of the average nucleon momentum, 

r = 2r]p + A, (84) 

characterized by the skewness parameter. 

In the DVGS kinematics, the remainder A in Eq. (jH^ is transverse to both p and 
g, (A-p) = —2rjp'^ and (A ■ g) = — 1/4 [82j. In the language of the 
parton model, the component 2rip of the momentum transfer would be identified 
with the longitudinal component ry while A would be a transverse component r± 
. Moreover, rj coincides with another scaling variable ^ defined as 
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For that reason, ^ in DVCS is frequently also referred to as skewness. It is easy to 
verify that C, = a;^/ (2 — x^) in the limit t/q"^ —>■ 0. In DIS, on the other hand, the 
asymptotic region, which is also characterized by the large virtuality of the incoming 
photon, corresponds to the situation in which qi = q2, Pi = P2, r = and accordingly, 
1] = and ^ = xb- 

Let us now turn to the formal light-cone expansion of the time-ordered prod- 
uct T {Jj^i^j {z /2) J^^j {—z 12)} in the coordinate representation. The expansion is 
performed in terms of QCD string operators, as discussed in Ref. The string 

operators have gauge links along the straight line between the fields (the gauge link 
connecting the two space-time points corresponds to the summation over twist-0 lon- 
gitudinal gluons, and it disappears, e.g. in the Fock-Schwinger gauge, z^A^ {z) = 0), 
however, for brevity we will not write them explicitly. The leading light-cone sin- 
gularity, z"^ —>■ 0, is given by the sum of two {s- and n-channel) handbag diagrams 
shown in Fig. ^1 The hard part of each of the diagrams begins at zeroth order in 

with the purely tree level diagrams, in which the virtual and real photons interact 
with the (massless) quarks. We have 

tT{r^,,{z/2)r^,,i-z/2)} = iY^Q){^s{zl2)^H ^{z)^^^^{-zl2) 

f 

{~z/2) ^ {-z) 7^^/ {zl2)\ , (87) 

where Qf denotes the electric charge of the quark with flavor / in units of |e|. Note 
that the vertices, in fact, contribute the factor (— z|e|Q/)^ rather than Qj. For 
convenience, we do not include an extra factor of — in the expression for the VGA, 
however, it will be included in the T-matrix. The free quark propagator between the 
initial and final quark fields in the coordinate representation is given by 



' ' 2^T^z^-i0f J (27r)^ + iO' ^ ' 

Using the 7-matrix formula, 

^i^YY = (s'^^^'' + ^£'^^^''75) 7^, (89) 

where s^^'^^ = gt^p g''^ -\- gi^^ gP'' — g'^'^ gP"^ is the symmetric and e^^*^^ is the antisymmetric 
tensor in the Lorentz indices fi and u, we express the original bilocal quark operators 
in Eq. (jHTjl in terms of operators with only one Lorentz index, 

iT{J>^j^iz/2)J'^j,,i-z/2)} = 
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27r2; 



'^^ {-z/2) 7,75^/ (z/2) + ijf (z/2) 7,75V>; (-z/2)] } . 



(90) 



It is customary to write these new QCD bilocal operators (the vector and axial vector 
string operators) as 



<(^|0) 



{-z/2) 7,^; {z/2) - ijj {z/2) {-z/2) 
{-z/2) 7^757/-/ {z/2) + {z/2) 7^75^-/ {-z/2) 



(91) 



and accordingly, Eq. ()9()p turns into 



iT{J%,{z/2)J_ 



u 

EM 



^/2)} 



^ E Q) {s^'-^'^Of^ (^ |0) + le^'^'^^Ol^ {z\Q)}. 



(92) 



The string operators in Eq. (^1]) do not have a definite twist. To isolate 
their twist-2 part one uses a Taylor series expansion of (-2 |0) and (-2 |0) in 
the relative coordinate z. This gives local operators, ■?/'/ (0) 7^D^^...D^^'?/'/ (0) and 
V^/ (0) 7»775-^Mi----^Mn'^/ where is the covariant derivative, which are not sym- 
metric in their indices. To get the twist-2 contribution, one should keep only the 
totally symmetric traceless parts of the coefficients in the expansion. As it was shown 
in Ref. [7H], the totally symmetric parts can be carried out by the following operation 



Oi{zK 



Oi,{z\Q] 



sym 



sym 



d„ 



d/3 

{-Pz/2) {Pz/2) - i,j {(5z/2) {-(5z/2) 
^ dp 

{-Pz/2) /75^/ {Pz/2) + i,j {(3z/2) /75^/ {-Pz/2) 



(93) 



where 9^ = d/dz^ is the derivative with respect to the relative coordinate. It becomes 
clear now why the term string is used: the argument oiipf and ipf takes all the values 
on the string from —z/2 to z/2. The subtraction of traces can be achieved by imposing 
the harmonic condition on the contracted vector and axial vector string operators. 



Of+{z) = z'^0{{z\'d 
{z) = z"0, 



i,f {-z/2) {z/2) + {z^ -z) 
ijf {-z/2) /75^/ {z/2) -{z^- 



(94) 
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which appear on the right-hand-side of Eq. (jHSJ- In other words, these two operators 
should satisfy d'Alembert equation with respect to z 



O^+iz)] = 0, (95) 



and similarly for the twist-2 part of (2;). Note that we have assigned an extra 
superscript ± to the operators (jMj) because they possess the symmetry with respect 
to the change z — > —z. 

To compute the amplitude (|HHj). the contracted twist-2 operators have to be sand- 
wiched between the initial and final nucleon states, and ultimately integrated over 
z. Thus we need to construct a parametrization for these nonforward nucleon ma- 
trix elements. The most convenient way is a decomposition into plane waves, i.e. a 
spectral representation, where the relevant spectral functions correspond to GPDs (to 
be specific, here and in the following, we use the off- forward parton distributions). 
In principle, we need to provide a parametrization valid everywhere in z since the 
coordinate z runs over the whole four-dimensional space. However, it turns out that 
the inclusion of the z^ terms in the matrix elements generate M'^ /q^ and t/q^ correc- 
tions to the amplitude (such corrections are analogous to the well-known target mass 
corrections in DIS [^011^]) and hence will be neglected. It is, therefore, sufficient to 
provide a parametrization only on the light-cone ^3]; namely. 



(iV(p2,S2)|0^+(2:)|iV(pi,si))^2=o = m(P2,S2) A(pi,si) 



(^^2, S2) — U {pi, Si, 



X 



{N{p2,S2)\0[ {z)\N {pi,Si))^2=Q = U{P2,S2) ^J5U{P1,S 



X y^^rfxe"(P")iJ/(a;,e,t) 

-U {P2,S2)^J^15U{P1,SI) 



X 



J^^ dx e^^(P-^)E+ (x, t) . (96) 



The flavor dependent OFPDs in the parametrization refer to the corresponding 
quark flavor / in the nucleon N. Apart from the scale = — g^, there are three 
variables necessary to specify OFPDs, namely, the usual light-cone momentum frac- 
tion X, the invariant momentum transfer t to the target and the skewness parameter 
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^, which corresponds to the hght-cone momentum fraction transferred to the target 
and characterizes the momentum asymmetry. Recall that the variables x and ^ only 
define the longitudinal momenta of the partons involved, that is their plus compo- 
nents. Schematically, the parton going out of the parent nucleon in Fig. carries 
the fraction {x + ^) of the average nucleon momentum p while the momentum of the 
returning parton is {x — ^) p. 

The OFPDs Hf (x, ^, t), Ej (x, ^, t), Hf (x, C,, t) and Ef (x, ^, t), introduced in Sec- 
tion III.5, parametrize the matrix elements of operators ipf {—z/2) ^ipj {z/2) and 
ipf {—z/2) ^'^^ipf {z/2), do not have symmetry with respect to the change x — x. 
In particular, the function Hf (x) in the forward limit, see Eq. (j69|) . corresponds 
for positive x to the quark distribution while for negative x, it corresponds to the 
minus antiquark distribution. On the other hand, the plus distributions, introduced 
through the matrix elements of (anti) symmetrized operators given by Eq. (jMj) . do 
have well defined symmetry properties with respect the scaling variable x. Simply by 
transforming z —>■ —z and x —>■ —x in Eq. (PU)) . we can establish the following 



It is easy to notice that the plus OFPDs are determined by the sum of quark and 
antiquark distributions. Since the quark distribution includes both the valence and 
sea contributions, and further the antiquark distribution only the sea contribution, 
the plus distribution turns into the sum of the valence and twice the sea quark dis- 
tributions. 

After substitution of Eq. ()96|) into the right-hand side of Eq. fl93|) . we first take 
the derivative with respect to z. Then we perform integration by parts over the 
parameter (3 and keep only the surface terms with the arguments ipf (±z/2) and 
ipf (±z/2). Finally, the integral over z is carried out with the help of the inversion 
formula for ^ {z) , 



Here the momentum / is given by / = {xp + g), so that P in the denominator of Eq. 



Hj (x) 
Ej (x) 
(x) 
{x) 



-HI (-X) , 
-Ef (-x) , 
Hi (-X) , 
El (-X) . 



(97) 




(98) 
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becomes P = 2 {p ■ q) {x — C,)- The expression for the reduced VGA reads 

1 ^ dx 

f 



-1 (x-^ + tO) 

X ^Hf {X, ^, t) [^M (p2, Ss) Yu (pi, Si) 
+rM (p2, S2) (pi, Si) - gf^^u {p2, S2) /u {pi, Si)] 

+E+(x,e,t) 



Fu {p2, S2) — u ipi, Si) 



AM 

+ru {p2, S2) 7l_tlplu {pi, Si) - g^'-u {p2, S2) ^^^rj^^ n [pi, si) 



AM 



AM 



-le 



Hj (X, t) IpU {p2, S2) (Pl, Si 

- T 
-Ef {X, t) IpU (p2, S2) 7^l5U (pi, Si 



(99) 



and can be further simphfied as follows. Let us perform the light-cone decomposition 
of 7-matrix, 



r 



The four-vectors rii and n2 in Eq. pOO|) are light-like, n1 



(100) 

?T,2 = 0, and satisfy the 
condition [rii ■77,2) = 1. Identifying rii — > p, ^2 — > q2 and neglecting the transverse 
component •j'^, since it corresponds to the higher- twist contributions, Eq. (|100|) takes 
the form 

1 



r 



(101) 



(P • 12) 

Using the above decomposition, the Dirac equation, j6i'u(pi,si) =^2u{p2,S2) = 
(recall that we neglect the nucleon mass), and the fact that and are odd 
functions in x, the reduced VGA in the leading twist transforms into 



rj-pu _ ^ ^2 

■'■EMtwist-2 ~ o „\ / ^ ^ f 



2{p-q) 



dx 



1 (x - ^ + iO) 



x\Hj{x,i,t) 



1 



+E}{x,i,t) 



{P ■ 92) 

{p'q2 + P^Q^) - 9'" 
1 



{p'q^ + P''q'2) - 9 



flU 



U{P2,S2) 42u{pi,Si) 



(P ■ 92) 

-Hl{x,U) 

1 



{P ■ ^2) 



-I Ai2i-ii2) , . 

U (P2, S2) — U (pi, SiJ 

U{P2,S2) ^2l5U {Pl, Si) 



+El{x,^,t) 



{p ■ 92) 



ie^''''^q2pP^ 



q2 ■ r) 
2M 



u {P2, S2) -f5U {Pl, Si) 



(102) 
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The amplitude has both real and imaginary parts. Since OFPDs are real functions, 
the imaginary part of the amplitude comes only from the singularity of the expression 
1/ [x — ^ + iO). The latter can be calculated by applying the formula 

^ = p(-)^ in5 (x) , (103) 



{x ±iO) \x 

where P denotes the principal value. Thus the real part of Tj^m is obtained using the 
principal value prescription. The imaginary part, on the other hand, generates the 
delta function 5 (x — ^), which assures taking OFPDs at the specific point, x = ^. In 
other words, the imaginary part of the amplitude is directly proportional to OFPDs 
evaluated along the line x = ^. Moreover, from the tensor structure of the VGA in 
Eq. (jl02p . it is easy to see that the amplitude exactly satisfies the electromagnetic 
gauge invariance with respect to the final real photon, 

<l2,,T^Mtwist-2 = 0. (104) 

The gauge invariant condition with respect to the initial virtual photon that has the 
momentum qi^, = — will only be satisfied if 'T^EMtwist-2''^i' = 0- However, for the 
part of '^EMtwist-2^ which is symmetric with respect to /i ^ z/, we find with the help 
of Eqs. dH and dHSl), 

EMtwist-2(symy f ' (^ ■ ^2) 

= A^-C(r2)p^'. (105) 



Since the components of A are all of the order ylt], see Ref. |83;, it follows that 
the leading-twist VGA is gauge invariant to accuracy 0(^\J\t\^. This violation of 
electromagnetic gauge invariance is a higher-twist (i.e. twist-3) level effect. 



IV.3 KINEMATICS 

In the generalized DVGS process, 

h{k)+N{p^) l^{k') + N'{p2)+i{q2), (106) 

there are a lepton li with the four-momentum k and a nucleon (pi) in the initial 
state, and a lepton I2 {k'), a nucleon A^' (^2) and a real photon 7 (^2) in a final state. 
The process has two contributions that are represented by two types of diagrams, see 
Fig. ^2 The DVGS diagram corresponds to the emission of the real photon from 
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the nucleon blob. In the Bethe-Heitler diagram, the real photon is emitted from a 
lepton leg. We denote the four-momenta as = (^uj,kj, pi = (^M,Oj, k' = (^uj',k'^, 
P2 = (-^2,^2) and q2 = (1^2,^2)- The differential cross section for lepton scattering off 
a nucleon, to produce a final state with a lepton, nucleon and a real photon is 

= 2;"^" (^^ ^^ + ^^-^-^^-^^^^2^2^' 

where T represents the scattering amplitude. It contains both the Compton and 
Bethe-Heitler contributions, 

T = Tc + Tbh. (108) 

In the laboratory frame, as the target rest frame, the total lepton-nucleon center-of- 
mass energy squared is s = {k + pi)"^ = 2ujM + M^. Note that in the following, we 
neglect both lepton masses. Integrating Eq. p07|) over the photon momentum gives 



1 



da = — ITI^ r 6 



UJ' 



'du'dn' d^p2 



{k + pr-k'- p2Y\ ^ (109) 



2s'"' (27r^^ 

The delta function in the cross section (jl(J9|) provides the constraint s + — 
2 [(yi + M) E2 — qi ■ P2] = 0, where the invariant s = {pi + qif' and qi = k — k' = 
(z/i,gi) is the four-momentum of the virtual photon. Furthermore, we choose the 
coordinate system so that the z-axis is in the direction of the incident lepton in the 
plane formed by the lepton momenta. After integrating Eq. (jl(J9|) over the magnitude 
of P2, the differential cross section reads 

|2 

16s (27r)^ \{Ul + M) \p2\ - |gi|^2COS0i2 



da = ^ .^M^i ill 7:. „„„ ^ ^ \Tfdu;'dn'dn2, (HO) 



where 0i2 is the angle between vectors gl and p2. 

Alternatively, one can choose the coordinate system depicted in Fig. El in which 
the virtual photon four-momentum has no transverse components, qi = {ui, 0, 0, \qi\), 
and the incoming and outgoing lepton four-momenta are k = u {1, sin0, 0, cos0) and 
k' = cij' (1, sin0', 0, COS0'), respectively. In this reference frame, the azimuthal angle 
of the recoiled nucleon corresponds to the angle between the lepton and nucleon 
scattering planes. Using now the delta function 6 {k + pi — k' — ^2)^ in Eq- (!109|) 
to integrate over the polar angle 02 of the outgoing nucleon, we find 



2JEI - J{u - ui'Y - ql 
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FIG. 13: Kinematics of the generalized DVCS process in the target rest frame. 



and the cross section assumes the form 



da 



-duj'dE2d (cos (f)') dif, 



(112) 



with the angle 0' denoting the polar angle of the scattered lepton. 

Instead of the kinematical variables cu', E2 and cos0', it is convenient to express 
the differential cross section in terms of the invariant variables, e.g. {y,t,XB) or 
{Ql, t, xb), where Q\ = —q\ and y = {pi ■ qi) / (pi • k). In the laboratory frame, these 
invariants are given by 



Ql = 2Mu;yxB, 



y 

t 



UJ — UJ 
UJ 

= 2M'^-2ME2. 



(113) 



The last invariant can also be written as t 



{q,-q2f = -Ql 



2v2 \ ijjy — \Joj'^y'^ + Q\ cos j , where the energy of the outgoing real photon is given 
by the energy conservation, 1/2 — M + coy — E2. By combining both expressions for 
t, we can express it as a function of the angle 9^^ between the incoming virtual and 
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outgoing real photon, namely, 

2M2 + Ql + 2{M + ujy) {ujy - ^u'^y^ + Qj cos 6^^^ 



t =2M2-M 



M + u;y- ^Ju'^y'^ + Ql cosO^^ 

(114) 



Moreover, from the three-momentum conservation, qi = k — k', and alternative 

expression for the invariant Ql in the laboratory frame, Qi = ~ {k — k')^ = 
2uiu}' (1 — cos 4>iep)i we get for the polar angles of the incoming and scattered leptons. 



cos 



and 



<t> = ^(7^^^vTgf + a;(l-y)cos(/)') (115) 



COS0 = , =, (lloj 

2uj (1 - y) ^u^y^ + Ql 



respectively The Jacobian is then equal to 



J = 



d (uj\ i?2, cos 6') 



uj\/ {uj + Mxb) 



2{l-y) (^2^2 + Q2) 



d{y,t,XB) 

and accordingly, the differential cross section turns into 



2\3/2 



(117) 



J 



dxBdydtd(p d (cos 0') duj'dE2dip 

M/f.A 

T\\ (118) 



1 1 l+XB(M/a;) ,^,2 



32s (27r)^ [y + 2xb (M/u;)]' 
Note that using the invariant Ql instead of y yields 

rfV 1 1 1 + xb{M/uj 



^|^ (119) 



dxBdQldtdcp 64s (27r)^ Mcuxb [y + 2xb iM/uj)f 

In the following, we determine the kinematically allowed region for the generalized 
DVCS process. In other words, one needs to find the upper and lower limits on xb, 
y and t. We require the following constraints: 

• The energy of the incoming lepton beam is fixed. In particular, wc consider two 
examples that are both relevant to Jefferson Lab, namely, a; = 5.75 GeV and 
a; = 11 GeV. 

• The invariant mass of the virtual photon-nucleon system should be above the 
resonance region, s = (pi + q^f > Smin — 4 GeV^. 
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• The virtuality of the incoming photon has to be large enough to secure the 
hght-cone dominance, Ql > Q\min = 2.5 GeV^. 

• The momentum transfer squared to the nucleon should be kept as small as 
possible, t/q^ <^ 1. 

The second constraint, M"^ — Q\ (1 — l/x^) > Smim together with Q\ = xbV {s — M^) 
implies that 

[s - M^) (1 - Xb) 

On the other hand, the scaling variable y reaches its maximum value, 

= (l + J3l|ij . (121) 

when the incoming lepton is aligned along the 2;-axis at the angle 0;ep = 180°. If one 
plots the region in the xbV plane, as illustrated in Fig. [T31 then Eqs. p20|) . p2H) 
and y = Qlmin/^B (s — M^) (the latter comes from the third constraint) correspond 
to its boundaries. Next, both lower and upper limits of the invariant t can be found 
if we go in the virtual photon-nucleon center-of-mass frame. In this particular frame, 
at the scattering angles between the initial and final nucleons equal to 0° and 180'', 
the invariant momentum attains, up to relative corrections of the order XbM"^ /Qf, 
its kinematical limits, 

-M'xl 



_ i-^-i J.B -^B -r v^i yu^B - ±j /^B (^nn\ 

T^max — 1 „ /I i\/r2 //n2\ ' \ ) 



l-XB{l-MyQj)' 
Wxl - 2M^xb + Ql {xb - 1) /xb 

i-xb{i- MyQi: 

They are presented, within the kinematically allowed regions, in Table HI Since we 
require small t, the upper limit becomes irrelevant. In addition, having small t yields 
a low-energy nucleon with E2 = M {1 — t/2M^) and a high-energy real photon in the 
final state. 

Finally, for each lepton energy we pick one kinematical point within the allowed 
region in the XBy plane. In particular, we take = 2.5 GeV^ for the virtuality of 
the initial photon and xb = 0.35. In Fig. ^[ the invariant t is plotted against the 



angle 0^^. It is customary to present the plot for the values of t up to —1 GeV , which 
happens at 6^^ ~ 12°, even though we require that, in principle, —t should be much 
smaller than Ql. The values of t vary from -0.15 GeV^ to -1.433 GeV^ for angles 
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TABLE I: Lower and upper limits of the invariant momentum transfer t within the 
kinematically allowed region illustrated in Fig. El for two different lepton beam 
energies u. 



u [GeV] 



GeV^ 



GeV^ 



5.75 
11 



0.058 
0.014 



10.019 
19.853 
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t [GeV^] 



Oyy [d^] 



2 



4 



6 



8 



10 



12 



14 



-1.4 



-0.6 



-0.2 



-0.4 



-0.8 



-1.2 



-1 




FIG. 15: Invariant momentum transfer t plotted as a function of the angle between 
the incoming virtual and outgoing real photon in the target rest frame for Q\ — 
2.5 GeV^ and xb = 0.35. 



TABLE II: Polar angles (j) and (f)' of the incoming and scattered leptons, respectively, 
in the target rest frame for Ql — 2.5 GeV^ and xb — 0.35 with two different lepton 
beam energies cu. 



uiGeY] 0[deg.] 0' [deg.] 



5.75 12.6 39.9 

11 18 28.2 
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between < 0^^ < 15°. Note that by fixing Q\ and x^, the product uy = QI/2Mxb, 
which appears in Eq. ()114|1 is also fixed and accordingly, we end up with the same 
6'^^-dependence for both beam energies. For convenience, we set the angle between 
the lepton and nucleon scattering planes to = 0, and arrive at the so-called in-plane 
kinematics. The polar angles of both incoming and scattered leptons are then fixed 
and given in Table ITTl 

IV.4 TOY MODEL 

Our simple model has the following properties: 

• We assume that the sea quark contribution is negligible. For that reason, the 
plus OFPDs in the parametrization (j^Hl) are equal to the valence OFPDs, H^^ = 

and = with the quark fiavor f = u, d, and similarly for E'j and 
distributions. 

• We assume that for all distributions the t-dependence factorize from the depen- 
dence on other two scaling variables. The dependence of GPDs on the invariant 
t is then characterized by the corresponding form factors. 

• The ^-dependence of OFPDs appears only in the E'j distribution. 

The parametrization of the unpolarized quark OFPDs is taken from Ref. |92j . 
Namely, 



In particular, for the proton target the unpolarized valence quark distributions in the 
proton are given by |^ 



ut (x) Fi„ it) /2, 

Ut {X) F2u (t) /2, 

< {x)F2d it). 



(123) 




1.89a;-°-^ (1 - x)^-^ (1 + ex) , 
0.54x-°-^ (1 - xf-^ (1 + 8x) . 



(124) 



They closely reproduce corresponding GRV parametrizations I^Sl at a low normaliza- 
tion point, -ql ~ 1 GeV^ 
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The u- and c?-quark form factors in Eq. (jl23|) can be extracted from the proton 
and neutron Dirac and Pauh form factors according to Fip(„) = QuFiu(d) + QdFid(u) 
and F2p(^n) = QuF2u{d) + QdF2d{u)- Furthermore, the proton and neutron form factors 
are related to the Sachs electric and magnetic form factors, see Eqs. and 
through 



^ lp(n) 



GEp(n) (t) — 



t 



F2p{n) (t) 



GMp{n) (t) — GEp(n) (t) 



GMp(n) (t) 
1 - 



t 



t 



4M2 

1 



-1 



4M2 



(125) 



where the nucleon mass is M ~ 0.94 GeV. In the region of small t, both Sachs form 
factors are well described by a dipole fit. 



Gep (t) 



Gmp (t) Gmu (t) 



1 + 



1 - 



t 

A2 



and Geu {t) = 0, 



(126) 



with the parameter = 0.71 GeV^ and Kp = 1.793 and «;„ = —1.913 are the proton 
and neutron anomalous magnetic moments, respectively. In the polarized case, we 
take for the valence distributions jHS] 



Hr{x,^,t) = Ad 



val 



t 



(127) 



with the mass parameter = 1.03 GeV. The t-dependence in Eq. (|127j) corresponds 
to the ratio qa (t) /g^it = O)- The polarized valence quark distributions in the proton 
can be expressed in terms of the unpolarized distributions in the following way |^ 

2 



Au;''' 



cos^^ (m^'^' 



3> 



Ad. 



val 



COS 6 



D 



^jval 

"3> 



'128) 



where cosOd = [1 + Hg (1 — x^) / y/x ] ^ with Hq = 0.06. Finally, for the Ef distribu- 
tion we accept the pion pole dominated ansatz, 

fx + ^\ 



FT {x,^,t) 



ipval 

Fd 



2e 



The function (t) is taken in the form valid for —t <^ ^7] , 

FAt) = 4gA{t = 0)M' ^ 



[nit 



t) /GeV2 (1 - t/2 GeV^y 



(129) 



(130) 
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where m^r — 0.14 GeV denotes the pion mass and qa (t = 0) = 1.267. For the pion 
distribution amphtude in Eq. ()129|1 we choose, for simphcity, its asymptotic form 

= Qu{l-u). (131) 

IV.5 CROSS SECTION 

Having defined the kinematics and described the simple model, in the following we 
compute separately the unpolarized cross sections for the Compton and Bethe-Heitler 
contributions to the DVCS process on a proton target using an electron beam. In 
fact, since the scattering amplitude consists of two parts, see Eq. fll08|) . its modulus 
squared (and, accordingly the cross section) is given by the sum of three terms, 
namely, the Compton term, the Bethe-Heitler term and the interference term. Here 
we only consider the first two terms. Their cross sections are plotted for the in-plane 
kinematics against the angle 6^^^ by taking the same kinematical point, Q\ = 2.5 GeV^ 
and xb = 0.35, however, for two beam energies, u = 5.75 GeV and 11 GeV. 

With the help of the momentum-space Feynman rules for QED we can immediately 
write down the T-matrix for the pure Compton process. Denoting the polarization 
vector of the final real photon by e* (^2), we have 

zTc = «(A:')0|e|7>(A:)(^)(-e^)e;(g2)(-zr^7/), (132) 

or 



e 



3 

where Tj^m is given by the expression ()102|) . Then one should average the square 
of Eq. ()133p over the initial proton and electron spins, and further sum it over the 
final proton and electron spins and photon polarizations. Note that this particular 
summation is performed using the Feynman gauge prescription, i.e. one can replace 

e*(g2)ea(g2) — > -Qt^a (134) 

7 polar. 

by virtue of the Ward identity. As a result, we get a factorized expression, 

——2 _ {^iraf c Tru/3 ..oc-N 
|lc| - — -4 — ^^/3^c ) (135) 
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in terms of the electron and hadron tensors. Neglecting the electron mass, the electron 
tensor simply reads 



kyk'p + kpk'^ - g^fi {k ■ k') 



(136) 



The hadron tensor requires some algebra, however, if we define the convolution inte- 
grals of OFPDs as a new set of functions. 



-1 (x - ^ + iO) ^ 



dx 



(137) 



and, in addition, neglect terms O (t/qf) and O (M^/qf), then H^^ can be written in 
a compact form as 

1 



H, 



uf3 
C 



2 EM \-^^iEM 



X 



9 



1-e 
t 

v/3 _ 



n- 



+ 



t 



4M2 



{P ■ (I2) 



M2 



4M^ 



'138) 



The initial overall factor of 1/2 in Eq. (|138|) comes from averaging over the initial 
proton spin. Substituting now Eq. p35|) into Eq. ()119|) gives the unpolarized dif- 
ferential cross section for the Compton contribution to the standard electromagnetic 
DVCS process 



d^ac 



a 



1 + xb (M/uj) 



dxBdQfdtd^ IGnQf M^u^ [2 + {M/u)] xb [y + 2xb {M/uj)] 



(139) 



Its dependence on the angle 9^^ is shown in Fig. ^1 

Apart from the Compton part, the Bethe-Heitler contribution can be fully calcu- 
lable in QED with the knowledge on the nucleon form factors. The amplitude, which 
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(fa-c/(dcBd^,ddp) [ifc/GeV*] 



2 4 



10 12 14 



0yy [deg] 



FIG. 16: Compton cross section ac plotted as a function of the angle 6^^ between the 
incoming virtual and outgoing real photon in the target rest frame for Q\ = 2.h GeV^ 
and xb = 0.35 with oj = 5.75 GeV (solid line) and u = 11 GeV (dashed line) electron 
beam. 



in contrast to Tc is purely real, emerges from two Bethe-Heitler diagrams shown in 
Fig. CH 



BH 



u{k') {^\e\r)e:{q2) 



{k' + gs) 



{i \e\ 7^") u (k) 



X {-i |e|) {p iP2, S2)\ Jem (0) \p {pi, si)) 



+u{k') it\e\Y) 



{z\e\r)e*(q2)u{k) 



(140) 



X {-i |e|) {p {p2, S2)\ J^M (0) \P (Pl, Si)) 

Using the relations, u {k') 7^ /(f' = 2u {k') k'^ and fi'-f^u {k) = 2k^u (p), we simplify the 
numerators of both fermion propagators in Eq. ()14()|1 . and obtain for the T-matrix 

I |3 



BH 



t 



< (^2) u {k') 



7^ ^21" + "^k'^Y ^ -Y 42Y + 27"^'' 



2 {k' ■ q^) 



-2 {k ■ ga) 



u{k) 



x{piP2,S2)\J^''{0) \piPuSi)). 



(141) 



The matrix element of the proton transition current, see Section III.2, is parametrized 
in terms of the usual Dirac and Pauli proton form factors. 



{piP2,S2)\JriO)\piPuSi)) = n(p2,S2) 



Fip (t) 7, - (t) 



2M 



u{pi,si 
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(142) 



The spin- averaged square of Eq. ()141|) is again written in terms of two tensors, 



— ^ _ (ATm) p 

\i-BH\ — L^jjJ^ulB ■ 

Here the hadron tensor reduces to a simple expression, namely. 



(143) 



ttBH 



1 



- E (p (P2, .2) I ^ (0) \p (Pl, .l)) {p {P2, S,) I Jj^ (0) \p {pu S,)r 



9uP 



[Fi, it) + F,, {t)f 



+4 



Piu + 



Pip + 



rp 



(144) 



and the lepton tensor is 



BH 



in. 

2 



X 



2 {k' ■ q^) ^ 2 {k ■ q^) 



2 {k' ■ q2) 



2 (fc • ga) 



(A:' ■ g2) 



k'ql + 



9''' {k ■ q2) 



+ 



{k ■ q2) 

{k ■ q2) h^k'^ + kf^k"" + 2k"'k"^ 



k"q^ + k'Pq]^ - g^l" {k' ■ q^) 



(k' ■ q2) {k ■ q2] 
- {k' ■ ga) [k^'k"^ + k^k"" + 2k''kf^' 

+ {k ■ k') [k^'q^ + k'^q'^ - k'^q^ - k"^q!^ + 2k''k'" + 2k'^k'' 
+2g''^{k-k') [{k'-q2)~{k-q2)-{k-k')] . 

Finally, the unpolarized differential cross section reads 
d'^asH 1 + xb{M/uj) 



(145) 



dxBdQldtdif 167rt2 m^uj'^ [2 + (M/u;)] xb [y + 2xb {M/uj)] 



jufi ttBH 
2^BH'^vl3 ■ 



(146) 



The Bethe-Heitler contribution is illustrated on a logarithmic scale in Figs. El and 
OH We plot the latter for a wider range in 9^^ to be able to see both poles at the 
angles (f) and 0', given in Table ITTl The poles corresponds to the situation, when the 
outgoing real photon is coUinear either with the incoming or scattered electron. In 
Fig. E[ a logarithmic plot of both contributions together is shown for = 5.75 GeV. 
The comparison of the results shows that the Bethe-Heitler signal is well above the 
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FIG. 17: Bethe-Heitler cross section asH plotted as a function of the angle 9.y^ between 
the incoming virtual and outgoing real photon in the target rest frame for = 
2.5 GeV^ and xb = 0.35 with u = 5.75 GeV (solid line) and u = 11 GeV (dashed 
line) electron beam. 



Compton one. This is easy to see just by comparing the factorized expressions fll35|) 
and ()143|) . and recalling that t <^ qf. The presence of in the Bethe-Heitler part 
enhances its contribution with respect to the Compton part, which is proportional to 

One way to minimize the contamination with the dominating Bethe-Heitler pro- 
cess is to find the kinematical regions, where the Bethe-Heitler contribution is sup- 
pressed, or at least comparable with the Compton contribution. Another way is to 
exploit the interference between the two processes. This approach is based on the fact 
that Tc has both real and imaginary parts while Tbh is purely real. By incorporating 
the interference terms between the Compton and Bethe-Heitler amplitudes, we can 
disentangle and '^Tc- In particular, using the positron beam in addition to the 
electron beam, one can measure the so-called beam-charge asymmetry, which is sen- 
sitive to the real part of T^^. Furthermore, measuring the single-spin (or alternatively 
the beam-spin asymmetry) by considering electrons with opposite helicities gives ac- 
cess to the imaginary part of Tc. First experiments of this kind were performed at 
Jefferson Lab jHE] and at Hermes [^H]- As a result, we can project out independently 
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FIG. 18: Bethe-Heitler cross section asn plotted as a function of the angle 6^^ between 
the incoming virtual and outgoing real photon in the target rest frame for Ql = 
2.5 GeV^ and xb = 0.35 with u = 5.75 GeV (solid hne) and a; = 11 GeV (dashed 
line) electron beam. 



10 12 14 



FIG. 19: Compton cross section ac (solid line) and Bethe-Heitler cross section asn 
(dashed line) plotted as a function of the angle between the incoming virtual and 
outgoing real photon in the target rest frame for — 2.5 GeV^ and xb — 0.35 with 
cu — 5.75 GeV electron beam. 
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both parts of the amphtude and accordingly, probe different hnear combinations of 
OFPDs. 
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CHAPTER V 

INCLUSIVE PHOTOPRODUCTION OF LEPTON PAIRS 

V.l INTRODUCTION 

By the inclusive photoproduction of lepton pairs, see Ref. |l(JUj . we refer to the 
reaction, in which a high-energy real photon 7 (with the four-momentum q) scatters 
inelastically from a nucleon (P) emitting a pair of leptons (electrons or muons) 
with momenta k and k', 

^iq) + N{P) I- {k) + l+ {k')+X. (147) 

The process is shown in Fig. 1201 where again X labels a system of hadrons produced 
through inelastic processes. The reaction ()147|) is a crossed channel to inclusive virtual 
Compton scattering, 

e' + N — > e-+7 + X, (148) 

which was originally studied in the parton model in Ref. |lUlj . 

In the framework of the QCD parton model introduced in Chapter II, the elemen- 
tary photon-parton scattering subprocess can be viewed in two different ways. The 
incident photon can either scatter off a parton or split into a lepton pair. Thus we 
have two contributions to the process at the amplitude level. According to the first 
scenario, known as the Compton contribution, a heavy time-like photon, 7* {q') with 
q'^ > 0, is produced and decays eventually into a pair of leptons, as shown in Fig. 




FIG. 20: Inclusive photoproduction of lepton pairs. 
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FIG. 21: Compton contribution to the inclusive photoproduction of lepton pairs in 
the parton model. 




(la(p) qAp') 

FIG. 22: Bethe-Heitler contribution to the inclusive photoproduction of lepton pairs 
in the parton model. 

1211 We discuss this particular subprocess (see also |102j ) in Section V.2. The second 
scenario is illustrated illustrated in Fig. |221 and corresponds to the Bethe-Heitler 
mechanism. It will be studied in Section V.3. The real and virtual photons can be 
interchanged, and hence both the Compton and Bethe-Heitler contributions consist 
of two Feynman diagrams. 



V.2 COMPTON SUBPROCESS 

The invariant matrix element for the Compton part comes from two Feynman dia- 
grams shown in Fig. |^ Neglecting the parton masses by taking = p''^ = 0, we 
have 
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X {-iQa\e\i'')eu {q)u{p) 



+u{k) [i\e\'y^jv{k') 
X {-iQa\e\Y)u{p) , 



(p - q') 



(149) 



where e,y (q) with = denotes the polarization vector of the initial real photon. 
By simplifying the numerator of each propagator in Eq. (jl49j) . we obtain for the 
T-matrix 

\efQl 



■c 



12 



(k) jf,v {k') (g) u {p') 



+ 



q 

-2 (p' ■ q) 



Y 4Y + 27>'' 
2 ■ g) 



u{p) 



(150) 



We will assume that both the nucleon target and the photon beam are unpolarized. 
Then the square of Eq. ()150p has to be averaged over the initial parton and photon 
polarizations, and further summed over the final parton and lepton polarizations. It 
can be written in a factorized form in the following way 

-CI = —J^L^H^c'- (151) 



— ^ (47ra)Q^ ^ 



The lepton tensor is 



4 



kf^k + /Cp/c g^p {^k ■ k ) -\- m 



(152) 



where m denotes the lepton mass. The hadron tensor, on the other hand, is more 
complicated. In the Feynman gauge, it reads 



Ac- 
Be 



Ac Be 



a 



c 



Xp-qf ' {p'-qf ' {p-q){p'-q). 

{p-q)[p'^q''+p''q^-g^'{p'-q)], 

{p ■ q') [p'^'pp + p V] + {p ■ q) [p^q^ + p^q'^ - '^p'^p"] 



/2 



c 



+ ij> ■ p) I'^p^p^ - p^q" - p^q^\ + Y [-pV - + 9^^{p ■ p)\ 

-f' {p ■ q') ip' ■ q') , 

[p ■ q) [p'^'q'f + p'Pq'^" - p'^'pP - p'Pp^"] + [p' ■ q) [2p>^ - p^'q'f - p^q'^"] 

+ {p ■ p) [-p^q" - p^q" - pW^ - p^q'^] 

+g^P [2 {p ■ p') {p ■ q') + (g ■ q) {p ■ p') + {p ■ q) {p' ■ q) - {p ■ q) {p' ■ q')] ■ 

(153) 
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At this point, it is convenient to introduce the Mandelstam variables for the scattering 
subprocess at the parton level, namely, 



(p + q) = {p' + q') . 
(p' - pf = (q' - qf 



S = 

i = 

u = {q' -pf = [p' -qf . (154) 

The four-momentum conservation implies that s + t + u = q'"^ . With the help of Eq. 
fjl54p one can immediately verify the gauge invariance of the hadron tensor, 

J2 a; , ^I2i , „,2„-. (g'2 _ ^) Ut + q'H) - q'H 



n' n' H^P 



sq 



si + q'H + q'% 



su is 



g'2) + (g'2 



+ 



u) — ut {q — u) 



su 



-t-t + 2t 



(155) 



In addition, by setting g'^ = in Eq. ()153p . in other words, the outgoing photon now 
becomes real, and also replacing the lepton tensor with the metric tensor —g^p, we 
easily recover the high-energy limit of the usual Compton scattering process. 



■c| 



(156) 



To calculate the unpolarized scattering cross section, we should integrate Eq. 
fllSip over the Lorentz-invariant phase space defined as 



(27r)^(5(^) {p + q-p' -k-k') 



d^p' 



{2TTf2Ep> 



d^k 



(27r)32cj 



d'^k' 



(27r)32w' 



(157) 



in a specific frame of reference, and divide the result of integration by the flux factor. 
The latter is written in the invariant form as 2s. Since the brute-force contraction of 
tensors and H^f is somewhat tedious, we use a little trick instead. One notices 
that L^p depends only on the momenta of the final leptons. For that reason, we can 
first integrate the lepton tensor over k and k', then contract it with H^'', and finally 
perform integration over the remaining momentum p'. Switching to the symmetric 
momentum variables (i.e. the lepton pair four-momentum, q' = k+k', and the relative 
four- momentum, k = {k — k') /2) and, with the help of the formula 

d^k 
2a7 



J d^k 5+ - 



m 



(158) 



61 



where {k"^ — m^) = 5 {k"^ — m?) 9 [ko) ^ writing the three-dimensional Lorentz- 
invariant volume in the momentum space as the four- dimensional integral, we find 

d?k rd^k' 



d\' ip + q-p'-q')\ I d'^K (5+ \{q'/2 + Kf - 



c 
tip 



x6^ 



{q'/2 - nf - 



L 



c 
pp 



(159) 



The lepton tensor is now equal to 



-c 
^pp 



1 o'^ 
2^M^P ■^'^p'^P ~ 9pp 2 



(160) 



The integral over the relative momentum on the right-hand-side of Eq. ()159j) has a 
tensor structure, and it can be constructed out of the tensors g^p and q'^q'p and scalar 
functions of the invariant mass of the lepton pair g'^. Moreover, the Ward identity 
implies that this integral must assume the form 



d^K6^ 



{q'/2 + - m^j 5+ Uq'/2 - - m 



■c 
'pp 



$ [q 



-12 



(i'^9pp - q'pq'p 



(161) 

where the scalar function $ (g'^) is regular (i.e. it has no pole) at g'^ = 0. Next we 
calculate $ (g'^). It is a function of a Lorentz scalar and thus frame independent. As a 
reference frame, it is convenient to choose the center-of-mass frame of the lepton pair. 
In this particular frame, the symmetric momenta are simply given by q' = (^y/q^, Oj 
and K = (0, k). After contracting both sides of Eq. p61|) with g^^^, then using 



dx / dy 5{x)5{y) = 2 dx dy 5{x + y)6{x — y) 



(162) 



and finally integrating over k, (the angular integration simply gives Air), we get for 
the scalar function 



$ (g'^) 



3g'2 
167r 
'3g^ 



d^K 6^ 



{q'/2 + nf - m^' 



6' 



2Jq'^ Ko 



{q'/2 - nf - m^l [3q'^/2 + 2k^ 
q'^/2 + 2Kl-2\K\^ -2m- 



X 1^1 3g'V2 + 2Ko-2|K 
27r g'2 + 2m' 



q 



12 



2 



q 



12 



(163) 
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Note that the scalar function $ (g'^) < 0. For the sake of completeness, it is worth 
showing that generalizing the calculation of $ (g'^) < to any reference frame indeed 
yields the same result. For instance, integration over hq gives 

2 



~ J d{cos^) J d{\Kf) \k\6 



X 



3q'y2 - 2 \k\^ 1 



q'y2-2\K\'[l 

COS^ '^/Qq 



where is the angle between vectors k and q' in the lepton scattering plane. Further- 

I — *i 2 

more, integrating over \k\ together with the substitution z 
the integral of the type 



dz 



l-l'l'H (1-^2)3/2' 

which leads to the same result for $ (g'^) given by Eq. 
Thus the double integral (jl59|) reduces now to 

d^k fd^k' 



(164) 
fther- 

cost^/^Q produces 
(165) 



/|^/|^5(^Up + g-p'-A:-A;')Lj, = Jd'q'6^'Hp + Q~p'-Q 



X 



/2 / / 



(166) 



Contracting the tensor in the integrand with Hff can be easily carried out. In terms 
of the subprocess Mandelstam variables, we obtain 

2g'2r 



TTlip 



Aq 



12 



U S 
- + - + 

s u su 



(167) 



Note that the second term on the left-hand side of Eq. ()167|) gives no contribution 
due to the condition After combining the result of contraction with Eqs. pSlj) 

and (|157|) and using Eq. (jl58j) . this time for the momentum p', the unpolarized cross 
section for the Compton subprocess reads 



ac = 1 d^q J d^p' 5^^) {p + q - p' - q') 6' 
x<l> 4g'- 



'2 



p- 



(47r«)^Ql 



2s {2Tif 2q 



/4 



u s 2q'H 

- + - + -h- 

s u su 



(168) 



The delta function S^'^^ (p + q — p' — q') can be trivially integrated over p'. The final 
integral over q' is performed in the photon-parton center-of-mass frame. Since we 
want the cross section ac to be differential in the invariant mass of the lepton pair, 
an extra substitution 



/A' = j d'q' j dMl„ 5 {q'^ - Ml^^, 



(169) 
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is being made. Then 



IT 



d (cos Or 



d 



X6 [q',' - 



-Ml 



pair 



1 1 



dqQ 5' 



Ml„ + s-2^/~s 



s (27r)^ 



M2 . 

■^'■^patr 



s u sit 



(170) 



where dcm is the angle between the directions of the incoming real and outgoing virtual 
photons in the photon-parton center-of-mass frame. The center-of-mass energy Ecm 
is equally distributed between the photon and the parton and hence the invariant 



s = -Ecm- differential form, after integrations over and 



we have 



dMp^^^d (cos 6c 



pair 



4s 

\ pair J 



1J_ {AnafQt 
s (27r) 



u s 

- + - 

s u 



M2 . 

pair 

2M2 . t 

pair 
SU 



(171) 



Finally, replacing cos 6', 
(q'o - 



by the invariant i through the relation i 



M2 . 

pair 



Mpair cos 6cmj, whcrc = ("^ + ^pair) /2V^, wc arrive at the par- 



tonic cross section 
d'ac 



1 + 



2m^ 

^'^pair , 



\ 



4m^ 



■^'■^patr ^ ^'^pair 



2M^ ■ t 

pair 



SU 



(172) 



We calculate now the cross section for photon-nucleon inelastic scattering with 
the help of the parton model master formula, see Section II. 3. The process cross 
section is obtained by summing Eq. p72|) over all types a of charged partons and 
all possible longitudinal momentum fractions x, where the single summands must be 
weighted with the proper PDFs fa/N (x). Namely, 



7(g)iV(P) (k) t{k')X 
^ I dx fa/N (x) a 7 (g) (xP) r {k) (A;') qa (p 



(173) 



In addition, the subprocess Mandelstam variables are expressed in terms of the process 
variables. The invariant i is identical to t = {q' — g)^ . Neglecting the nucleon mass, 
the invariant s is simply given by s = xs, where s = {P + q)^ . Accordingly, the 
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Mandelstam invariant u = M^^-^. — t — xs. Recall that by virtue of the vanishing 
mass of the scattered parton, p'"^ = [xP + q — q')^ = 0, the longitudinal momentum 
fraction coincides with the Bjorken scaling variable, 

t 



Substituting Eq. p72p into Eq. ()173|) ultimately yields the Compton differential 
cross section for the inclusive photoproduction of lepton pairs in the parton model 



d^ac {-fN 1-1+ X) 
dMp^i^dtdxB 



^ I dx S{x - Xb) fa 



/iV [x) i — - — 1 + 



M2 . 



\ 



4m^ 



•pair 



[xsf M2 



pair 



[XS 



f + [Ml,,,-t-xs)\2Ml,„t 



xs 



(xs + t- Mj 



2 

pair 



' 1 + * 



4^2 



X 



{xBsf + {M^,^ - t - XB-sf + 2Ml^,^t 



xbs (xbs + t — 



2 

pair 



xT.Qafa/N (Xb) 
a 



(175) 



V.3 



BETHE-HEITLER SUBPROCESS 



The Bethe-Heitler amplitude is calculated from the Feynman diagrams shown in Fig. 
Adopting the notation of the preceding subsection, we write 



iT 



BH 



u {p') i-iQa |e| 7 )u (p) 



^U(A:) i^\e\r) 



{q — k') — m? 



x(^|e|7'^) {q)v{k') 
+u{p')(~iQa\e\^^)u{p)( 
X (i \e\ 7^) V {k') 



u{k) {i \e\Y) (g) 



{k — qf' — rn? 

(176) 



or 



BH 



1^1 Qa — I I 



u {p ) 7^M (p) ey (q) u (k) 



Y 4Y - '^7^k" 
-2 {k' ■ q) 



-J 



-2 {k ■ q) 



v{k') 



(177) 
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Again, after averaging |Tb//| over the polarization states of the initial particles and 
summing over the polarization states of the final parton, we have for the spin-averaged 
square of the amplitude 



BH\ 



2^2 



3 ^2 



BH- 



:i78) 



A straightforward calculation of traces gives 

Cbh 



r/J-P 
^BH 



4 



^BH BbH 
72 ~' 



A 



BH 



B 



BH 



c 



BH 



[{k-qf ■ [k'-qf {k-q){k'-q)_ 
{k ■ q) [k'^'qP + k'Pq^' - g^P {k' ■ q)] - [k^k'P + k^k'^ - g^^ {k ■ k')] 

[k'^qP + k'Pqf" - g^P {k' ■ q)] - {k ■ q) g^P + m^^^", 
{k' ■ q) [k^'qP + kPqP - g^^P {k ■ q)] - [k^k'P + kPk'P - g^'P {k ■ k')] 

[kPqP + kPq" - g^P {k ■ q)] - {k' ■ q) g^P + m^g^P, 
2 {k ■ k') [kPk'P + kPk'^" - g^'P {k ■ k')] - {k ■ q) [k^k'P + kPk'p - 2k'Pk'P] 
- {k' ■ q) [k^k'P + kPk'f" - 2kPkP] - {k ■ k') [k^'qP + kPqP + k'PqP + t-'^-'^l 



+2g^'P{k-k') {k-q) + {k' ■q)-m^ -2m\^'qP 



for the lepton tensor and 



ttBH 
^^^p 



P,j^P'p + PpP'p - 9t.p {P ■ P') 



(179) 



(180) 



for the hadron tensor. 

Recall now that the expression (|178p has to be integrated over the final-state 
momenta constrained by a four-momentum delta function, see Eq. ()157|) . However, 
before contracting tensors, we perform integration of over both lepton three- 
momenta. One defines this particular integral as 

d^k fd^k' 



'-BH 



which, with the help of the formula ()158|1 . turns into 

C'^P^ = Jd^kJ d*k' {q - k- k') 5+ {k^ - m") 5+ (k'^ - m") L 



;i8i) 



MP 
BH- 



;i82) 



We evaluate C^bh center-of-mass frame of the lepton pair. Unfortunately, due 

to presence of a large number of terms in the lepton tensor, the calculation requires 
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more algebra. Let us demonstrate the method by calculating only the integral of the 
first term in the tensor ()179|1 . namely, 



J d^k J d'^k' {q'-k- k') 6+ - 



Ak'f'qP 
JTq) 



(183) 



Integration over k' leads to 

Aq'^qf J d^k 5+ _ ^+ J^^/ _ _ ^ 

-AqP I d^k 5+ - nv') 5+ \{q' - kf - 



(k-q) 
kf 



;i84) 



For simplicity, we choose the incident photon momentum along the z-axis so that 
q = (?0)0,0,go)- The four-momentum of the outgoing lepton is A; = (ko^kj_,kz^. 
Accordingly, the calculation of the first integral in Eq. (jl84j) (we will call it the 
function Jq for later convenience) goes as follows 



Iq = Tc J dk^, I d 
1 

X 



dko 6'^ i /c. 



ki 



go (^0 - k^) 



dky 



2v^ go J-^q'VA-m^ k, - 
^ In 



2 (g ■ q'] 



1 + Jl- 4mVg'2 



1 - Jl -4m2/g' 



/2 



(185) 



Note that we have expressed the initial photon energy in the invariant form, go = 
(g ■ g') / Vq^- The integral over k in the second term of Eq. ()184j) has a Lorentz 
vector-like structure and can be, in general, written as 



\2 2 

m 



hqf + hq'^- (186) 



d^k S+ [k^ - m^) 6+ [{q' - kf 
The scalar functions Ji and I2 are determined by solving the system of two equations, 
^k 6+ - m^) 6+ [{q' - kf - m^] = {q ■ q') h, 
d^kS+ -m^)6+[{q' ~kf -m^ 



(k-q) 



(g ■ q') h + q"h 



;i87) 
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obtained from contracting both sides of Eq. 
some algebra we find 



with and q' respectively. After 
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h 
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7T [/o - 2/2 
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/2 



1 — Am? /q' 
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2a 1 



4^2 



The term p8Hj) now becomes equal to 



Ahq^'qP + 4 (Jo - /a) g'^^g^ 



;i89) 



We repeat this procedure for all the remaining terms in the lepton tensor (all the 
integrals with the corresponding scalar functions are presented in Appendix B), and 
combine them together according to their structure. The final expression for the 
lepton tensor then reads 

27r 



BH 



(g ■ g') 



(190) 



where the coefficients 
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V 
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2 (g ■ q') - 2g'2 + 



2^/2 
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2 (g ■ q') - Aq" + 



(g ■ q') (g ■ g') (g ■ g') 

'4 2m2g'2 
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(g ■ g') (g • g') 



'1 - 



4^2 



12 



(191) 



are scalar functions of the invariants m^, g'^ and (g ■ g') = (g'^ — t) /2. 

For the sake of completeness, it is worth noting that the same method can be 



applied to compute the integral of the lepton tensor, L 



c 
tip 



given by Eq. (jl52j) (note that {k ■ k') + = g'^/2), over the final lepton momenta 
for the Compton subprocess, namely. 



= j d^k j d^k' {q'-k- k') 5+ - m") 5+ {k'^ - m") L^,. 



(192) 



Using 5^'^'> {q' — k — k') to integrate over k' and then writing each term separately we 
get 

+4g; J d^k (5+ (k^ - m^) 5+ [(g' - k)^ - 
-8 j d^k 5+ - m") 5+ [(g' - kf - rr? 
-2gppq'' [ d'k 5+ - m') 5+ [(g' - kf 



kpkp 



m 



(193) 



These integrals can only be constructed out of a four-vector g' and a metric tensor, 
accompanied with the proper scalar functions of g'^. One can immediately discard 
the first two integrals. They are both of the form q'^q'p and, since Eq. p55|) holds, 
they give vanishing contribution, when contracted with Hq'. Moreover, we write the 
third integral as 

ld''k5+{P-m'')5+[iq'-kf-m'']kpkp = N.q'^q'p + N^g^p. (194) 

Contracting both sides with q'^^q'P and g^^ and evaluating the integral, we find the 
solutions for the scalar functions 
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12 
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12 



(195) 
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Again, due to the gauge invariance of H^^ only the term N2g^p survives. Lastly, the 
final integral in Eq. p9Hj) can be read off from Eqs. ()187|1 and p88|l . 



d'^k 6+ - m^) 6+ [{q' - kf - 
By collecting terms, C^p reduces to 



TT 



r,l2 ' 



(196) 



c 
MP 
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(197) 



-8iV2 - 2g'^j g^, 
27r 

~T 

which is identical, apart from the vanishing contribution g^g^, to the expression ()161|) . 

Next we contract Eq. p90|) with the hadron tensor ()180p and follow the steps 
from Eq. (|168|) to Eq. (jl72j) leading to the Bethe-Heitler differential partonic cross 
section, expressed in terms of the subprocess Mandelstam invariants 
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The new coefficients are given by £bh 
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(199) 



Finally, with the help of the master formula (|173p . we obtain the Bethe-Heitler con 
tribution to the inclusive photoproduction of lepton pairs in the parton model 

c^Vbh ir^N 1-1+ X) , 1 
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(200) 
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V.4 



INTERFERENCE TERMS 



We turn our attention to the interference between the Compton and Bethe-Heitler 
subprocesses. Four Feynman diagrams give rise to eight interference terms. They all 
are calculated in a similar way. For illustration, let us consider the first two terms. 
Up to a common factor we have 
1 



|Ti| 



{k' ■ q) 
1 

Wo) 



Tr (^+ i) 7" ^75] Tr ( ^ + m) 7^ ( ^' - m) 7, ( ^- ^' + m) 7^ 



Tr [/7^ ^) 7*" ^75] Tr + m) 'y^{fi' ~ m) 'y^ ^ + m) 7^, 



(201) 



After evaluating both second traces in the above expressions, we observe that they 
differ, apart from the sign, by the interchange k ^ k' . Thus integrating the sum of 
|Ti|^ and |T2|^ over the final lepton four-momenta yields zero. 



j d^k J d^k' (|Ti|' + IT2I') = 0. 



(202) 



By the same token, the remaining six interference terms mutually cancel each other, 
after being integrated over k and k' . 

In fact, there is a more elegant way (without cumbersome and explicit calculation 
of traces) to see the cancellation of the interference terms. Note that one needs to 
compare the following two double integrals 

1 



d^k j d^k' {q'-k- k') ■ 
d^k I d^k' {q'-k- k' 



{k' ■ q) 

. 1 



{k-q) 



Tr 

Tr 



(/(f + m) 7^ {}i' - m) 7^ {4- ^' + m) 7' 
/{f + m)7^(/^'-m)7^(/{f- ^ + m) 7^ 



(203) 



Due to the presence of the delta function, one of the integrals, e.g. the integral over 
fc', can be trivially carried out. In addition, using the fact that the trace does not 
change under the cyclic permutation, Eq. ()203|) can be rewritten into 
1 



d^k 



[{q' -k)-q] 



d^k 



Tr 



1 



+ m) 



{k-q) 



Tr 



' + m) 'y^ {^ + m) 



(204) 



Transposing now the trace in the first integrand, and then performing the momentum 
shift, k = q' — k, this integral turns into 



d^k 



1 



(k-q) 



Tr 



(/-/(f + m) 7^^(^-/if + m) 7j(^/{f-m) 7^ 



(205) 
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Using the property C^^C^^ = — 7J, where C = 27^7" represents the charge conjuga- 
tion operator, and further permuting the last two factors 7^C~^ under the trace to 
the left side, we convert the integral finally into the form 

1 



7^ (/- ^ - m) 7« (^- ^ + m) 7, + m)] . (206) 



- / rfU -j^ ^Tr 

J (k-q) 

The latter is, apart from the sign, identical to the second integral of Eq. ()204|) and 
therefore, the two integrals exactly cancel each other. This result of cancellation is, in 
general, known as the Furry's theorem. It states that Feynman diagrams containing 
a closed fermion loop with an odd number of photon vertices can be omitted in the 
calculation of physical processes |103j . It should be pointed out, however, that the 
cancellation of the interference terms happens only after integration over the momenta 
of the final leptons. These terms can still be accessed if one measures, for example, 
the angular distribution of the lepton pair. 

V.5 KINEMATICS 

As the Bjorken scaling variable xb tends to zero, both the Compton and Bethe-Heitler 
differential cross sections, see Eqs. fll75|) and ()200p . become singular. Nevertheless, 
the following kinematical constraint 

s-M^,^. >{q~q' + Pf> (207) 

holds, where M stands for the nucleon mass. Thus, in the laboratory frame, where 
the real and virtual photon momenta are k = (00, k^ and k' = (uo', k'^, respectively, 
and the nucleon momentum is P = (m, 0), the variable xb satisfies 

Xb > 7 ^ — . 208 

1 - {2Mu - Ml„) It 

For the invariant momentum transfer t we have 



t = Ml^„ - 2uj (uj' - _ M2^.^ cos^^^) , (209) 

where is the scattering angle, i.e. the angle between the directions of the initial 
and final photons. 

We consider a proton target and the muons as the outgoing pair of leptons and 
thus we have M ~ 0.94 GeV and m ^ 0.106 GeV. In Fig. |2S1 the Bjorken scaling 
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FIG. 23: Bjorken scaling variable xb plotted as a function of the angle 6^^ between the 
incoming real and outgoing virtual photon in the target rest frame for uj' = 10 GeV 



and M, 



pair 



3 GeV with uj = AO GeV photon beam. 



variable given by Eq. p74|) . is plotted against the angle 6^^ for fixed values of the 
incoming photon beam energy, u = 40 GeV, the energy of the outgoing virtual 
photon, u' = 10 GeV, and the lepton-pair mass, Mpair = 3 GeV. We find 0.49 < 
Xb < 0.96 for the scattering angles < 6^/.y < 15°. In this region of xb, the valence 
quarks play the dominant role, in particular the ^-flavor contribution, as shown in Fig. 
I2H1 Next we plot the angular dependence of the invariant t for the same kinematics, 
see Fig. |^ Due to high energy of the photon beam, rather large values of t are 
expected, namely, 27.85 GeV^ < — ^ < 53.85 GeV^. It is important to note that 
our description is not limited only to the small invariant momentum transfer, like in 
the case of DVCS kinematics discussed in Chapter IV. In contrary, it works for any 
kinematically allowed value of t. 

To estimate the cross sections, we use a simplified parametrization of the unpo- 
larized valence quark PDFs in the proton, given by Eq. ()124|) . together with the sea 
quark distribution, 

sea{x) = 0.5x"°-^^ (1 - x)^ (210) 

taken from Ref. These are, however, the valence distributions for the u- and 

(i-flavor components in the proton, respectively, and the sea distribution in the proton 
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FIG. 24: Invariant momentum transfer t plotted as a function of the angle Q^^ between 
the incoming real and outgoing virtual photon in the target rest frame for lo' — 10 GeV 
and M 



pair 



3 GeV with a; = 40 GeV photon beam. 




FIG. 25: Unpolarized valence quark distributions for u-fiavor (bold solid line) and 
d-flavor (solid line), and the sea quark distribution (dashed line) in the proton. 
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FIG. 26: Compton cross section ac (solid line) and Betlie-Heitler cross section asH 
(dashed line) plotted as a function of the angle 9^^ between the incoming real and 
outgoing virtual photon in the target rest frame for uj' 
with uj = GeV photon beam. 



10 GeV and M 



pair 



3 GeV 



whereas the distributions, which appear in the cross section formulas ()175|) and ()20m) 
are the quark and antiquark proton PDFs. In fact, what one really needs is the sum 
of the quark and antiquark distributions, 

E^/.(^) = E [/p (^) + /p (^)] • (211) 
/ 

As noted in Section IV. 2, sum ^2\\\ corresponds to adding the valence PDFs and 
twice the contribution from sea quarks. 

The unpolarized cross section results are presented in Fig. 1211 We observe that 
the Compton contribution dominates by a factor ~ 3 in the forward direction, i.e. 
for 9^^ < 3*^, and falls off rapidly, even below the Bethe Heitler contribution for the 
scattering angles 9^^ > 10°. 

In summary, we have demonstrated the application of the QCD parton model to 
the high-energy photoproduction of lepton pairs. The Compton and Bethe-Heitler 
cross sections were estimated for a given simple parametrization of usual parton 
distributions in the proton. Without explicit calculation of the interference terms, 
we have shown their mutual cancellation, after being integrated over the final lepton 
momenta. 
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CHAPTER VI 

EXCLUSIVE PHOTOPRODUCTION OF LEPTON PAIRS 

VI.l INTRODUCTION 

Deeply virtual Compton scattering is considered to be theoretically the cleanest and 
simplest process, in which GPDs can be accessed. We can now slightly modify the 
process, simply by interchanging the two photons at the subprocess level, 

7(gi) + iV(pi) 7*fe) + iV(p2). (212) 

Then one arrives at time-like Compton scattering (TCS) ^104J which is, in principle, 
the inverse process to DVCS. In a specific kinematical regime, where the virtuality of 
the final-state photon is large, g| oo, TCS becomes a handbag dominated process 
and accordingly, it can be studied in same way as DVCS was studied in Chapter IV. 

Time-like Compton scattering can be accessed through the physical process known 
as the exclusive photoproduction of lepton pairs. At this point, we need to draw a 
distinction between the exclusive and the inclusive photoproduction of lepton pairs. In 
the former case we detect, in addition to the lepton pair, also the scattered nucleon N 
(with the four-momentum as illustrated in Fig. EZI Another important difference 
is that unlike the inclusive photoproduction, we are limited now to small invariant 
momentum transfer t = {pi — ^2)^ to the nucleon. Recall that the DVCS kinematics 
requires small t compared to the invariant g| (the DVCS scale Qf is now replaced 
by the virtuality and thus the terms O {—t/q^) in the amplitude can be safely 
dropped out. On the other hand, the common feature of both inclusive and exclusive 
photoproduction processes is the presence of two types of diagrams. In particular, for 
the exclusive photoproduction of lepton pairs, they are shown in Fig. EHl 

In Section VI. 2, we define the kinematics. In the next section, we compute the 
leading- twist amplitude for the Compton part. It can be easily derived from the 
result for VCA established in Section IV. 2. The relevant integrals in the Compton 
and Bethe-Heitler cross sections can be evaluated simply by following the steps of 
Sections V.2 and V.3. Finally, the cross sections are estimated by using the same 
model for the nucleon GPDs as in the DVCS process. 
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FIG. 27: Exclusive photoproduction of lepton pairs. 




a b 



FIG. 28: DVCS (a) and Bethe-Heitler (b and c) diagrams contributing to the exclusive 
photoproduction of a lepton pair. 



VI.2 KINEMATICS 

In the laboratory reference frame, as the target rest frame, the photon and nucleon 

four-momenta arc denoted by qi = (z/i,g!), q2 = (^'2,^2), Pi = (-^,0^ and p2 = 
(£"2,^2) where M stands, as usual, for the nucleon mass. One rotates the coordinate 
system so that the 2;-axis is in the direction of the incident photon in the photon 
scattering plane, formed by the photon thrcc-momcnta qi and q2. Similarly, the final 
leptons define the lepton scattering plane. Wc write their four-momenta as k = (^u, 
and k' = (uj',k'^, respectively. In this reference frame, the invariant momentum 
transfer can be expressed as a function of the angle 9^^ between the incoming real 
and outgoing virtual photon, namely, 

t ={q2-qif= ql - 2z/i [1/2 - \] i^l - ll cosO^^ . (213) 

On the other hand, the invariant t is also equal to 



t = 2M(i/2-i^i), 



(214) 
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which follows from the expression t = 2M^ — 2ME2 and the energy conservation, 
M + Ui = E2 + i'2- Combining Eqs. ()21H|1 and ()214|1 gives the quadratic equation in 
the energy of the virtual photon z/2 with the solutions 



The coefficients are 



AB ± JA^ - ql (B"^ - 1) 



A ^ ^1±4^>M, 

cos 9yy 

B = ' > 1, (216) 



cos 9^^ 



where the scaling variable, 



Y = ^ = , (217) 

was introduced. It is an analog of the Bjorken scaling variable, xb = Qif^ {pi ■ Qi), in 
the DVCS process. Since in the DVCS kinematics the invariant momentum transfer t 
is kept as small as possible, only the plus solution of Eq. ()215|) should be considered. 

To illustrate the kinematics, we set the beam energy of initial photons to ui = 
5 GeV, and the invariant mass of the lepton pair to = 3 GeV^. Accordingly, the 
scaling variable (j217|) is fixed at x = 0.32. The invariant momentum is plotted as a 
function of the scattering angle in Fig. |^ It covers the range 0.15 GeV^ < — ^ < 
1.19 GeV^ for 6*^^ up to 10°. For the same region in 9^^, we also plot the energy E2 
of the final-state nucleon, see FigEHl and obtain 1.019 GeV < E2 < 1.574 GeV. 

VI.3 CROSS SECTION 

In general, the differential cross section for the photoproduction of lepton pairs is 

where the invariant s = {pi + gi)^ = 2Mz/i + M^. Similarly to the DVCS case, the 
invariant matrix element consists of two contributions. Let us first investigate the 
Compton part. 



78 



t [Ge\2] 

-A Oyy [d^] 




FIG. 29: Invariant momentum transfer t plotted as a function of the angle 6^,^ between 
the incoming real and outgoing virtual photon in the target rest frame for the invariant 
mass of the lepton pair g| = 3 GeV^ and x = 0.32 with ui — 5 GeV photon beam. 
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FIG. 30: Energy of the scattered nucleon E2 plotted as a function of the angle 6^^ 
between the incoming real and outgoing virtual photon in the target rest frame for 
the invariant mass of the lepton pair = 3 GeV^ and x = 0.32 with z/i = 5 GeV 
photon beam. 
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VI. 3.1 Compton Process 



The T-matrix of the Compton contribution is given by the diagram (a) in Fig. EHl 
Using Ey (gi) to denote the polarization of the initial photon, we have 

(219) 



Q2 



where T^^cs reduced VGA for time-like Compton scattering, namely. 



^TCS 



I / d'z e^(«-) (iV s,) I T { J^A, {z/2) r^,, i~z/2)} \N {p„ s,)) . 



(220) 



In the twist-2 approximation, this amplitude is calculated from two {s- and w-channel) 
handbag diagrams illustrated in Fig. El Unlike the standard electromagnetic DVCS 
process, the incoming photon is now real, qf = 0, and the outgoing photon is time- 
like, g| > 0. The final result for T^j^^. at the twist-2 level in terms of OFPDs can be 
derived from the expression (|lU2j) for the reduced VGA in the DVCS process simply 
by replacing the final photon momentum with the initial one, q2 —^Qi, and changing 
the sign in the skewness parameter rj. Recall from Section IV.2 that the latter is 
defined as = [r ■ q) /2 {p ■ q), where r = pi — p2 is the overall momentum transfer 
to the nucleon, and p = {pi + P2) /2 and q = (gi + ^2) /2 are the average nucleon 
and photon momenta, respectively. As for DVCS rj = xb/ {2 — xb) in the DVCS 
kinematics, we have for TCS a similar relation, i.e. r] = (2 — x). Moreover, the 
scaling variable ^ = — g^/2 {p ■ q) now coincides with —77. 

Substituting x — x and using the symmetry properties of OFPDs given by Eq. 
(P7j). we obtain the reduced VCA for TCS at the leading- twist level in the form 

1 ■s-^ „o /"^ dx 



^TCStwist-2 



Y.Q} 



-1 [X 



1] — iO) 



x{Hf{x,r],t) 



+El {x,r],t) 



{P ■ Qi) 

{p'qi+p''q'i)-g''' 
1 



U{P2,S2) 4lu{pi,Si) 



{P-Qi 

+Hfix,ri,t) 
1 



(P ■ (I2] 



W (P2, S2) — U [pi, Si) 

uip2,S2) 4ll-,u{pi,Si) 



^e^'"^qlpP^ 



Et {X,7],t) 



(p-qi 



^^^M (P2,S2)75m(Pi,Si) 



(221) 



Since the photons in the amplitude ()22H) are interchanged compared to the standard 
VCA, the electromagnetic gauge invariance of '^TCstwist-2 is now satisfied with respect 
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to the initial photon, but violated in O (r^) with respect to the final one. Similarly 
to the standard DVCS process, we introduce the following integrals 

dx 
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[x 


-V- 


iO) 




dx 




[x 
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iO) 




dx 




[x 


-V- 
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iO) 



H 
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x,r],t) , 
x,r],t) , 
x,r],t), 
E+{x,r],t). 



H 



H 



(222) 



They can be calculated with the help of Eq. ()103|) . For illustration, we plot 
in Figs. 1^ and 221 the real and imaginary parts of the convolution integral 
'H^(jg{rj,t) = 'H^^jgj^^-j {ri,t) + 'H^^jgj^^-^ {ri,t). We present the u- and c?-quark con- 
tributions separately, using a GPD model from Section IV. 4. To get rid of the 
^-dependence, we divide by appropriate factors, Fiuit) /2 and Fid it). Analogous 
plots for Ti-Tcs (^)^) = '^tcsm iv^t) + '^TCs(d) iv^t) (here we divide by the factor 



9a (t) /gA(t = 0) to remove the dependence on t), are given in Figs. ESI and EH All 
four figures clearly show that, for the proton target, the valence u-quark contribution 
to the real and imaginary parts of the integrals 'H'^cs iv^ and 'H^cs ili is much 
larger compared to the valence c?-quark contribution. 

We will restrict ourselves to the unpolarized incoming photon beam and the un- 
polarized nucleon target. Then the spin-averaged square of the T- matrix in the fac- 
torized form is 



-c 



2g| 



r C TjlJ-a 



(223) 



where the lepton and hadron tensors are 



'C 
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(224) 
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FIG. 31: Contributions from u quarks (solid line) and d quarks (dashed line) to the 
real part of H^cs- They are divided by Fi„ (t) /2 and Fi^; (t), respectively. 




FIG. 32: Contributions from u quarks (solid line) and d quarks (dashed line) to the 
imaginary part of H-tcs- They are divided by Fi„ {t) /2 and Fi^ (t), respectively. 
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FIG. 34: Contributions from u quarks (solid line) and d quarks (dashed line) to the 
imaginary part of Ti-Tcs divided qa if) /gA{t = 0). 
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with m denoting the lepton mass. Note that if^" has the same structure as its DVCS 
partner given by the expression p38|l . To calculate the cross section, we integrate 
Eq. ()223|1 over the Lorentz-invariant phase space. This is done exactly in the same 
way as in Section V.2 for the case of the inclusive photoproduction of lepton pairs. 
Thus first we integrate Lj^^ over the part of the phase space associated with both 
final leptons, 



2^ 



2uj' 



7 5''^ iPi 



qi -P2 



X 



k - k') L% = I d% (pi + q,-p2 
27rg| + 2m2 / 



52 J 



q29fia 



q2/iq2c 



ql 



Am? 
If, 



(225) 



and then contract the tensor [q2g^j.a — q2ij.q2a] with i?^". It is easy to see that, after 
contraction, the first term is of the order g| while the second term is (9(M^). The 
latter can be therefore ignored in the limit M^/g| — > 0. We are still left with two 
more integrations. One of them, e.g. integration over p2, is trivial due to the presence 
of the delta function S^^^ (pi + qi — P2 — q2) on the right-hand side of Eq. (|225|) . The 
cross section now turns into 

'l-^ fd%S+ ip, + q,-Q,?-M'^ 



37T 



sq2 



ql 



(226) 



where the contraction in the square brackets reads 
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The remaining integral in Eq. ()226|) is performed in the photon-nucleon center-of- 
mass-frame, in which the nucleon mass can be safely neglected. The final result for 
the unpolarized Compton differential cross section in the invariant form is 

+ Im^ 



d^ac 
dM^^i^dt 



a 



pair 



6 



A 



1 - 



■'"pair 



(228) 



VI.3.2 



Bethe-Heitler Process 



The Bethe-Heitler amplitude. 



BH 



-ey {qx)u{k) 
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x{Nip2,S2)\J,{0)\Nip,,si)) 



(229) 



comes from the diagrams (b) and (c) in Fig. EHl Dirac and Pauli nucleon form factors 
parametrize the nucleon transition current matrix element, 



{N{p2,S2)\J^,{0)\N{p^,S,)) = U{P2,S2) 



2M 



(230) 



After averaging and summing |Tb/^| over the polarizations of the initial and final 
particles, respectively, one has 

^^^^ (4vra)' 



BH\ 



J IjBH 



(231) 



where both tensors have already been calculated. The hadron tensor is obtained by 
replacing the Lorentz indices z/ — /i and (3 ^ a m the expression 
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and the lepton tensor by replacing the momenta q qi and indices p 
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(232) 
a in Eq. 
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The integral of L'^^j over the final lepton momenta is known from Section V.3, 
namely, 



BH 



fj.a 
BH 
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The unpolarized Bethe-Heitler cross section is then 
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and finally in the differential invariant form, 

d?aBH CK^ 1 



At: sH'^ 



pixa ttBH 
'-'BH^^Ixa 



(236) 



(237) 



86 



The contraction in the square brackets of Eq, 
ttBH _ 47r I Arm I . t 




t M2 

1 + + ^ 

s s 2s 











pair / 



X 




M2 , 
+AVbh— ( 1 + 



'f^ (t) - -^Fi (t) 



(238) 



Note that and (gi ■ ^2), which appear in the coefficients ()235j) . should be replaced 
by M^^i^, and [M^air - t) /2, respectively. 

In Figs. EHl and EHl we plot the Compton and Bethe-Heitler contributions using 
a toy model for the proton valence OFPDs described in Section IV.4. Both cross 
sections are shown together on a logarithmic scale in Fig. |23 We notice that, 
similarly to the standard DVCS, the Bethe-Heitler process is again the dominating 
contribution (i.e. ^ 0.95) to the cross section. 

We conclude with the remark that the cross sections presented here were inte- 
grated over the angles of the final-state leptons. One may expect to get additional 
information of GPDs, by measuring the cross sections that are sensitive to the an- 
gular distribution of final leptons. In particular, the interference terms between the 
Compton and Bethe-Heitler processes, together with the use of the polarized photon 
beam, may provide a new insight. This is an interesting subject by itself, and requires 
more attention in the future. 
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FIG. 35: Compton cross section oc plotted as a function of the angle 6*^^ between the 
incoming real and outgoing virtual photon in the target rest frame for M^^^j, = 3 GeV^ 
and X = 0.32 with i/i = 5 GeV photon beam. 




FIG. 36: Bethe-Heitler cross section gbh plotted plotted as a function of the angle 
^-y-y between the incoming real and outgoing virtual photon in the target rest frame 
for = 3 GeV^ and x = 0.32 with i/i = 5 GeV photon beam. 
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FIG. 37: Compton cross section ac (solid line) and Bethe-Heitler cross section asH 
(dashed line) plotted as a function of the angle 9^^ between the incoming real and 
outgoing virtual photon in the target rest frame for M^^^^^ — 3 GeV^ and x = 0-32 
with 1/1=5 GeV photon beam. 
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CHAPTER VII 
WEAK DEEPLY VIRTUAL COMPTON SCATTERING 
VII. 1 INTRODUCTION 

Complementary to the standard electromagnetic DVCS process, different combina- 
tions of quark flavors can be accessed by utilizing the weak current, which couples 
to quarks with strengths proportional to the quark weak charges. In analogy with 
DVCS, these can be studied in neutrino- induced virtual Compton scattering, 

p{k) + N{p{) zy(A;')+V(p2)+7fe) (239) 

for the neutral current, and 

v{k) + N{p,) f,^(k') + N'{p2)+i{q2) (240) 

for the charged current reactions (and similarly for antineutrinos). Note that in the 
case of the latter the initial and final nucleons will be different. Because of the V-A 
nature of the weak interactions, one can probe C-odd combinations of GPDs as well as 
C-even, where C is the charge conjugation operator, and thus measure independently 
both the valence and sea content of GPDs. The weak current also allows one to 
study flavor nondiagonal GPDs, such as those associated with the neutron-to-proton 
transitions in charged current reactions in Eq. ()240|) . The use of weak currents can 
thus provide an important tool to complement the study of GPDs in more familiar 
electron-induced DVCS or exclusive meson production processes. 

In this chapter, we study the weak deeply virtual Compton scattering processes 
given by Eqs. ()240|) and ()239|) . and present a comprehensive account of the ampli- 
tudes and cross sections in the kinematics relevant to future high-intensity neutrino 
experiments |lU5j . In Section VII. 2, we derive both the weak neutral and weak 
charged VCAs using the nonlocal light-cone OPE, and introduce an appropriate set 
of off-forward parton distributions, which parametrize the weak DVCS reactions. The 
weak DVCS processes are analyzed in Section VII. 3. Using a simple model for nu- 
cleon OFPDs from Section IV.4 (recall that the model does not include the sea quark 
contribution), we estimate the cross sections and compare the respective rates in neu- 
trino scattering with those in the standard DVCS process. In addition, for the sake 



90 



of completeness, we also discuss the electron-induced DVCS process associated with 
the exchange of the weak boson Z^. 

Deeply virtual neutrino scattering has been discussed recently in Ref. |10fij for 
neutral currents and in Ref. |107j for charged currents. Also a preliminary report 
containing some of the formal results from Ref. |109j appeared in Ref. |108j . 

VII.2 WEAK VIRTUAL COMPTON SCATTERING AMPLITUDE 

In this section, we discuss the amplitudes for the weak virtual Compton scattering 
process. Before turning to the specific amplitudes for the neutral and charged current 
cases, we review some general aspects that have been already discussed in Chapter 
IV. 

At the subprocess level, in analogy with the photon-induced DVCS amplitude, 
the weak virtual Compton scattering amplitude can be obtained by simply replacing 
the incoming virtual photon with the weak boson B, 

B{q,) + N{p^) 7(g2) + Ar'(p2). (241) 

where B is either or W^. Note that in the electromagnetic and weak neutral 
cases both the incoming and outgoing nucleons are the same, N = N'. Similarly to 
the electromagnetic VCA introduced in Section IV.2, for the weak process with an 
incoming or Z^ boson and outgoing photon, one has 

= i j J d^y e-*(''i-)+*{g2-y) 

X {N' (p2, 52)1 T { J^,, (y) (x)} \N (pi, s,)) , (242) 

where J^f corresponds either to the weak neutral current or the weak charged current. 
We will denote the currents by and J^q^ respectively. Again, due to the current 
conservation 

T^V = and q2^.T^ = Q. (243) 

In terms of symmetric coordinates, X = {x + y) /2 and z = y — x, and symmetric 
momentum variables, q = {qi + ^2) /2 and p = {pi + P2) /2, the weak VCA takes the 
following form 

T/^r = (2vr)^5(^)(pl + gl-p2-g2)r4^^ (244) 
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FIG. 38: Handbag diagrams (s- and w-channel) in the weak virtual Compton scat- 
tering amplitude. 



with the expression for the reduced weak VGA 

W = rf';^e^(^-^)(iV'(p-r/2,S2)|T{J^^,(z/2)J^(-z/2)}|iV(p + r/2,si)). 

(245) 

Recall that, in general, the initial and final nucleons may be different, p\ = Mf and 
pI = M|. Accordingly, Eq. (jHTj) is now replaced by 

p2 = l(M2 + M2-t/2) and {p ■ r) = ^ [m', - M^) , (246) 

however, it is recovered by neglecting the mass difference between the proton and 
neutron (we set Mi = M2 = M). 

Our starting point is the coordinate representation of the time-ordered product of 
the weak and electromagnetic currents. The leading light-cone singularity is contained 
in the handbag contribution illustrated in Fig. EHl Note that the weak current couples 
to the quark fields through two types of vertices, qqZ^ and qqW^. Therefore, the 
quark fields at coordinates ±z / 2 can carry either the same or different flavor quantum 
numbers. We study these two cases separately. 



VII. 2.1 Weak Neutral Amplitude 

We start with the expansion of the time-ordered product of the weak neutral and 
electromagnetic current. Apart from an overall factor — \e\g/ cos^iy (recall that g is 
the coupling constant of the weak interaction), we have 

^T{r^^{z/2)J^A-zm} = I E Qf (^/2) ^ {z) 7- (^ '^^ ~ 



/ 

x^^(-z/2)+^;(-z/2)7'^ 



Cy - 75CA 



^(-^)7^^;(^/2) , (247) 
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where Cy and c;^ are the weak vector and axial vector charges, respectively. They are 
given in the standard model by 



u,c,t 
d,s,b 



u,c,t 
d,s,b 



l/2-2g„,e,tsin2^ 



-1/2 - 2Qd,s,bSm^ Ow, 

1/2, 

-1/2, 



(248) 



where 9w is the Weinberg angle with sin^ 6]y ~ 0.23. After substituting the explicit 
expression for the free quark propagator in the coordinate representation, see Eq. 
()88|) . and using the formula ()89|) . the expansion ()247|) can be written in a compact 
form in terms of the string operators 



tT{r^,,iz/2)J^^i-z/2)} 



Att'^z' 



f 



{4[. 



i,f{-z/2)^,i,f {z/2)-{z^-z) 



+ie 

Ca 
+ie 



fipur] 



tPf {-z/2) -f^-f^iJf {z/2) + {z^ -z) 



ppuri 



(-^/2) 7ry75^/ {.z/2) - ^ -z) 

^f{-z/2)^^,ljf{z/2) + {z-^-z)\]]. 

(249) 



In contrast to the standard electromagnetic DVCS process, see Eq. fl90|l . we end 
up with two additional terms. Namely, the presence of the axial part 750;^ of the 
V-A interaction gives rise to a vector current symmetric in the Lorentz indices /i, v 
and to an axial vector current antisymmetric in /i, v. The string operators here are 
accompanied either with QfC^ or Qfc{ rather than Q). Also the denominator of Eq. 
carries an extra factor of 2 due to the nature of the vertex qqZ^ . 
Accordingly, to obtain the twist-2 part of Eq. ()249p we need, in addition to 



contracted string operators given by Eq. ()94|) . two extra operators. One can conve- 
niently write all four contracted vector and axial vector string operators under the 
same footing as 



Of^{z) = 
Ot(z) = 



i,f {-z/2) {z/2) ±{z^-z 
ijf{-z/2) i^^i,f{z/2)±{z^-z) 



(250) 



They all should satisfy the harmonic condition (jHSl)- Furthermore, the parametriza- 
tion of the nonforward matrix elements of these operators is performed by means of 
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the relevant nonperturbative functions, namely, the off-forward parton distributions. 
A straightforward generalization of the parametrization ()9(j|l from Section IV. 2 can 
be performed by introducing a new set of GPDs, i.e. the minus OFPDs. Then 

{N {p2, S2)\ {z) \N (pi, si))^2=o = w (P2, S2) iu (pi. Si) 



X 



y dx e'^'^P-'^Hf {x, e, t) 



X 



j'^dxe'''^P-'^Ef{x,^,t) 



{N {p2,S2)\Ol^ {z)\N {pi,Si))^2^Q = u{p2,S2) ^-f5u{puSi) 

X rfxe"(^")i7j(x,e,t) 



(r • z) 

-U (P2,S2)-^^75m(Pi,Si 



X 



l'^dxe'''^P-'^Ef{x,^,t). (251) 



Note that (z) has a superscript opposite in sign with respect to the corresponding 
tilde OFPD. While the standard DVCS process gives access only to the plus distribu- 
tions (recall that they correspond to the sum of quark and antiquark distributions), 
scattering via the weak virtual boson exchange probes also the minus distributions. 
It can be shown that the latter correspond to the difference in quark and antiquark 
distributions, i.e. the valence configuration. Like the plus OFPDs, the minus distri- 
butions have also symmetry with respect to the change x —>■ —x. Let us summarize 
the symmetry properties in x for both plus and minus OFPDs, 

Hf{x) = THf{-x), 

Ef{x) = TEf{-x), 

Hf{x) = ±Hf{-x), 

Ef{x) = ±Ef{~x). (252) 

After following the steps described in Section IV.2, we arrive at the expression for 
the reduced weak neutral VGA in the leading-twist approximation 

dx 



1 /•! 

'Twm^^si-'i = E Qf J ^ _ ^ + ,0) 



x{4h+ (x,e,t) 



U{P2,S2) 42u{pi,Si 
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1 



{p ■ (12) 



1 



-ie'''''q2pPr, 



^ VP2, S2) ^^^^^ U {pi, Si) 

U{P2:S2) !i2l5U {Pl, Si) 

{<12 ■ r) 



te'""^q2pP, 



2M 



-U {P2: S2) {Pu Si) 



-c{Hj{x,tt) 



+ciEj{x,^,t) 



1 



{p'q2+P''q^)-9"' 



ip-q2 

1 



U{P2,S2) i2lbU{pi,Si) 



-ciEf{x,^,t) 



(p ■ 52) 

+c{Hj{x,^,t) 

1 



(g2 • r) 
2M 



u{p2,S2)j5uipi,Si) 



{P • Q2) 



te'"'''^q2pP, 



m(P2,S2) £uipi,Si) 



{P ■ (I2) 



te'"""^q2pPr, 



« (P2, S2) — — — u [pi, Si) y 



(253) 



VII.2.2 Weak Charged Amplitude 

Skipping an overall factor — |e| g/V^, the expansion of the time-ordered product of 
two currents in the weak charged sector reads 



iT{Jj^^{z/2)J^ai-^/m 



E {Qf^ijf, {z/2) YYl'' (1 - 75) V'/ {-z/2) 



fj' 



-Qf^f {-z/2) Y (1 - o'5) YYy^f (^/2)} 



fJ' 

-if^f\z^-z)] 
+ie^P^^[Qfi,f{-z/2)^,^,,Pf, {z/2) 
+ {f^f,z^-z)] 
-s^'^'^^[Qfi^f{-z/2)^,^,^f. {z/2) 
-{f^f',z^-z)] 
-ie^P'^r,^Qf^f{-z/2)j,^j, {z/2) 
+ {f^f,z^-z)]}. 

In principle, one has to include the CKM matrix since for the weak charged couplings 
the quark mass eigenstates do not coincide with the weak eigenstates. However, we 
shall ignore all Cabibbo angles. Here the sum over quark flavors is subject to an 
extra condition, — Qf, — 1 or —1, due the fact that the weak virtual boson 



(254) 
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carries an electric charge ±1 in units of |e|. For that reason, the initial and final 
nucleons are not the same particles anymore. Hence we are dealing either with the 
neutron-to-proton transition via the exchange of W^, or with the proton-to-neutron 
transition via the exchange of W~ . One notices that the vector and axial vector 
string operators in the expansion ()254|) are not diagonal in quark fiavor, i.e. they are 
accompanied by different quark fiavors as well as by different electric charges. The 
corresponding contracted string operators, which appear when extracting the twist-2 
part of Eq. ()254|) . can be expressed as the linear combinations. 



Q;^; {-z/2) {z/2) ±(f^f\z^ -z) 
Qf^f {-z/2) /75^^. {z/2) ±{f^f\z^ -z) 



Q^Off'+ {z) + Q^Off- {z) , 

(255) 



Q±0(^'+ {z) + Q^Qf/- {z) , 



of operators. 



O^f'^ {z) 
Oi^'^ {z) 



i>f {-z/2) {z/2) ±{f^f\z^ -z) 
'i,f {-z/2) /75^; {z/2) ±{f^f,z^ -z 



(256) 



with the coefficients Q± = {Q f ± Q f/) /2. The matrix elements of these newly intro- 
duced operators are parametrized in terms of the fiavor nondiagonal OFPDs, 



{N'{p2,S2)\0^f'^{z) \N{p„s,)) 



{N' {p2,S2)\Ol^'^ {z)\N{p^,s^)) 



z2=0 



U{P2,S2) ^u{pi,Si) 

X J dxe'''^P-'^Hff,{x,^,t) 

(^2, S2) — U {pi, Si) 

X J^^dx e'^^P-'^Eff, (x,e,t), 

U{P2,S2) /75u(pi,Si) 

X J'^dx e'^'^P-'^ Hff, {x,^,t) 
(t ' z) 

~U {P2,S2) -17V^l5U {pi,Si) 



2M 



X 



J^^ dx e'^^P-'^Ejj, {x, ^, t) . (257) 



Finally, with the help of Eqs. ()255|) and ()257|) . the leading- twist reduced weak charged 
VGA can be easily obtained from the result for the weak neutral case with the proper 
replacements. We get 



dx 
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X • 



Q+H+, {x, ^, t) + Q.Hjf, (x, t) 



X 



X 



+ 



{P ■ (12, 



{p'q2+P''<l^)-r 



Q+Ep {x,^,t) + Q^Ejf, (a;,e,t) 

^ (^2, S2) — U (pi, Si) 



Q+Hp {x, t) + Q.Hjf, {x, e, t) 



X 



U{P2,S2) ^2l5uipi,Si) 

\P ■ (I2) J 

+ [Q+Ejf> (x, e, t) + g_^7^, (x, e, t) 



X 



^e^''''y2pPn 



-^^m(P2,S2)75M(Pi,Si, 



X 



(P • <l2) 



{P'q2 + ^^^92^) - 9 



X 



(P • Q2) 



{p'q2+P''q^)-9''' 



U{P2,S2) i2lbU{pi,Si) 

+ [Q-Ejf, (x, e, + (x, i)" 

<l2-r)_ 

-m(P2,S2)75«(Pi,Si) 



2M 



X 



' m(P2,S2) ^2m(Pi,Si) 



(P • 52) 



X 



(P • 52) 



+ 



U (P2, S2) — U (pi, Si) 



AM 



(258) 



VII.3 WEAK DVCS PROCESSES 



In the following section, we discuss three examples of the weak DVCS processes. Their 
kinematics has ahcady been analyzed in detail in Section IV.3. The only difference 
is that the incoming virtual photon 7* of the standard DVCS process is now being 
substituted by the virtual weak boson B. We use the OFPD model from Section 
IV.4. Note that by neglecting the sea quark contribution, the plus OFPDs become 
equal to the minus ones, and hence we write Hj^ = Hj = HJ"'^ H~j — Hj = 
E^ — Ej = EY^ and E~^ — Ej = EJ'^K For the weak neutral current scattering 
process, we discuss neutrino-nucleon and electron-nucleon scattering. We present 
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FIG. 39: DVCS (a) and Bethe-Heitler (b and c) diagrams in the weak deeply vir- 
tual Compton scattering process. The incoming leptons can be either neutrinos or 
electrons. 



the calculation by taking an unpolarized proton as a target nucleon. In the case 
of an unpolarized neutron target, one can simply use the isospin (i.e. the charge) 
symmetry to express the neutron GPDs in terms of the proton ones. In the weak 
charged current interaction sector, we consider neutrino scattering off a neutron via 
the exchange of a boson producing a proton in the final state. All three examples 
are illustrated in one figure, see Fig. EHl As in Section IV.5, the relevant cross sections 
are plotted against the scattering angle Os-y between the directions of the incoming 
weak virtual boson and the outgoing real photon for one kinematical point with 
Qi = 2.5 GeV^, xb = 0.35 and = 0, however, for two different lepton beam energies. 
We take u = 5 GeV and 20 GeV rather than uj = 5.75 GeV and 11 GeV. Similarly 
to the standard DVCS process, the Bethe-Heitler cross sections are presented on a 
logarithmic scale as well as the plots with both contributions. In addition, for given 
kinematical regions, 0.28 < xb < 0.3 with u = 5 GeV and 0.26 < o^b < 0.28 with 
u = 20 GeV, we estimate the orders of magnitude for the total cross sections, coming 
purely from the DVCS diagrams and compare them with the Compton contribution 
to the standard electromagnetic DVCS process. They are listed in Table |TV1 at the 
end of this chapter. 



VII. 3.1 Weak Neutral Neutrino-Proton Scattering 



The photon cannot be emitted by the neutrino, and thus in neutrino scattering from 
a proton via the exchange of Z^, we only measure the Compton part. The T- matrix 
of the process is given solely by the DVCS diagram, i.e. the diagram (a) in Fig. 



iT 



up 



u{k') 



cos 6 



w 



1 



\u{k) 



-i {qux - qiuQix/Mlo) 



Ml, 
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X 



cos 9 



w , 



(259) 



with T-^^ computed in the handbag approximation from the diagrams shown in Fig. 
^ The vector and axial vector couphngs of the vertex vvZ^ are Cy = = 1/2. 
Furthermore, since M|o ^ Ql and cos 9w = Mw/Mzo (recall that the Fermi constant 
IS Gp = 9^ I^^My^), the amplitude turns into 



- vp 



V2 |e| Gfu {k') 7, (1 - 75) u{k)tl (gs) U^n- 



Its spin-averaged square then reads 



- vp I 



1 ri2 T {v) TT^P 



The neutrino tensor is simply 

8 k^k'f^ + kpk'^ - g^p {k ■ k') + ieyp„rk"k"' 



(260) 



(261) 



(262) 



On the other hand, the weak neutral hadron tensor has a rather complicated expres- 
sion. Nevertheless, in the DVCS kinematics, it reduces to 



rrvlB 
J^WN 



1 



e+ 



t 



n 



WN 



4M2 



£ 



WN 



+ 



+ 



n 



WN 
2 



+ 



n 



WN 



+ 



n 



WN 



WN 



WN 



£ 



WN 



-2^ 3? (T^WN^WN ~^ '^WN^WN ~^ '^WN^WN ~^ '^WN^^ 



4M2 



+2 



-e 



4M2 



? Jv \^n.i4/;vC;y^ -t- c^^^n^^ -t- riy^j^o-i^j^ -t- c^^^n^^^vj 
1 



X 



(263) 



where the integrals of OFPDs are conveniently defined as 



/I 

( ) 

^WN (^5 ^) = X! Qf^V{A) J ^ 



(x — 4 + 20) 



-i(x-^ + iO) 
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FIG. 40: Weak neutral DVCS cross section a^p plotted as a function of the angle Os-y 
between the incoming weak virtual boson and outgoing real photon in the target rest 
frame for Q\ = 2.5 GeV^ and xb = 0.35 with tu = 5 GeV (sohd hne) and u = 20 GeV 
(dashed line) neutrino beam. 



dx 



-H 



+(-) 



dx 



-E 



+{-) 



x,^,t) 



(264) 



After the substitution of Eq. p61|) into Eq. ()119|) . the unpolarized differential cross 
section for the weak neutral DVCS process on a proton target using the neutrino 
beam takes the following form 



up 



1 aGl 



1 + Xb{M/uj) 



dxBdQfdtd^ {2ny 16 M^u^ [2 + {M/u)]xb [y + 2xb (M/u)] 



T M ttvI3 



(265) 



The cross section is illustrated in Fig. EOl 



VII. 3. 2 Weak Neutral Electron-Proton Scattering 

In the electron scattering, one has to include both the Compton and the Bethe-Heitler 
contributions, Tep = Tcep + ^BHep- The T-matrix for the Compton part can be easily 
obtained from Eq. ()259p by replacing the neutrinos with the electrons, 

Tcep = 2V2\e\GFuik')^,{c'y-^,c':,)uik)e;iq2)r^%. (266) 
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FIG. 41: Compton cross section acep for the weak neutral DVCS process pfotted as a 
function of the angle Ob-^ between the incoming weak virtual boson and outgoing real 
photon in the target rest frame for Q\ = 2.5 GeV^ and xb = 0.35 with u = 5 GeV 
(solid line) and a; = 20 GeV (dashed line) electron beam. 



The couplings are now equal to Cy = — 1/2 + 2 sin^ 9w and Cj^ = —1/2. Note that the 
hadron tensor is the same as in the neutrino case. Due to averaging over the spin of 
incoming electrons, the electron tensor has an extra factor of 1/2, 



-(e) 



Qup {k ■ k') 



k^k'^]. 



(267) 



The unpolarized Compton differential cross section for the weak neutral DVCS process 
with a proton target and an electron beam reads 

d^ocep _ 1 aG\ 1 + xb{M/uj) 



dxBdQldtd^ {2T^f 4 M'^u'^ [2 + {M /uj)]xB[y + 2xb {M /uj\ 



T (e) 
72^vl3^WNi 



(268) 



and it is plotted in Fig. |^ We notice that the cross section in the weak neutral 
current sector is larger with neutrinos than it is with electrons. This difference is a 
consequence of the structure difference in leptonic tensors given by Eqs. (j262p and 
(EHZD- 

Next we investigate the Bethe-Heitler contribution to the process. Since incoming 
and scattered leptons are electrons, we have to include both Bethe-Heitler diagrams. 
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see the diagrams (b) and (c) in Fig. 1211 The Bethe-Heitler amphtude can be imme- 
diately written down with the help of the expression p4()j) along with the following 
modifications: 

• The vertex ee7 is replaced by the vertex eeZ^, 
i \e\ Y -igY (cy - 75Ca) /2cos6'h/. 

• The photon propagator is replaced by the weak boson propagator, 



The coupling 7p is replaced by the coupling Z^p, —i\e\ 



-ig/ cosily. 



The proton electromagnetic transition current is replaced by the proton weak 
neutral transition current, J'^^ (0) J'^q (0). 



Hence 

iBHep = 2V2\e\GFe*{q2)u{k') 



7^^+ ^2)7"(4-75ci) 



7"(4-75C^)(^- #2)7^ 



u{k) (0) 



(269) 



where the transition current, 

jf^(o) = 1 [v;^'^ (0) - (0) 

1 
2 



\ E (0) 7.^/ (0) - c^Ai^f (0) 7.75^/ (0) 



(270) 



has both the vector V^^ and axial vector A^'~^ contributions. Their proton matrix 
elements are parametrized in terms of the four form factors |42j . 



(p(P2,S2)|l^j^^(0)|p(pi,Si)) = U{P2,S2) 
(p(p2,S2)|A^^(0)b(pi,Si)) = n(p2,S2) 



NC u\ 



2M 



uipi,si), 

u{pi,si). 
(271) 



Each of the vector form factors can be further expressed as a linear combination 
of flavor triplet, octet and singlet form factors (they correspond to the proton 
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matrix elements of the appropriate linear combinations of quark vector currents, 
(u'jiyU — d'ji.dj /2, (u'ji.u + d'ji.d — 2s^i,s^ and {u'-)i,u + ^71^0? + s^vS^, respectively), 



^1^2"?^ = (l-2sin2^^) 



^f(2) {t) + -Fi%) {t) 



- -A) it) 



(272) 



Moreover, the axial and pseudoscalar form factors are expressed as the difference 
between the isovector and strangeness form factors (they correspond to the proton 
matrix elements of current combinations, (u'jty'-f^u — d'-)y'^^d}j /2 and (571^755), respec- 
tively). 



G 



(273) 



If the sea quarks are neglected, they all simplify to 



(274) 



and 



9 A (t = 0) 



1 - 



2 ml-t' 



(275) 



For the unpolarized Bethe-Heitler differential cross section we have 
d^asHep 1 cyGp 1 + xb{M/uj) 



dxBdQfdtdip {2TTy 4 M^u;''[2 + {M/u;)]xB[y + 2xb{M/u;)] 



where the electron tensor is 



TUp ttBH 



(276) 



(7'^ i2 + 2A;'^) Y (4 - 75CI) , Y (4 - 75C^) (#27^ - 2A;'^) 



2 {k' ■ q^) 



+ 



2 {k ■ gs) 



X }i 



7^ (4 - 75C^) {421, + 2/c;) (^^ ^2 _ 2k,) 7^ (4 - 754) 



2 (fc' • ^2) 



+ 



2 (A; • gs) 



(277) 



and the hadron tensor is given by 



1 



H^p = -E(p(P2,52)|J.^^(0)b(pi,5i))(p(p2,52)|Jr(0)|p(Pi,^i)r- 



S\,S2 



(278) 
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TABLE III: Polar angles (f) and (j)' of the incoming and scattered leptons, respectively, 
in the target rest frame for Q\ = 2.5 GeV^ and xb = 0.35 with two different lepton 
beam energies uj. 



uj [GeV] 


<P [deg.] 


<P' [deg.] 


5 


10.3 


47.9 


20 


20.2 


25.3 



The angular dependence of the cross section is shown in Figs. 021 and EHl The two 
poles at the locations (j) and cj)' from Table 11111 are clearly seen. In Fig. El both 
contributions are plotted together. Contrary what one might expect from the results 
on the standard DVCS, the Compton cross section in the weak neutral interaction 
sector is larger, or at least comparable, with the corresponding Bethe-Heitler cross 
section for the scattering angles Oz-y < 8*^, which gives the region 0.15 GeV^ < -t < 
0.572 GeV^ in the invariant momentum transfer, see Fig. and hence the ratio 
-t/Ql < 0.23. 

VII. 3. 3 Weak Charged Neutrino-Neutron Scattering 

The invariant matrix element for the Compton contribution is 

Tcun = V2\e\GFu{k')-f^{l--f5)uik)e;iq2)T^c, (279) 

where the amplitude T^'^ is, in general, given by Eq. ()258j) . However, in our simple 
model the quark flavors / and / ' in Eqs. ()258|) become d and u, respectively, and the 
coefficients are therefore equal to (5+ = 1/6 and = —1/2. In the following, we 
relate the flavor nondiagonal GPDs to the flavor diagonal ones. This can be achieved 
by using the isospin symmetry. Through the isospin symmetry, the nucleon matrix 
elements, {p {p2, S2)\0'^'^^ (z) \n {pi, si)) and {p {p2, S2)\Of^^ (z) \n {pi, si)) given by 
Eq. fj256p that are nondiagonal in quark flavor, can be expressed in terms of the 



104 



d'o-aiplidcB^idtdp) [ib/GbV*] 



ezy [c3eg.] 



FIG. 42: Bethe-Hcitlcr cross section asHep for the weak neutral DVCS process plotted 
as a function of the angle Os-y between the incoming weak virtual boson and outgoing 
real photon in the target rest frame for Ql — 2.5 GeV^ and xb — 0.35 with uj — 5 GeV 
(solid hne) and a; = 20 GeV (dashed line) electron beam. 




0Zy [deg.] 



FIG. 43: Bethe-Heitler cross section asHep for the weak neutral DVCS process plotted 
as a function of the angle Os-y between the incoming weak virtual boson and outgoing 
real photon in the target rest frame for Ql = 2.5 GeV^ and xb = 0.35 with u) = 5 GeV 
(solid hne) and a; = 20 GeV (dashed line) electron beam. 
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d'a-/(dxgd^,ckdp) [rb/GsV*] 



0zy [cleg.] 



FIG. 44: Compton cross section acep (solid line) and Bethe-Heitler cross section asHep 
(dashed line) for the weak neutral DVCS process plotted as a function of the angle 
Ob'j between the incoming weak virtual boson and outgoing real photon in the target 
rest frame for Ql = 2.5 GeV^ and xb = 0.35 with lu = 5 GeV electron beam. 



flavor diagonal nucleon matrix elements |110j . Namely, 



{piP2,s,)\0''-^{z)\n{p^,s,)) = {p {p2, s,)\0^^ (z) \p {p,, s^)) 



u± 



{pip2,s,)\0'^{z)\pip^,s^)), (280) 



and similarly for {p {p2, S2)\ O^"^"^ (z) \n {pi, si)). Accordingly, the reduced weak 
charged VGA assumes the form 



A{p-q) 



{P ■ (I2) 



{p'q'2 + P'Q^) - 9 



X 



^WCU {P2, S2) 42U (pi, Si) + £^cU (p2, S2, 



{42 j- i i2 

4M 



■U (pi. Si) 



-'^WC^ iP2, S2) 42lbU {Pl, Si) + SwC ^^lZ^ ^ iP2, S2) I5U {Pl, Si) 



1 



ip-q2 



-ie^'''''q2pPn 



X 



2M 



(?2 ■ r) 



T-^wcU (P2, S2) 42I5U (pi, si) - S^c 2M ^ 
-T-(-wcU (P2, S2) i2U (pi, si) - E^c"^ (p2, S2) ^^^-^j^/^" '^^i' 
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(281) 



with the convolution integrals 

dx 



-i{x-i + zO) 

-1 {x-i + iO) 

-1 (x - ^ + iO) 

1 (ia; 

-1 (x - ^ + zO) 



Q+(_) (Et - Et) + Q-(+) (ez - E^) 



Q+(_) (e+ - Ej) + Q_(+) (e- - E, 



(282) 



The weak charged hadron tensor, 



'^■'■wc y-'-fiWc) 



t 

' UP 



n 



wc 



+ 



wc 



wc 

2 



+ 



WC 



'^WC 

t 



+ 



n 



wc 



4M2 



wc 



-2,^^3fJ (ji-wc^wc + '^wc^wc + '^wc^h/c + '^wc^^ 



+ 



wc 



wc 



+2 



{p-q2) ■ ' ■ (p- g2 

(l — 3^ ('^w^c'^VKc + '^wc'^wc 



V /9 



+ 4M2y ^^"^"^'^"^^'^ ^ 4M2 
<; ji \^ 'T-Tyc'^VKC "T ^vKc'T-vyc ''■VKC^Tyc ^wc'^wc ) 
1 



X 



(p ■ 92) 



^283) 



has identical structure as the expression ()263p for the weak neutral hadron tensor. 
Moreover, by neglecting the mass of the outgoing muon, one can substitute the so- 
called muon tensor l!]'^^ by the neutrino tensor L^^p, see Eq. (j262j) . The unpolarized 
Compton differential cross section for the weak charged DVCS process on a neutron 
target using the neutrino beam is 



d^acv 



1 aGl 



1 + Xb (M/lu) 



dxBdQidtdifi {2ny 16 M'^cu^ [2 + (M/lu)] xb [y + 2xb {M/cu)] 



T (m) Tjt^P 

2^vl3 ^WC- 



(284) 
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d'a-o,rJ{dxBd^,ddp) [rb/GeV*] 
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FIG. 45: Compton cross section acun for the weak charged DVCS process pfotted as a 
function of the angle Os-y between the incoming weak virtual boson and outgoing real 
photon in the target rest frame for Q\ = 2.5 GeV^ and xb = 0.35 with u = 5 GeV 
(solid line) and a; = 20 GeV (dashed line) neutrino beam. 



The cross section is presented in Fig. |33 

We turn our attention now to the Bethe-Heitler background. It is given only 
by the diagram (b) in Fig. EHl At the angle 0', when the outgoing real photon is 
collinear with the scattered muon, one should expect a pole. The amplitude of the 
Bethe-Heitler contribution reads 

V(r+ ^2)7^(1-75)' 



BHun 



V2|e|Gpe* {q2)u{k') 



u{k) 



X (P(P2,S2)| J,^^(0)|n(pi,si)). 
The spin-averaged square of Eq. ()285|) is 



■ BHun 



[285) 



(286) 



with tensors 



rZ//3 

J^BH 



ik' ■ q2) 



[k-q§ + kPq^, - g"^ {k ■ q^) - ze'^^'^^fc.ga. 



BH 

ul3 



E iP2, S2)\ J^"" (0) \n {p,, s,)) {p {p2, S2)\ J^^^ (0) \n {p,, s,))* • 



(287) 
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The matrix element of the weak charged transition current J^^ between the nucleon 
states is defined as 

(P (P2, S2)\ J^^ (0) \n (pi, Si)) = {p ip2, S2) I iJp (0) 7,^ (1 - 75) iJn (0) \n (pi, Si)) . 

(288) 

Using the same isospin symmetry relation between the fiavor nondiagonal and fiavor 
diagonal nucleon matrix elements, see Eq. ()280|1 . we parametrize the vector current 
part of the matrix element in Eq. ()288|) in terms of the Dirac and Pauli form factors 
for each quark fiavor, 

{p{p2,S2)\'ipp{0)lu'^n (0) |n(pi,si)) = {p{p2,S2)\lpu (0) lu^Pd (0) \n{pi,Si)) 

= {p{p2,S2)\'ipu (0)7i.V^« (0) \p{puSi)) - {p{p2,S2)\lpd {0)lu'ipd (0) \p{pi,Si)) 

{1(7 \ 1 
[Flu (t) - Fi, (t)] 7. - [F2U (t) - F2d (t)] I u (pi, si) . 

(289) 

Furthermore, for the axial vector current part we have 

{P (P2, S2)\i>p (0) -f^'y^ipn (0) \n (pi, si)) = u (ps, S2) {qa (t) 7^,75 

m(pi,si), (290) 



-9p [t) 



2M 

where the t-dependence of the axial and pseudoscalar from factors is given by 



gpit) = g^it)^- (291) 

Thus the unpolarized Bethe-Heitler differential cross section for the weak charged 
DVCS process is written in the form 

tLijdijIJ,, 



dxBdQ'idtd^ (27r)' 16 M^u^2 + {M/u)]xB[y + 2xb iM/u)f ' 

(292) 

Its dependence on 9wy is plotted in Figs. |iniand|17[ and together with the Compton 
cross section in Fig. |3H1 In particular, the plots in Fig. EHl reveal that, in the forward 
direction up to 4°, the Bethe-Heitler contamination is suppressed compared to the 
Compton contribution. 
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FIG. 46: Bethe-Heitler cross section asHun for the weak charged DVCS process pfot- 

ted as a function of the angle Os-y between the incoming weak virtual boson and 
outgoing real photon in the target rest frame for = 2.5 GeV^ and Xb — 0.35 with 
LU = 5 GeV (solid line) and ui = 20 GeV (dashed line) neutrino beam. 



To summarize, we gave a comprehensive review on neutrino-induced DVCS pro- 
cesses. Through the weak interaction currents, these processes give access to two 
different sets of GPDs, namely, to the sum of quark and antiquark distributions and 
to their difference. Thus we can measure independently both the valence and sea con- 
tent of GPDs. In addition, the weak charged current interaction probes GPDs that 
are nondiagonal in quark flavor, such as those associated with the neutron-to-proton 
transition. 
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FIG. 47: Bethe-Heitler cross section asHun for the weak charged DVCS process plot- 
ted as a function of the angle 9bj between the incoming weak virtual boson and 
outgoing real photon in the target rest frame for Ql — 2.5 GeV^ and xb — 0.35 with 
CO — 5 GeV (sohd line) and u; — 20 GeV (dashed hne) neutrino beam. 
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FIG. 48: Compton cross section acun (sohd line) and Bethe-Heitler cross section 
CBHun (dashed line) for the weak charged DVCS process plotted as a function of the 
angle Os-y between the incoming weak virtual boson and outgoing real photon in the 
target rest frame for Ql = 2.5 GeV^ and xb = 0.35 with uj = 5 GeV neutrino beam. 
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TABLE IV: Orders of magnitude (in nbarns) for the unpolarized total weak neutral 
(T;^p, Compton weak neutral crc^p, Compton weak charged Ccun and Compton elec- 
tromagnetic acep cross sections for two different kinematical regions of the Bjorken 
scaling variable with two different lepton beam energies cu. 



kinematical region u [GeV] a^p (Z°) acep (2"°) acun {W^) crcep (7) 



0.28 <xb< 0.3 



5 



10-11 



10-13 



10-11 



10-5 



0.26 <xb < 0.28 



20 



10-9 



10-10 



10-9 



10-=^ 
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CHAPTER VIII 
WIDE-ANGLE REAL COMPTON SCATTERING 

VIII.l INTRODUCTION 

Standard electromagnetic Compton scattering provides a unique tool for studying 
hadrons. We have demonstrated in the previous chapter how the process can be ex- 
tended into the weak interaction sector, in which the relevant Compton amplitude 
probes the hadronic structure by means of the quark electromagnetic and weak cur- 
rents rather than two electromagnetic currents, and promise therefore a new insight. 
In the lowest QCD approximation, these currents couple to photons and weak bosons 
through point-like vertices. Moreover, in a specific kinematical regime, featuring the 
presence of a large momentum transfer, the behavior of both the standard and weak 
Compton amplitudes is dominated by the light-like distances. As a result the am- 
plitudes are given in terms of the handbag diagrams. There are several situations, 
where the handbag contribution plays an essential role ^H] : 

• Both initial and final photons are highly virtual and have equal space-like vir- 
tualities. This situation corresponds to the forward VCA. Its imaginary part 
determines the structure functions of DIS. 

• The condition on photon virtualities may be relaxed in a sense that the ini- 
tial photon is still far off-shell but the final photon is real and the invariant 
momentum transfer to the hadron is small. The situation corresponds to the 
nonforward VCA. It is accessible through DVCS. 

• The configuration, in which both photons in the initial and final states are real 
but the invariant momentum transfer is large. The physical process correspond- 
ing to this situation is known as wide-angle real Compton scattering (WACS). 

The Compton amplitude of WACS is a subject of the study in the present chapter, 
which contains some of the formal results from Ref. |lllj . An efficient way to study 
this amplitude is again to use the light-cone expansion of the time-ordered product of 
two electromagnetic currents in the coordinate representation. At moderately large 
values of Mandelstam invariants (s, —t, —u < 10 GeV^), the dominant reaction 
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mechanism corresponds to the handbag contribution just as in DIS and DVCS. To 
the leading order, the amphtude is given by the s- and ^-channel handbag diagrams, 
see the diagrams (a) and (b) in Fig. IIHI On the other hand, for extremely large 
—t ^ 10 GeV^ the purely hard contributions involving two-gluon exchange, such as 
the diagram (c) in Fig. 0^1 and having power-law behavior (i.e. ~ (— t)^), become 
dominant. In particular, at large values of t, the handbag diagrams are accompanied 
by the Sudakov form factor S (t) ~ exp — Q;sln^(— t) /Sn that drops faster than 
any power of —t. At moderate t, however, the hard contributions are numerically 
suppressed by a factor (a^/Tr)^ c:^ 1/100, and will not be considered here. 

In Section VIII. 2, we start from the formal light-cone expansion in terms of QCD 
string operators (with gauge links along the straight line between the fields which, 
for brevity, we do not write explicitly). In addition to the twist-2 string operators, 
we also include now the operators that are algebraically of twist-3 but given by the 
total derivatives of twist-2 operators. One refers to these as the kinematical twist- 
3 contributions, in order to distinguish them from the so-called dynamical twist-3 
contributions involving quark-gluon operators. The latter cannot be reduced to total 
derivatives jHSl- Furthermore, for simplicity, we consider a spin zero (i.e. pion) target. 
For that reason, only the matrix element of a vector-type string operator is needed. 
Its spectral representation is given in terms of double distributions. It is worth noting 
at this point that even though we restrict ourselves to the pion, our approach can be 
used in more realistic situation involving the nucleon target |112j as well as can be 
applied to the time-like region, in particular, to study the exclusive production of two 
pions in the two-photon collisions |113j . In Section VIII. 3, we compute the Compton 
amplitude for the pion up to the twist-3 level, and inspect its transversality, in other 
words, the electromagnetic gauge invariance with respect to both initial and final 
photons. The physical amplitudes (or alternatively the photon helicity amplitudes) 
of the WAGS process are obtained by projecting the Compton amplitude on different 
the polarization states of physical photons. We estimate these amplitudes and their 
corresponding cross sections by taking a simple model for the double distribution of 
the pion. The results are compared with the QED calculation. 
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FIG. 49: Handbag contribution (a and b) and the configuration (c) with the exchange 
of two hard gluons. 



VIII.2 



TWIST DECOMPOSITION AND PARAMETRIZATION 



We recall from Section IV. 2 that the leading light-cone singularity of the time-ordered 
product of two electromagnetic currents is contained in the handbag contribution, 



2n^z^ 



J2 Q) {s^'^'^Of {z\0) + it'^P'^'^Oi^ |0)} 



(293) 



in which a highly virtual quark propagator connecting the two photon vertices is 
convoluted with the blob describing the long-distance dynamics. In Eq. fj293|) the 



tensor s^^""^ = g^^pg''P -\- g^^'^gp'^ — g'^'^g^ and the QCD bilocal operators with the flavor 
index / and one open Lorentz index are 



i,f {-z/2) ^^i,f {z/2) - i,f {z/2) 7,V^^ {-z/2) 
xPf {-z/2) 7^75^/ (z/2) + (z/2) 7^75^/ i-z/2) 



(294) 



The Taylor expansion of these operators in the relative distance z involves covariant 
derivatives and yields local operators that are not symmetric in their indices. Never- 
theless, it was shown in Refs. gllHS] that O^{z\0) and {z\0) can be expressed 
in a compact form in terms of the contracted operators {z) = z^O{^{z\Q) and 
Ol {z) = z^Ol^ {z |0) that produce only symmetric operators, and quark-gluon string 
operators. Namely, the fully deconstructed vector and axial vector operators take the 
following form 



dv I cos 



IVK 



+ y sm 



IVK 



[Z ■ V) Vr, - V^Zr, ' 



K 
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X 



X 



ijj {-vz/2) ^ijj {vz/2) - i,j ivz/2) {-vz/2] 



dv 




'iVK 




sin 




Jo K 




[ 2 \ 



i)f {-vz/2) /75V^/ {vz/2) + {vz/2) /75V'/ {-vz/2) 



+ (quark— gluon operators) 



5r] 



{z\0) 



dv < cos 



tVK, 



+ — sm 



[Z ■V)Vr,- V^Zr, 



K 



X 



X 



ijf {-vz/2) /75^/ {vz/2) + t/jf {vz/2) ji-f^i^f {-vz/2) 



/■I dv 


'iVK' 


/ — sin 




Jo K 


. 2 . 



i,f {-vz/2) {vz/2) - i,} {vz/2) {-vz/2] 



+ (quark— gluon operators) , 
where we use the notation v = 1 — d^^ = d/dz^ and 



K 



z-vy 



(295) 



(296) 



with T> = d/dX being the total derivative. Neglecting terms and ignoring quark- 
gluon operators (note that they correspond to the dynamical twist-3 contributions, 
which are not considered in our approach) the decomposition ()295|) gives jHH lEHl 



of{z\o: 



dv < cos 



iv {z ■ T>) 



+ — sm 



IV (z 



■V) 



Vr 



'ijf {-vz/2) ^ijf {vz/2) - ijf {vz/2) {-vz/2)\ 



1 dv 



{z-V) 



sm 



iv {z ■ V) 



tfjf {-vz/2) /75^^ {vz/2) + {vz/2) ^^^tfjf {-vz/2) 



oi{z\o: 



dv < cos 



iv {z ■ T>) 



+ — sm 



iv {z ■ V) 



{-vz/2) ^-i^iPf {vz/2) + {vz/2) ^-f^tpf {-vz/2) 



1 dv 



{z-V) 



sm 



iv {z ■ V) 



'ipf {-vz/2) {vz/2) - {vz/2) (-f 2/2)] . (297) 



For the pion, the matrix element of the contracted axial vector operator vanishes and 
accordingly, the matrix element of the time-ordered product (j293j) reduces to 



(tt (P2) I iT {JU {z/2) Jl^, {-z/2)} Ivr (pi)) = 
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dv I cos 



V [r ■ z] 



+ie 



1 dv 
Jo [r ■ z) 



sm 



V [r ■ z) 

2 

V (r ■ z) 



x{n{p2)\tPf{-vz/2) ^ijf{vz/2)-^f{vz/2) {-vz/2) {tt (p,)) . 



(298) 



We have used the fact that under the matrix element, the total derivative of the string 
operator turns into the momentum transfer, id/dX^^ — > r'^ (recall that r =pi — p2 is 
the overall momentum transfer to the pion target). 

The next step is to parametrize the nonforward matrix element of the twist-2 part 
of the contracted vector string operator in Eq. ()298|) . This is performed in terms of the 
spectral functions, namely, the double distributions. Thus one needs to construct the 
parametrization for the matrix element of {z) satisfying the harmonic condition 



{7,{p2)\Of (^)lvr(pO) 



twist— 2 



0. 



(299) 



In general, the matrix element can be regarded as a function of the variables {p ■ z), 
(r ■ z) and z^, where p = {pi + P2) /2 is the average pion momentum. The stan- 
dard way to parametrize it, is by the decomposition into plane waves, where one 
explicitly separates the {p ■ z) components from the (r ■ z) components, and hence 
use two spectral functions, ff{P,a,t) and gf{P,a,t). Alternatively, one can write 
the parametrization in the form 



(tt {p - r/2)| (z) \7T {p + r/2)) 



twist— 2 



-2% {z''d,)T{z,p,r) 



with 



J^{z,p,r) 



^ d(3 C da e'^^-'^ \hf ((3, a,t) + O (z^) 
-1 J-i+W\ ^ ^ ^ 



(300) 



(301) 



where the four-momentum k = j3p + ar/2. The latter has the same 

structure as that parametrizing the twist-2 part of the matrix element 

{p — r /2\ (j) {—z/2) (j) {z/2) \p + r /2) in the scalar case. In particular, it satisfies 
d'^J-'{z,p,r) = 0. It is easy to see that for that reason 



{z''d^)J^{z,p,r) 



i2 + z''d„)d^J^{z,p,r) 
0. 



(302) 



Let us note at this point that in Appendix C, we consider the scalar parametrization 
JF {z,p, r), which includes terms up to order O (z^). The inclusion of z"^ and z"^ terms 
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in the expansion is due to the fact that the fermion propagator has Xj singularity, 
and besides the operator (f z) in Eq. ()298j) is differentiated with respect to z. 
Now, by taking only the first term in the z^ expansion in J-'{z,p,r), the expression 
PIH]) turns into 



(tt {p-r/2)\Of (z) |7r(p + r/2)) 



twist— 2 



df3 



1-1/31 



1+1, 



da e^^'^-") 



X {k ■ z) hf {(3, a, t) . 



(303) 



The representation ()303|) is given in terms of only one double distribution hj {/3,a,t) 
and corresponds to the so-called single-DD parametrization |88j. We observe that the 
original distributions are then expressed in terms of the single one as 



ff{p,a,t) = phf{P,a,t), 
gf{(3,a,t) = ahf{(3,a,t). 



(304) 



Moreover, we establish the symmetry properties of hf{(3,a,t) with respect to pa- 
rameters P and a. Since // {P,a,t) is odd in /3 and even in a while for gj {P,a,t) 
the situation is the other way around, it follows from Eq. ()304p that the double 
distribution hf {(3,a,t) is an even function in both variables f3 and a. 

Finally, by substituting the single-DD parametrization (|303p into the matrix ele- 
ment (|298p . one arrives at 

(vr (p2)NT { J^^, iz/2) r^,, {-z/2)} Ivr {p,)) = 



HQ] 

f 



dp 



1-1/31 



1+1/31 



da hf {p, a, t) / dv e^^'^'") 









V {r ■ z) 




COS 








2 



[kjj + ivkr, {k ■ z)] + -r^ [k ■ z) sin 



v[r ■ z] 



(r ■ z) 



sm 



v{r ■ z) 



k^ + ivk^ {k ■ z) 



(305) 



To isolate the leading twist-2 contribution in Eq. ()305|) . we integrate by parts over v 
the first fc^ term in the symmetric part (with respect to the Lorentz indices and v) 
of the matrix element, and the first k^ term in the antisymmetric part. We find 



(tt {p,)\tT {JP,, {z/2) r^,, {-z/2)} In {p,)) = 

Z /"l /•1-I/3I 



X I s^^'^^ I kr,e'^^-'^ + ^ [r^ {k-z)- kr, {r ■ z)] £ dv t;e™(^ ") sin 



V [r ■ z] 
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+U'"''"'e„5Kxz^r^k^ C dv ve'"^''-'^ cos 
2 7o 



V {r ■ z) 



(306) 



where only the surface term in the symmetric part corresponds to the twist-2 con- 
tribution (note that the surface term in the antisymmetric part vanishes). All the 
remaining terms are of the twist-3. We explicitly write 



(vr {p,)\iT {r^,, {z/2) r^,, {-z/2)} |7r {p,)) 



twist— 2 



Zp 



Qjsf'""'^ J^^ d(3 j^^ da hf {(3, a, t) kr^e'^^''^ (307) 



Tl'^Z^ ^ ' J -I "'-1+1/3 



for the twist-2 part. Now for the twist-3 part, after expressing the sine and cosine 
functions in Eq. ()30(i|l in terms of the exponentials and combining them with the 
overall exponential factor exp [iv {k ■ z)], we obtain 



/ 

X \ s^''^" [r^ (k-z)- kr, (r ■ z)] I dv v 



TT {p2)\tT {r^M {z/2) J^EM (-^/2)} k (Pl)),^,,,_3 = 



IZp 







^iv{ki-z)—i{r-z)/2 ^iv{k2-z)+i(r-z)/2 



where ki = k + r/2 and ^2 = k — r /2 are the momenta of the incoming (with respect 
to the interaction with two photons) and outgoing quarks, respectively. 

VIII.3 COMPTON SCATTERING AMPLITUDE 

As a result of the twist decomposition of the matrix element, see Eqs. (|307|) and 
), the Compton scattering amplitude. 



' f 



J d'z e^^"-) {7r{p,)\iT{JUAzmrEM{-z/2)}\7^{p,)), 

(309) 



where q = {qi + q2) /2 is the average photon momentum, can be written in our 
approximation, as the sum of the twist-2 and twist-3 contributions, 

^ — ■'■twist-2 ^ ■'■twist-3- K'^^^) 
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VIII.3.1 Twist-2 Amplitude 

After integrating over z with the help of Eq. the twist-2 part of the amplitude 

reads 



J-t, 



twist— 2 



dp 



1-1/31 



^ ■' J ■n^z^ 7-1 

X {2A;^r + A;^g^ + Fg'^ - g^"" + (A; ■ g)] } 



da hf {P, a, t) kr, 



(311) 



If the pion mass is neglected then k"^ = — (/?^ — a^) t/A and the denominator for 
the quark propagator results into 



{k + qf = {k, + q,f 

= {k2 + 



(312) 



where s = 2{p-q) — t/2 = 2 (pi • gi) = 2 (p2 • ^2) or, after writing (1 — /?)^ 
(1-1/31)^ + 2(1/31-/3), 



(fc + g)' 



ps-[{l- p'f -a'^]t/4 for/3>0 
-/?M - [(1 + Pf - a^] t/4 for /3 < 



(313) 



VIII.3.2 Twist-3 Amplitude 

We compute now the twist-3 amplitude, which contains two contributions. First, the 
contribution from the s^^"^^ term in Eq. ()308|) is 



■^twist-3{l) 



1 



2 „tipvr] 



X J dv V I d'^z 



1 .1-1/31 

dp I da hf {P, a, t) 

1 "'-1+1/31 



tZr 



2ti^z^ 



[r^ {k- z)- kr, (r ■ z) 



(314) 



where the combined exponentials, /i = gi + vki and I2 = q2 + vk2, correspond to 
the effective momenta flowing through the quark propagator, and gi = g — r/2 and 
q2 = Q + f/2 are the momenta of the initial and final photons. We have two types of 
quark propagators here. Their denominators are given by 



{ki ■ gi) + vkl 
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{k + qY — vk\ 
{h ■ qi) + vkl 
{k + qf — vkl 



where the incoming and outgoing quark virtuahties are 

kl = + (1 + t/4, 



^2 



-P^ + (1 - a)^] t/A. 



After the integral over z is carried out using of the formula 

J 27r2 (z2 _ ^0)^ 

the expression ()314|) turns into 



(/2 + iO) 



2 ' 



1-1/31 



/ / da hf {(3,a,t) / u 



-1 ^-1+1/3 



(315) 



(316) 



(317) 



X 1^ [- (/i • A:) (/fr'^ + l\rn + (/i ■ r) (/^^r + /^fc^)] 
[- (^2 ■ A;) (/^r'^ + /^V) + (/2 • r) (/^^F + /^^F)]! 



(318) 



Next we study the second contribution to %^ist-3- It is due to the e^'''^^ term in the 
matrix element ()308t and comes from the axial vector operator. We have 

/•1-I/3I 



\Y.Q) t dp da hj {P,a,t) 

2 ^ ' J-i J-i+\l3\ 



xe^^'^f'erjSKxr^k^ f dvv f d^z 



Jih-z) I i{l2-z) 



IZpZ 
2^ 



da hf (3, a, t) I dv v 

-1+1/31 Jo 



+ Y, [{h ■ k) {l',r^ - r,rn - {k ■ r) {ll^k'' - F.k^] 



Finally, by adding Eqs. 



■if. 



twist— 3 



and 

■1-1/3 



the total twist-3 result reads 



(319) 



rl rl-\P\ /-l 
2yQ] d(3 dahf{l3,a,t) dvv 

J-i J-l+l/31 Jo 

X U [- (/i ■ k) l\r^ + (/i ■ r) l\k^] + i [(/2 ■ k) l^r'' - {k ■ r) l!^k''] 
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2 y Q] d(3 da hf (/?, a, t) 



/ 

X <; / dvv^ 




rl 

+ / (if f 



^ [- (/i ■ A;) + (/i ■ r) F] + ^ \{h ■ A;) - (Z2 ■ r) 







f [- (Zi ■ A;) r'^ + (Zi ■ r) A;^] + f [(/s ■k)r^- {k ■ r) 



with the following scalar products 



(320) 



Hi-k) = {k + qf /2 + {v-l/2)kl + {y-va)t/A, 

(k-k) = {k + qf /2 + {v- 1/2) kl + {v + va)t/4, 

(l^.r) = -{v-va)t/2, 

ik-r) = {v + va)t/2. (321) 

Note that the second integral over v in the amplitude ()320p is logarithmically di- 
vergent, however, the coefficient for the divergent integral vanishes, when %'^ist-3 is 
contracted with the incoming and outgoing photon polarization vectors, (gi) and 
e^. (g2)- 

VIII. 3. 3 Electromagnetic Gauge Invariance 

The total Compton scattering amplitude ()309|) is transverse with respect to incoming 
and outgoing photon momenta 

r^'qiu = 0, 

52;.^'^" = 0. (322) 



Let us check whether the gauge invariant condition ()322|) also holds for the ampli- 
tude ()310p . We expand the contracted twist-2 and twist-3 amplitudes in terms of t/s 
and, for simplicity, keep only the lowest (zeroth) order corrections. A straightforward 
contraction, together with imposing the symmetry properties of the double distribu- 
tion hf{P,a,t) with respect to the variables a and /3, yields for the leading-twist 
amplitude 

.1 i-l-\l3\ 



d(3 / dahf{(3,a,t)r^/2, 

rl /-l-l/^l 

q2>.TZst-2 = d(3 dahf{l3,a,t)r-/2 (323) 

J-i J-i+\0\ 



122 



while for the contracted next-to-leading twist amplitude we get 



/i /"J-— IPI 

df3 / dahf{(3,a,t)r^/2, 
J -1 "'-1+1/31 

Q2,TZst-3 = J2Q} t d/3 da hf a, t) r72. (324) 

f "'-1 "'-i+Ls 



/ 

Clearly, even at the level 0{t/sf, the amplitude Tl^^g^_2 alone is not gauge invari- 
ant. However, the electromagnetic gauge invariance at this level can be recovered by 
including the contribution from T^l^isi-s- The latter exactly cancels Eq. ()32H|) . 

VIII. 3. 4 Helicity Amplitudes 

We shall project the Compton scattering amplitude onto the polarization states eiy = 
(gi) and = (^2) of the physical photons. In the Compton-like processes, these 
two states can be written as the linear combination of two basic polarization four- 
vectors |114j . 



n 



712 



A 

(325) 



-nl 



where 



is in the reaction plane and 



= p^-^^lSlq^ (326) 



= e^^^Pp^r^qp (327) 

is perpendicular to it. Then, by construction {q ■ ei) = (r ■ ei) =0 and similarly for 

€2- The contracted twist-2 amplitude is 

ri ri-\P\ 1 
%w^st-2^*2^,^lv = '^Y.Q) dp da hf{P,a,t)-——^ 

^ f -'-^ J-^+m (k + q) 

x{2(3' {p.e,){p-e;)-{e^-e;) + {k ■ q)]} 



f 



X 



-1 + 



iP^ -a^-1) t/A 



2 



(ei ■ ^2) + Tr^-T2 (p ■ ei) (p ■ £2) 

[k + q) 



(328) 
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and for the twist- 3 part we find 

V + va 



X / dv 
Jo 



V — va 



+ 



{k + q) — vk\ {k-\- q) — 



(329) 



Furthermore, we introduce twist-2 and twist-3 photon hehcity amphtudes by pro- 
jecting the photon polarizations onto the basic vectors ei and 62, 



H 
H 



(1-1) 

twist-2{3) 



(1,2) 

twist-2{3) 



twist- 2{3)^lH^'^'^^ 
-'imst-2(5)^lAi^2i/, 



H 



(2,1) 

twist-2{3) 



t-r-pi,V * 



rj-{2,2) 

^■^ twist- 2 {3) 



-^twist-2(3)^2tj.^2v 



Out of eight hehcity amphtudes, only three are nonzero, namely. 



H 



(1,1) 

twist— 2 



f 



d(3 



rl-l 
-1+ 



da hf (/3, a, t) 



a 



l)t 



+ 



us 



H 



(2,2) 
twist— 2 



H 



(1,1) 

twist— 3 



{k + qf 4 ' {k + qf t 

Vg^. / d(5 dahf(p,a,t) 
f J-i J-i+Wl 

{(3^ -a'-l) t\ 
{k + qf 4/' 

/•I /-i-l/SI 
-us y Ql / d/? da hf (p, a, t) 



f 

xll- 



Jo 



dv < 



V — va 



V + va 



{k + qf — vki {k + q) 



vkn 



where we have used 



(P • eiY 



t / t \2 

= 4-7(^ + 2;) 

us 

~ T' 



(330) 



(331) 



(332) 
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VIII.3.5 Model 



For the pion double distribution hf {j3,a,t), we use a simple model in a factorized 
form |115j . 

1 



f{(3)h{(3,a)T{(3,t). 



The profile function is chosen to be infinitely narrow 



(333) 



h a) 



6(a) 



(334) 



Accordingly, the integral over a in the helicity amplitudes ()331|) is trivial. Next, for 
the forward distribution we take 



31- 

4 



and the t-dependence is described by the following Regge ansatz 



(335) 



(336) 



where the parameter a' = 1 GeV^^. The summation of quark flavors in Eq. 1)3311) 
gives an overall charge factor of 5/9. 

We calculate separately the s- and u-channel contributions in the helicity am- 
plitudes. Using Eq. ()313|) and making the change in the variable /3 — >■ — /? in the 



w-channel, we have for the helicity amplitude H^^'^^ 



twist — 2 



H, 



(1,1) 

twist— 2 



1 + 



9 Jo pi ^/p 

iP' - 1) \t\ /4 



Ap^us/ \t\ 



Ps + {l-pf\t\/4 Ps + {1- Pf\t\/A\ 



^9 Jo ^ PA ^ ^ 



1 + 



Ap'^us/ \t\ 



iP' - 1) \t\ I A 

pu\{\- pf\t\lA Pu\{\- pf\t\lA_ 



(337) 



and after some manipulation 



i^t-2 = lll'dl^ a-f3)P 



~l/2+a'\t\{l-t3) 



X 



{2a'^-l)p + a {2a^-l)P-a 
p^ + 2ap + l ^ /52 -2ap + l 



(338) 
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where the parameter a = 2s/ |t| — 1 > 1 depends only on the ratio s/ \t\. The first 
term in the square brackets of Eq. ()338j) corresponds to the 5-channeI contribution. 
For the moderate invariant momentum transfer, the term has no singularities in the 
region < /3 < 1, in other words, the 5-channel contribution is real. The situation 
is, however, different with the second term in the square brackets, which represents 



the u-channel contribution. Its denominator has the roots /5i_2 = a ± — 1 (in 
fact, the denominator is equal to — 2a/? + 1 + ie = {P — Pi + ie) {P — P2 — i^)) 
and, unlike the s-channel, there is one relevant pole at /? = P2 + The u-channel 
contribution has therefore both the real and imaginary parts. They can be evaluated 
in the complex plane by making use of the residue theorem in the following way: one 
introduces a new variable, 7 = /3 — 1/2, which shifts the integration region from [0, 1] 
to [—1/2, 1/2]. Then by closing the contour in the lower half-plane we find 

r.1/2 



/■I/2 1-2-K 

/ ^ ^7 {...} = / ^7 {■••} 

J~l/2 Jtt 



(339) 



with 7 = (1/2) exp (i0). In summary, the helicity amplitudes H^tmil-s both the 
real and imaginary parts, where the latter is generated solely by the corresponding 
M-channel contribution. In the similar way, we compute the remaining two helicity 
amplitudes. 



H 



(2,2) 
twist— 2 



p + a 



+ 



-l/2+a'\t\{l-(3) 

P — a 



P^ + 2aP + l P^-2aP + l 



(340) 



and 



H 



(1,1) 

twist— 3 



^3(a2-l) J'dp {1-P)P 



l/2+a'\t\{l-f3) 



1 



{p^ + 2ap + if ^ (/32 - 2ap + 1)". 



(341) 



For the sake of simplicity, we have neglected terms and in the denominators of 
see Eq. (13311) . Then the integral over v simply gives 1/ {k + g)^. At present. 



(1,1) 

twist— 3^ 



we only consider the twist-2 helicity amplitudes. In Figs. EOl and |^ their real and 
imaginary parts are plotted against the center-of-mass scattering angle 6cm rather 
than against —t (recall that t = — s sin^ (6'cm/2)) in the region of sufficiently large 
invariant s = 5 GeV^. We divide the amplitudes by 5/9 to remove the overall charge 
factor. Analogous plots are presented in Figs. |S21and|SSlfor s = 10 GeV^, and in Figs. 
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FIG. 50: Real part of the helicity amplitudes Hljjf_2 (solid line) and Hj:jjf__2 (dashed 
line) plotted as a function of the center-of-mass scattering angle Ocm for the invariant 
s = 5 GeV^. The amplitudes are divided by the charge factor 5/9. 



l5ll and l55l for s = 20 GeV^. Due to the presence of the term 4/3^ (us/t) / {k + in 
the amplitude Hj.l^l]^_2, one might naively expect that the real part of would 
dominate over the real part of -f^t^'jsj_2- Both the s-channel and the real part of the 
w-channel contributions in -fft^'/st-2 ^^-^^ indeed large, however, they come with the 
opposite signs, and as a result, the real part of the amplitude h[1^1}^_2 diminishes, as 
oppose to the amplitude Hlljf}^_2j in which both terms add up. Another observation 
is that the imaginary parts of the helicity amplitudes are suppressed with respect to 
their real parts. 

We can now estimate the polarized differential cross sections in our simple model 
for the pion DD hf {i3,a,t). In the center-of-mass frame, we write 

1 



2Ei2uJiVj:cl 

dp2 



(27r)^5(4) {p, + q,-p2-q2) 



dq2 



(342) 



(27r)^2E2 (27r)^2cj2 

Note that we multiplied the helicity amplitude with e^, which has been conveniently 
skipped in our analysis. Partially evaluating the integrals over the final-state mo- 
menta, e.g. integration over q2 kills S^^^ (pi + qi — P2 — ^2) and sets q2 = —p2, yields 
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FIG. 51: Imaginary part of the helicity amplitudes H^l]^_2 (solid line) and H^]^_2 
(dashed line) plotted as a function of the ccntcr-of-mass scattering angle 9cm for the 
invariant s = 5 GeV^. The amplitudes are divided by the charge factor 5/9. 



Be 



FIG. 52: Real part of the helicity amplitudes h\^"^^ 



twist— 2 

(sohd line) and 



hue) plotted as a function of the center-of-mass scattering angle 9cm 
s = 10 GeV^. The amplitudes are divided by the charge factor 5/9. 



2 (dashed 
for the invariant 
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FIG. 53: Imaginary part of the helicity amplitudes H'^1]^_2 (solid line) and H\yj^g^_2 



(dashed line) plotted as a function of the center-of-mass scattering angle 9cm for the 
invariant s = 10 GeV^. The amplitudes are divided by the charge factor 5/9. 




FIG. 54: Real part of the helicity amplitudes H^^'^^ 



twist— 2 

(solid line) and //£fj 



line) plotted as a function of the center-of-mass scattering angle 9cm 
s = 20 GeV^. The amplitudes are divided by the charge factor 5/9. 
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FIG. 55: Imaginary part of the helicity amplitudes tl\wist-2 (solid line) and ntwist-2 
(dashed line) plotted as a function of the center-of-mass scattering angle dcm for the 
invariant s = 20 GeV^. The amplitudes are divided by the charge factor 5/9. 



at high energies (and thus frei = 2) 



1 1 1 



2s4(27r)^ 



b {2Ei - 2E2) dE2dn 



g2^(Ai,A2) 



(343) 



where s = E^^ and Ei = i?cm/2. The differential cross section can be further simpli- 
fied into 



dn 

which can be written in the invariant form as 

^^(Ai,A2) g4 



QAtt'^s 



(344) 



(345) 



dt IGvrs^ 

The 6'cm- dependence of the twist-2 cross sections cr*^^'^^ and a'-^'^^ are illustrated on a 
logarithmic scale for all three values of the invariant s in Figs. EHlandEZl respectively. 
Similarly to the situation with the nucleon target (see, e.g. p58j), the plots reveal a 
slight increase in both polarized cross sections at large scattering angles, dcm > 130°. 
Moreover, by comparing their orders of magnitude, we find that the cross section 
a^^'^^ is significantly suppressed compared to a*^^'^-*. Next, the angular dependence of 



the combination s^'^fia'-^'^^ / dt is shown in Fig. 
i.e. they basically coincide up to 6cm — 100°. 



The curves obey scaling behavior. 
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FIG. 56: Twist-2 cross section a^^'^^ plotted as a function of the center-of-mass scat- 
tering angle 9cm for the invariant s = 5 GeV^ (bold solid line), s = 10 GeV^ (solid 
line) and s = 20 GeV^ (dashed hne). 




6U[de&] 



FIG. 57: Twist-2 cross section cr*^^'^^ plotted as a function of the center-of-mass scat- 
tering angle 9cm for the invariant s = 5 GeV^ (bold sohd hne), s = 10 GeV^ (solid 
hne) and s = 20 GeV^ (dashed hne). 
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FIG. 58: Combination s^-^da^'^''^y dt plotted function of the center-of-mass scat- 
tering angle 9cm for the invariant s = 5 GeV^ (bold solid line), s = 10 GeV^ (solid 
line) and s = 20 GeV^ (dashed line). 



VIII.4 COMPTON SCATTERING ON A PION IN QED 



We compare the QCD results of Section VIII. 3 with the QED calculation. Assuming 
that vr"*" is an elementary particle we have, in the lowest order in uem-, three relevant 
Feynman diagrams, see Fig. EHl namely, the s- and n-channel tree diagrams and the 
four-point contact interaction diagram. Adding them coherently, one finds for the 
T-matrix 



iT 



i \e\ {2p2 + ■ ^2 T 

[Pi + Qi) - mi 



+i\e\ (2p2 
+2i\efg^' 



{P2 - qif - ml 



i \e\ {2pi + qiY 
-i \e\ (2pi - 



(346) 



The expression in the brackets corresponds to the Compton scattering amplitude T^'^ 
in QED. We can immediately check its gauge invariance and find 



r,u, _ _, . / i2p2 + q2r[2 {pi ■ qi) + qf] ^ {2p, - q^f [2 {P2 ■ Qi) - ql] 



2(pi-gi) 

-2%e^ {P2 + q2-Pi- qiY 
0. 



-2 (Pi ■ qi] 



-2gr| 
(347) 
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FIG. 59: Feynman diagrams for Compton scattering on a pion in QED. 



Similarly, we prove that q2iiT^^ = 0. 

Neglecting now the pion mass in Eq. ()34(ij) (recall that s, —t, —u ^ and hence 
s + t + u = 2ml — 0), we get 



-2e' 



2 (Pi ■ ei) {p2 



(348) 



By writing the initial and final pion momenta as = p + r /2 and p2 = p — r/2, 
respectively, and imposing (r ■ ei) = (r ■ Cg) = 0, we can express the T-matrix in 
terms of the average pion momentum. 



T = -2e2^2(^- + -j(p.ei)(p-e;)-(ei-e;)J. (349) 

Finally, the projection of polarization vectors ei and €2 onto ei and 62 gives the 
following photon helicity amplitudes in QED 



^(1,1) 

^(1,2) 
^(2,2) 



-2e' 

0, 
0, 

-2e' 



1 1 
- + - 

s u 



us 

T 



+ 1 



(350) 



The QED helicity amplitudes if^^'^^ and i/^^'^) are opposite in sign but constant. By 
comparing them to the real parts of the QCD helicity amplitudes -f^tVist-2 Hf:JJ^_2 
(and, as usual, ignoring — e^), we can establish the effective form factor, which solely 
describes the underlying QCD dynamics. 

In summary, we presented an improved treatment of the DD formalism to the case 
with the pion target. In addition to deriving the amplitude in the twist-3 approxima- 
tion, the helicity properties were considered. We gave predictions for the polarized 
cross sections with the use of the simple model for the pion double distribution. 
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CHAPTER IX 
CONCLUSIONS AND OUTLOOK 

In this dissertation we have studied inclusive and exclusive Compton processes in 
quantum chromodynamics. Let us summarize our main conclusions. 

First, we briefly reviewed the crown example of the hard (light-cone dominated) 
inclusive scattering process in QCD known as deeply inelastic lepton-nucleon scatter- 
ing. We introduced the forward virtual Compton scattering amplitude and discussed 
its relation to DIS. We presented the cross sections results in terms of structure func- 
tions for both the electromagnetic and weak probes. The last part of Chapter II 
was devoted to the QCD parton model. This simple model, which naturally emerges 
from the QCD operator product expansion, successfully describes the scaling prop- 
erty of structure functions, and further allows to express them in terms of the parton 
distributions functions. 

In Chapter III, we started with the systematic description of old phenomenologi- 
cal functions (form factors, parton distribution functions and distribution amplitudes) 
that were used for years in the studies of hadrons. New phenomenological functions 
(generalized parton distributions) became necessary in order to generalize the Comp- 
ton scattering amplitude in the region of nonforward (skewed) kinematics. They 
were introduced in various ways and under different names as nonperturbative func- 
tions describing the soft (long-distance) part of the factorized scattering amplitudes 
in QCD. Generalized parton distributions accumulate the most complete informa- 
tion about the hadronic structure, and as such combine the features of form factors, 
usual parton distributions and distribution amplitudes. We presented some of their 
theoretical aspects. 

Generalized parton distributions have numerous applications to hard processes, in 
particular, to the exclusive ones. The pedagogical example in this respect is deeply 
virtual Compton scattering. In Chapter IV, we demonstrated the derivation of the 
amphtude at the leading twist-2 level in the DVCS kinematics, using the light-cone 
expansion in terms of QCD string operators in coordinate space. We investigated 
the electromagnetic gauge invariance of the amplitude. The latter is violated with 
respect to the initial virtual photon and therefore, higher twist corrections have to 
be taken into account. Furthermore, we discussed the kinematics and introduced a 
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toy model for the nucleon off-forward parton distributions, which did not include the 
contribution from the sea quarks. We estimated the unpolarized cross sections in 
the kinematical region relevant to Jefferson Lab, and observed that unfortunately the 
contamination from the Bethe-Heitler process predominated the pure DVCS signal. 
The way out is to exploit the interference between the two processes. This approach 
allows to project out independently both the real and imaginary parts of the DVCS 
amplitude and probe different linear combinations of GPDs. 

In Chapter V, we studied the inclusive photoproduction of massive lepton pairs. 
The process was considered in the framework of the parton model at rather high pho- 
ton beam energy. Again one deals with two types of subprocesses. We picked only 
the lowest order electromagnetic contributions and calculated both the Compton and 
Bethe cross sections. We found that in the forward direction the Compton contri- 
bution was slightly larger than the Bethe-Heitler contribution, however, it decreased 
much faster. In addition, we illustrated that the interference terms cancel in pairs, 
after being integrated over the momenta of final leptons. 

The inverse process to DVCS is known as time-like Compton scattering. It can 
be accessed through the photoproduction of a heavy lepton pair. In Chapter VI, we 
considered this particular reaction in the DVCS kinematics. The unpolarized cross 
section results revealed that the Compton contribution was significantly suppressed 
with respect to the Bethe-Heitler contribution. In analogy with the DVCS process, we 
expect that the interference terms between the Compton and Bethe-Heitler processes, 
in particular, from the experiments with the polarized photon beam, may play an 
important role in extracting new information on GPDs. Thus more future studies are 
needed in this direction. 

By utilizing the weak currents to probe the nucleon structure we are able to mea- 
sure, due to the V — A nature of interaction, a different combination of generalized 
parton distributions as well as the distributions that are nondiagonal in quark fla- 
vor. Hence neutrino-induced virtual Compton scattering, which has been studied in 
Chapter VII, provides an important tool to complement the study of GPDs in more 
familiar electron-induced DVCS or exclusive meson production processes. We derived 
the twist-2 amplitudes for the weak neutral and weak charged current interactions, 
and gave predictions for the cross sections in the kinematics relevant to future high- 
intensity neutrino experiments. At small scattering angles, we observed that, unlike 
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the standard electromagnetic DVCS process, the Compton signal was enhanced com- 
pared to the corresponding Bethe-Heitler contribution. Prom the theoretical side, it 
is expected to use a more realistic model for the nucleon GPDs, and further check 
separately the contributions from the plus and minus distributions, in particularly, 
the contribution coming from the pion pole, i.e. the Ef distribution. Moreover, 
one may estimate the contribution from the interference of the Compton with the 
Bethe-Heitler process, and extend the approach in order to include the twist-3 terms. 

In the final chapter, the GPD formalism was applied to real Compton scattering. 
Again we used the same light-cone expansion technique as in Chapters IV and VII, 
however, in addition to the twist- 2 operators we have also included a set of kinematical 
twist-3 operators. They appear as total derivatives of twist-2 operators. Unhke the 
previous chapters, we considered, for simplicity, the pion target. Instead of off-forward 
parton distributions, more general objects, namely, the double distributions were used 
to describe the nonperturbative stages of interaction. We obtained the expressions for 
the twist-2 and twist-3 Compton scattering amplitudes. It was shown, in the lowest 
order in the invariant momentum transfer, that the kinematical twist-3 contribution to 
the amplitude for pion was required to restore the tranversality. Next we found that a 
divergent part of the twist-3 amplitude had zero projection on the polarization vectors 
of initial and final photons, and calculated the photon helicity amplitudes, using a 
simple model for the pion double distribution. They were compared to QED helicity 
amplitudes. Finally, we estimated the relevant polarized cross sections, observed 
their approximate scaling behavior, and established the noncanonical powers. A 
straightforward generalization to the case with the nucleon target as well as the 
study of the exclusive production of two pions in two-photon collisions (the latter is 
a crossing process to Compton scattering on a pion) are in progress. 
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APPENDIX A 
GLOSSARY 

The following abbreviations are often used in a text: 

DA distribution amplitude 

DD double distribution 

DIS deeply inelastic lepton scattering 

DMP deeply exclusive meson production 

DVCS deeply virtual Compton scattering 

GPD generalized parton distribution 

OPE operator product expansion 

PDF parton distribution function 

QCD quantum chromodynamics 

QED quantum electrodynamics 

TCS time— like Compton scattering 

VGA virtual Compton scattering amplitude 

WACS wide — angle real Compton scattering 
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APPENDIX B 

LIST OF INTEGRALS AND SCALAR FUNCTIONS 
THAT APPEAR IN THE BETHE-HEITLER 

SUBPROCESS OF THE INCLUSIVE 
PHOTOPRODUCTION OF LEPTON PAIRS 

In Section V.3, we compute the double integral ()182|1 . The calculation is carried out 
through the following set of integrals: 

„,.„_1_ . 

+h {q^q" + q'q'n + hg"', 

I'''''^^''^ Uq.q')-ik.q) ^ 
r kf^ 

r k^kP 

d'k6+06+{)- = J.q^qP + J.q'^q'" 

J [q ■ q') -{k-q) 

+J5{q^q"' + q'q'n + Jeg^', 

/c^^fc^+O^O \ = -^(/o + Jo), 

J [k-q) [{q ■ q') - [k ■ q)\ {q ■ q') 

J [k-q) [{q ■ q') - [k ■ q)] 2 [q ■ q') 

I'''''^^'' ^k.qnq.':tik.q)^ ^ ^ ' 

I'''''^^''^ \q.q'-lq) ^ " ^'^ 

fd'kS^OS^O = K,q^ + K,q'P, 

J [k-q) [{q ■ q') - (A; ■ q)\ 

(d'k5^{)5^{) /r^^^lw, ,1 = K.q^^q" + K.q'Pq'P 
J [k-q) [{q ■ q') - [k ■ q)\ 

+K,{q^q'^ + q^q'n + K,g^', 
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d^k 5+ 5+ 

d^k 5+ 6+ 



1 



kf" 

{k-qf 
kt'kP 

{k-qf 



(k ■ q) [{q ■ q') - (k ■ q)] 
k^'kf 

{k ■ q) [{q ■ q') - {k ■ q)] 



L,q>^ + L,q"^, 
L.q'^qP + Uq'Pq'P 

Mig^ + M^q'^, 

M^q^'qf + M^q'^q'P 

+M5 iq^q'" + qPq'^") + M^^g^P . 



(351) 

The abbreviation 5+ () 5^ () = 5^ (A;^ — m?) 5^ {q' — kf — is being used together 
with the notation 6^ {k"^ — m^) = 6 {k"^ — m?) 9 (ko), where m labels the lepton mass. 
The coefficients on the right-hand side of Eq. 1)3511) are scalar functions, which depend 
upon three invariants, i.e. m^, g'^ and (g ■ q') = {q'^ — t) /2, and can be derived by 
contracting both sides of Eq. ()35ip with the four-momenta q^ and q'^ of the initial 
real and final virtual photons, respectively, and with the tensors q^qp, q'^q'p, q'^qp and 
Qpp. These functions can be further reduced to the expressions written in terms of 
only two scalar functions, i.e. Jq and l2- After some algebra we obtain 



TT 



2 {q ■ q>) 

-,12 



In 



1 + Jl - Am?/q''^ 



I — Jl — Am?/q 



/2 



h 
h 
h 
h 
h 

h 
Jo 



2iq-q, 

TT 

2 (g ■ q') 
g'2 



TV [^0 - 2/2] 



Am? 



q'Y 



3q' 



12 



2 (q ■ q') 



-.12 



2 T ^ T 

m Iq- —h 



1 
2 



^12 



2r 1 T 



Jl 
J2 
J3 

J4 

K2 

K4 

Ke 
Lo 
Li 
L2 
L3 

Ml 

M2 



12 



2 (q ■ q') 

lo - h, 
q'^ 

2 {q ■ q'] 
3 

1 

'2(g-g') 



[^0 - 2/2] , 



J2 



^12 



1 

2 

0, 



m Iq- —I2 



-+m I0-—I2 



J2 



J2 



2 (g • q') 



'0, 



,/4 



2(?-gr LV2 

„I2 



Q'\ 2\r 3g'% 
— + /o — /2 



2 

~2(g-g' 

2 (g • g') 



TV [^0-/2], 



^/2 



— + m /o — /; 



,/2 



2r Q T 

m Iq - —I2 



q- 



rn? {q ■ q 
g'2 



q 



12 



{q . q'f [ ' 2m? '\ ' 
1 



/4 



6/, 



~ ~2 



+ 8/2 



(g ■ g') 



q'- 



{q ■ q') 



J2 



2{q-q 



n [^0 - 2/2] , 



0, 



(g • g') 
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Ma 
M4 
M5 
Me 



q' 



{q-q' 
1 



1 /, 



[q ■ q') 



1 



J2 



7:~J-2 



,/2 



q'^ , 2\ r 3g' 



,/2 



m /o — — /2 



(352) 
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APPENDIX C 

MODIFIED PARAMETRIZATION OF THE 
NONFORWARD MATRIX ELEMENT IN THE SCALAR 

TOY MODEL 

We modify the scalar parametrization J-'{z,p,r) of the nonforward matrix element 
(p — r/2| {—z/2) (-2/2) \p + r/2), introduced in Section VIII. 2, by adding the sub- 
leading O [z"^) and O (z^) terms. Namely, 

(p - r/2\ (f) (-z/2) (z/2) \p + r/2) = df3 C da e'^^''^ 



X 



z\ ,^ .. z^ 



hf {f3, a, t) + —h2f {(3, a, t) + —h^j (/?, a, t) 



(353) 



where k = jSp + ar /2. The function hf (/5, a, t) is the original double distribution or 
the so-called parent DD. The new functions, known as the daughter DDs, /12/ (/5, a, t) 
and h4^f {f3,a,t), are not independent but rather to be determined in such a way 
that the matrix element {p — r /2\(f) {—z/2) (p {z/2) \p + r /2) has only terms of the 
certain twist. In particular, the twist-2 part of the matrix element should satisfy the 
d'Alembert equation with respect to z, 

a2(p-r/2|0(-z/2)0(z/2)b + r/2),_,_2 = 0. (354) 



By imposing this condition on the right-hand side of Eq. ()353|) . one can determine 
both h2f {P, a, t) and h^j {P, a, t). Thus in the lowest order in z"^ we find 

d(3 / da e*(^-^) 

-1 J-1+\I3\ 

-k'^hf{f3,a,t)+i{k-z)h2f{f3,a,t) + 2h2f{l3,a,t)] = 0. (355) 

The second term in the brackets can be written as 

d(3 dah2f{(3,a,t)[(3— + a—\e'^''-'\ (356) 

-1 J-i+\p\ \ op da J 

and further integrating it by parts and using the boundary conditions for /12/ (/?, a, t) 
(note that, like the parent DD, both daughter distributions vanish at the boundaries 
of the support region), it takes the form 

- dp dae'^''Up— + a— + 2]h2f{P,a,t). (357) 

J-i J-i+\i3\ \ op oa I 
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Substituting Eq. (j357|) back into Eq ()355|) leads to the partial differential equation 

for h2f {(3,a,t), 

( Q Q\ 

+ /i2/(/3,a,t) = -k^hf{(3,a,t) (358) 

with the solution 

h2f{/3,a,t) = dr rhf {T/3,Ta,t) . (359) 



It is easy to see that this solution indeed satisfies Eq. ()358p since 

/oo ( d d \ 1'°° 

= -hf{(3,a,t) 

/oo 
dr rhf {t/3, ra, t) . 

(360) 

Similarly, in the next order in the expansion, the harmonic condition ()354p gives 
the partial differential equation for the daughter DD {P,a,t), 

+ = -k'h2f{f3,a,t). (361) 

The solution is 

/oo 
dr h2f {T(3,Ta,t) . (362) 

It is worth noting at this point that to evaluate the twist-2 and twist-3 contribu- 
tions to the Compton scattering amplitude using a modified parametrization ()353|) . 
the inclusion of and terms requires, in addition to Eqs. and ()317|) . four 

more integrals over z. Here we present the list of all six integrals 



27r2 (^2 - zO) (/2 + zO) 

27r2(z2_i0)' P + iO' 



2 ' 



d'z e^('-^) 



d'z e^('-^) 



27r2 (^2 _ ^0) (/2 + ^0) 



4 ' 



27r2(z2_^o)' {P + i^f ' 



,4 i{i.z) zsZpZx _ {gsph + gsxlp + Qpxh) l"^ — hlpl\ 
27r2 (^2 _ ^0) (/2 + ,0)^ 

74 i{hz) zsZpZx _ {gsph + gs\lp + gpxh) — ^hlphl'^ 

(Xi Z G TT ^ < . 

27r2(z2_i0)2 (/2 + iO)^ 



(363) 
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